#INTRO: LOADING PACKAGES ETC----
library(foreign)
library(reshape2)
library(data.table)
library(stargazer)
library(MatchIt)
#library(nonrandom)
library(ggplot2)
#library(gdata1)
library(gmodels)
library(gridExtra)
library(haven)
library(plyr)
#library(car)
#library(xlsx)
library(tidyr)
library(reshape)
library(countrycode)
#library(dGlyr)
library(readr)
library(gnm)
library(nls2)
library(nlstools)
library(stats)
library(dplyr)
library(matrixStats)
library(minpack.lm)
library(lubridate)
library(AER)
library(reshape)
library(reshape2)
library(tictoc)
library(xtable)
library(tidyverse)
library(openxlsx)
library(logr)
rm(list=ls())
cat("\014")

#READ DATA----

d1=read_csv("cses_imd.csv")#read CSES Integrated Module Dataset (IMD)
d11=read_csv("cses1.csv")#read CSES Module 1: 1996-2000
d12=read_csv("cses2.csv")#read CSES Module 2: 2001-2005
d13=read_csv("cses3.csv")#read CSES Module 1: 2006-2010
d14=read_csv("cses4.csv")#read CSES Module 1: 2011-2015
d5=read_csv("cses5.csv")#read CSES Module 5: 2016-2021
#CSES Module 1: 1996-2001
#CSES Module 2: 2001-2006
#CSES Module 3: 2006-2011
#CSES Module 4: 2011-2016
#CSES Module 5: 2016-2021


#ELECTION EXCLUSIONS----
#elections exclusion because no data on leader likability:
d1=d1[!(d1$IMD1008_MOD_2==1&d1$IMD1004!="PRT_2002"&d1$IMD1004!="AUS_2004"&d1$IMD1004!="BRA_2002"),]#in CSES#2
d1=d1[d1$IMD1004!="BELF1999",]#BELF1999: IMD3009_ no data on leader likability
d1=d1[d1$IMD1004!="BELW1999",]#BELW1999: IMD3009_ no data on leader likability

#elections exclusion for other reasons (alphabetical):
d1=d1[d1$IMD1004!="AUT_2008",]#no leader likability question
d1=d1[d1$IMD1004!="BLR_2001",]#only presidential election
d1=d1[d1$IMD1004!="BLR_2008",]#>50% NAs on vote choice
d1=d1[d1$IMD1004!="BRA_2002",]#>50% NAs on vote choice
d1=d1[d1$IMD1004!="BRA_2006",]#>50% NAs on vote choice
d1=d1[d1$IMD1004!="BRA_2010",]#>50% NAs on vote choice
d1=d1[d1$IMD1004!="BRA_2014",]#>50% NAs on vote choice
d1=d1[d1$IMD1004!="CHL_1999",]#only presidential election
d1=d1[d1$IMD1004!="CHL_2009",]#binominal system
d5=d5[d5$E1004!="CHL_2017",]#>50% NAs on vote choice
d1=d1[d1$IMD1004!="EST_2011",]#double/multiple voting ###
d1=d1[d1$IMD1004!="FRA_2012",]#only presidential election
d1=d1[d1$IMD1004!="HKG_1998",]#non free Freedom House
d1=d1[d1$IMD1004!="HKG_2000",]#non free Freedom House
d1=d1[d1$IMD1004!="HKG_2008",]#non free Freedom House
d1=d1[d1$IMD1004!="HKG_2012",]#non free Freedom House
d5=d5[d5$E1004!="HKG_2016",]#non free Freedom House
d1=d1[d1$IMD1004!="HUN_1998",]#parallel voting ###
d5=d5[d5$E1004!="HUN_2018",]#parallel voting ###
d1=d1[d1$IMD1004!="KEN_2013",]#two party elections
d1=d1[d1$IMD1004!="JPN_1996",]#parallel voting ###
d1=d1[d1$IMD1004!="JPN_2007",]#only upper house election ###
d1=d1[d1$IMD1004!="JPN_2013",]#only upper house election ###
d1=d1[d1$IMD1004!="KOR_2000",]#parallel voting ###
d1=d1[d1$IMD1004!="KOR_2008",]#parallel voting ###
d1=d1[d1$IMD1004!="KOR_2012",]#parallel voting ###
d5=d5[d5$E1004!="KOR_2016",]#parallel voting ###
d1=d1[d1$IMD1004!="LTU_1997",]#only presidential election
d5=d5[d5$E1004!="LTU_2016",]#parallel voting ###
d1=d1[d1$IMD1004!="MEX_1997",]#non free Freedom House
d1=d1[d1$IMD1004!="MEX_2012",]#non free Freedom House
d1=d1[d1$IMD1004!="MEX_2015",]#non free Freedom House
d5=d5[d5$E1004!="MNE_2016",]#non free Freedom House
d1=d1[d1$IMD1004!="PER_2000",]#non free Freedom House
d1=d1[d1$IMD1004!="PHL_2010",]#only presidential election
d1=d1[d1$IMD1004!="PHL_2016",]#non free Freedom House
d1=d1[d1$IMD1004!="ROU_1996",]#coalition coding not reliable
d1=d1[d1$IMD1004!="ROU_2009",]#only presidential election
d1=d1[d1$IMD1004!="ROU_2012",]#coalition coding not reliable
d1=d1[d1$IMD1004!="ROU_2014",]#only presidential election
d1=d1[d1$IMD1004!="RUS_1999",]#non free Freedom House
d1=d1[d1$IMD1004!="RUS_2000",]#non free Freedom House
d1=d1[d1$IMD1004!="CHE_1999",]#double/multiple voting ###
d1=d1[d1$IMD1004!="CHE_2007",]#double/multiple voting ###
d1=d1[d1$IMD1004!="CHE_2011",]#double/multiple voting ###
d1=d1[d1$IMD1004!="TWN_2008",]#only presidential election
d1=d1[d1$IMD1004!="TWN_2012",]#parallel voting ###
d5=d5[d5$E1004!="TWN_2016",]#parallel voting ###
d1=d1[d1$IMD1004!="THA_2001",]#missing vote choice data
d1=d1[d1$IMD1004!="THA_2007",]##non free Freedom House
d1=d1[d1$IMD1004!="THA_2011",]##non free Freedom House
d1=d1[d1$IMD1004!="TUR_2011",]#non free Freedom House
d1=d1[d1$IMD1004!="TUR_2015",]#non free Freedom House
d1=d1[d1$IMD1004!="UKR_1998",]#non free Freedom House
d1=d1[d1$IMD1004!="USA_1996",]#two party elections
d1=d1[d1$IMD1004!="USA_2008",]#two party elections
d1=d1[d1$IMD1004!="USA_2012",]#two party elections
d5=d5[d5$E1004!="USA_2016",]#two party elections
d1=d1[d1$IMD1004!="URY_2009",]#compulsory primaries to choose leader

#exclude for same reasons elections in second release CSES wave 5:
d5=d5[!(d5$E1004 %in% c("BRA_2018","CHL_2017","FRA_2017","HKG_2016","HUN_2018","KOR_2016","LTU_2016","MNE_2016","TUR_2018","TWN_2016","USA_2016")),]

d5$E1004[d5$E1004=="GRC_2015"]="GRC_2016"#temporarily modifying 2015s year to code all the variables with
#unique election-year, it will be modified before the analysis


#CODING ELECTIONS DETAILS:----
#compensatory voting: DEU; HUN
#double voting for local candidate but for party for share of seats: NZL
#single voting for local candidate also valid for party for share of seats: TWN
#single transferrable vote: IRL

#general:
d1$NP=NA
d1$CoalCode=0
d5$NP=NA
d5$CoalCode=0
#Argentina:#
d1$CoalCode[d1$IMD1004=="ARG_2015"]=1
d1$NP[d1$IMD1004=="ARG_2015"]=6
#Australia:#
d1$CoalCode[d1$IMD1004=="AUS_1996"]=1
d1$CoalCode[d1$IMD1004=="AUS_2004"]=1
d1$CoalCode[d1$IMD1004=="AUS_2007"]=1
d1$CoalCode[d1$IMD1004=="AUS_2013"]=1
d5$CoalCode[d5$E1004=="AUS_2019"]=1#
d1$NP[d1$IMD1004=="AUS_1996"]=4
d1$NP[d1$IMD1004=="AUS_2004"]=6#scaled up since E not coded
d1$NP[d1$IMD1004=="AUS_2007"]=4
d1$NP[d1$IMD1004=="AUS_2013"]=4#scaled up E since D not coded
d5$NP[d5$E1004=="AUS_2019"]=6#
d1$IMD3008_E[d1$IMD1004=="AUS_1996"]=NA#party likability: E,F not coded
d1$IMD3008_F[d1$IMD1004=="AUS_1996"]=NA
d1$IMD3009_E[d1$IMD1004=="AUS_1996"]=NA#leader likability: E,F not coded
d1$IMD3009_F[d1$IMD1004=="AUS_1996"]=NA
d1$IMD3008_E[d1$IMD1004=="AUS_2004"]=d1$IMD3008_F[d1$IMD1004=="AUS_2004"]#party likability: scaled up since E not coded
d1$IMD3008_F[d1$IMD1004=="AUS_2004"]=d1$IMD3008_G[d1$IMD1004=="AUS_2004"]
d1$IMD3009_E[d1$IMD1004=="AUS_2004"]=d1$IMD3009_F[d1$IMD1004=="AUS_2004"]#leader likability: scaled up since E not coded
d1$IMD3009_F[d1$IMD1004=="AUS_2004"]=d1$IMD3009_G[d1$IMD1004=="AUS_2004"]
d1$IMD3008_E[d1$IMD1004=="AUS_2007"]=NA#party likability: E,F not coded
d1$IMD3008_F[d1$IMD1004=="AUS_2007"]=NA
d1$IMD3009_E[d1$IMD1004=="AUS_2007"]=NA#leader likability: E,F not coded
d1$IMD3009_F[d1$IMD1004=="AUS_2007"]=NA
d1$IMD3008_D[d1$IMD1004=="AUS_2013"]=d1$IMD3008_E[d1$IMD1004=="AUS_2013"]#party likability: scaled up E since D not coded
d1$IMD3009_D[d1$IMD1004=="AUS_2013"]=d1$IMD3009_E[d1$IMD1004=="AUS_2013"]#leader likability: scaled up E since D not coded
d1$IMD3008_E[d1$IMD1004=="AUS_2013"]=NA#party likability: E,F not coded
d1$IMD3008_F[d1$IMD1004=="AUS_2013"]=NA
d1$IMD3009_E[d1$IMD1004=="AUS_2013"]=NA#leader likability: E,F not coded
d1$IMD3009_F[d1$IMD1004=="AUS_2013"]=NA
d1$IMD3007_E[d1$IMD1004=="AUS_1996"]=NA#LRP: E,F not coded
d1$IMD3007_F[d1$IMD1004=="AUS_1996"]=NA
d1$IMD3007_E[d1$IMD1004=="AUS_2004"]=d1$IMD3007_F[d1$IMD1004=="AUS_2004"]#LRP: scaled up since E not coded
d1$IMD3007_F[d1$IMD1004=="AUS_2004"]=d1$IMD3007_G[d1$IMD1004=="AUS_2004"]
d1$IMD3007_E[d1$IMD1004=="AUS_2007"]=NA#LRP: E,F not coded
d1$IMD3007_F[d1$IMD1004=="AUS_2007"]=NA
d1$IMD3007_D[d1$IMD1004=="AUS_2013"]=d1$IMD3007_E[d1$IMD1004=="AUS_2013"]#LRP: scaled up E since D not coded
d1$IMD3007_E[d1$IMD1004=="AUS_2013"]=NA#LRP: E,F not coded
d1$IMD3007_F[d1$IMD1004=="AUS_2013"]=NA
#Austria:#
d1$NP[d1$IMD1004=="AUT_2013"]=6
d5$NP[d5$E1004=="AUT_2017"]=6
#Bulgaria:#
d1$NP[d1$IMD1004=="BGR_2014"]=6
#Canada:#
d1$NP[d1$IMD1004=="CAN_1997"]=5
d1$NP[d1$IMD1004=="CAN_2008"]=5
d1$NP[d1$IMD1004=="CAN_2011"]=5
d1$NP[d1$IMD1004=="CAN_2015"]=5
d11$A3020_E[d11$A1004=="CAN_1997"&d11$A2019!=24]=NA#set party like Bloc Quebecois as NA outside Quebec
#NB: this choice is what was done in the Canadian NES, check for example highlights in NES 1997, 1993;
# which are supposedly collected by the EV dataset. Also, since 1997 appears in both CSES amd EV
# I checked the ratio of NAs for the Bloc Quebecois after setting it NA in Quebec for the CSES,
# and the proportion is the same. Therefore I set NA for Bloc Quebecois outside Quebec,
# and analogously I set NA party Reform inside Quebec.
d1$IMD3008_E[d1$IMD1004=="CAN_1997"]=d11$A3020_E[d11$A1004=="CAN_1997"]
d13$C3009_C[d13$C1004=="CAN_2008"&d13$C2027!=24]=NA
d1$IMD3008_C[d1$IMD1004=="CAN_2008"]=d13$C3009_C[d13$C1004=="CAN_2008"]
d14$D3011_D[d14$D1004=="CAN_2011"&d14$D2028!=24]=NA
d1$IMD3008_D[d1$IMD1004=="CAN_2011"]=d14$D3011_D[d14$D1004=="CAN_2011"]
d14$D3011_D[d14$D1004=="CAN_2015"&d14$D2028!=24]=NA
d1$IMD3008_D[d1$IMD1004=="CAN_2015"]=d14$D3011_D[d14$D1004=="CAN_2015"]
d11$A3020_B[d11$A1004=="CAN_1997"&d11$A2019==24]=NA#set party like Reform as NA inside Quebec
d1$IMD3008_B[d1$IMD1004=="CAN_1997"]=d11$A3020_B[d11$A1004=="CAN_1997"]
d11$A3021_E[d11$A1004=="CAN_1997"&d11$A2019!=24]=NA#set leader like Bloc Quebecois as NA outside Quebec
d1$IMD3009_E[d1$IMD1004=="CAN_1997"]=d11$A3021_E[d11$A1004=="CAN_1997"]
d13$C3010_C[d13$C1004=="CAN_2008"&d13$C2027!=24]=NA
d1$IMD3009_C[d1$IMD1004=="CAN_2008"]=d13$C3010_C[d13$C1004=="CAN_2008"]
d14$D3012_D[d14$D1004=="CAN_2011"&d14$D2028!=24]=NA
d1$IMD3009_D[d1$IMD1004=="CAN_2011"]=d14$D3012_D[d14$D1004=="CAN_2011"]
d14$D3012_D[d14$D1004=="CAN_2015"&d14$D2028!=24]=NA
d1$IMD3009_D[d1$IMD1004=="CAN_2015"]=d14$D3012_D[d14$D1004=="CAN_2015"]
d11$A3021_B[d11$A1004=="CAN_1997"&d11$A2019==24]=NA#set leader like Reform as NA inside Quebec
d1$IMD3009_B[d1$IMD1004=="CAN_1997"]=d11$A3021_B[d11$A1004=="CAN_1997"]
d11$A3032_E[d11$A1004=="CAN_1997"&d11$A2019!=24]=NA#set LRP Bloc Quebecois as NA outside Quebec
d1$IMD3007_E[d1$IMD1004=="CAN_1997"]=d11$A3032_E[d11$A1004=="CAN_1997"]
d13$C3011_C[d13$C1004=="CAN_2008"&d13$C2027!=24]=NA
d1$IMD3007_C[d1$IMD1004=="CAN_2008"]=d13$C3011_C[d13$C1004=="CAN_2008"]
d14$D3013_D[d14$D1004=="CAN_2011"&d14$D2028!=24]=NA
d1$IMD3007_D[d1$IMD1004=="CAN_2011"]=d14$D3013_D[d14$D1004=="CAN_2011"]
d14$D3013_D[d14$D1004=="CAN_2015"&d14$D2028!=24]=NA
d1$IMD3007_D[d1$IMD1004=="CAN_2015"]=d14$D3013_D[d14$D1004=="CAN_2015"]
d11$A3032_B[d11$A1004=="CAN_1997"&d11$A2019==24]=NA#set LRP Reform as NA inside Quebec
d1$IMD3007_B[d1$IMD1004=="CAN_1997"]=d11$A3032_B[d11$A1004=="CAN_1997"]
#Croatia:#
d1$CoalCode[d1$IMD1004=="HRV_2007"]=1
d1$NP[d1$IMD1004=="HRV_2007"]=6
d1$IMD5001_A[d1$IMD1004=="HRV_2007"]=36.6#recoding HRV_2007 percent vote (error)
d1$IMD5001_B[d1$IMD1004=="HRV_2007"]=31.3
d1$IMD5001_C[d1$IMD1004=="HRV_2007"]=6.7
#Czech Republic:#
d1$NP[d1$IMD1004=="CZE_1996"]=6
d1$NP[d1$IMD1004=="CZE_2006"]=6
d1$NP[d1$IMD1004=="CZE_2010"]=6
d1$NP[d1$IMD1004=="CZE_2013"]=6
#Denmark:#
d1$NP[d1$IMD1004=="DNK_1998"]=6
d1$NP[d1$IMD1004=="DNK_2007"]=6
#Finland:#
d1$NP[d1$IMD1004=="FIN_2007"]=6
d1$NP[d1$IMD1004=="FIN_2011"]=6
d1$NP[d1$IMD1004=="FIN_2015"]=6
#France:#
d1$NP[d1$IMD1004=="FRA_2007"]=6
#Germany:#
d1$CoalCode[d1$IMD1004=="DEU_1998"]=1
d1$CoalCode[d1$IMD1004=="DEU_2005"]=1
d1$CoalCode[d1$IMD1004=="DEU_2009"]=1
d1$CoalCode[d1$IMD1004=="DEU_2013"]=1
d5$CoalCode[d5$E1004=="DEU_2017"]=1
d1$NP[d1$IMD1004=="DEU_1998"]=5
d1$NP[d1$IMD1004=="DEU_2005"]=6
d1$NP[d1$IMD1004=="DEU_2009"]=5
d1$NP[d1$IMD1004=="DEU_2013"]=5
d5$NP[d5$E1004=="DEU_2017"]=6
d1$IMD5001_A[d1$IMD1004=="DEU_1998"]=40.9#recoding DEU_1998 percent vote
d1$IMD5001_B[d1$IMD1004=="DEU_1998"]=28.4#(they coded LH instead than PL vote share)
d1$IMD5001_C[d1$IMD1004=="DEU_1998"]=6.7
d1$IMD5001_D[d1$IMD1004=="DEU_1998"]=6.7
d1$IMD5001_E[d1$IMD1004=="DEU_1998"]=6.2
d1$IMD5001_F[d1$IMD1004=="DEU_1998"]=5.1
d11$A3020_B[d11$A1004=="DEU_1998"&d11$A2019==9]=d11$A3020_C[d11$A1004=="DEU_1998"&d11$A2019==9]#set party like CDU as CSU in Bavaria
#this is also how the TEV is coded
d1$IMD3008_B[d1$IMD1004=="DEU_1998"]=d11$A3020_B[d11$A1004=="DEU_1998"]

d1$IMD3008_C[d1$IMD1004=="DEU_1998"]=d1$IMD3008_D[d1$IMD1004=="DEU_1998"]#party likability scale up letters after removing CSU
d1$IMD3008_D[d1$IMD1004=="DEU_1998"]=d1$IMD3008_E[d1$IMD1004=="DEU_1998"]
d1$IMD3008_E[d1$IMD1004=="DEU_1998"]=d1$IMD3008_F[d1$IMD1004=="DEU_1998"]
d1$IMD3008_F[d1$IMD1004=="DEU_1998"]=NA
d13$C3009_B[d13$C1004=="DEU_2005"&d13$C2027==9]=d13$C3009_F[d13$C1004=="DEU_2005"&d13$C2027==9]#set party like CDU as CSU in Bavaria
d1$IMD3008_B[d1$IMD1004=="DEU_2005"]=d13$C3009_B[d13$C1004=="DEU_2005"]
d1$IMD3008_F[d1$IMD1004=="DEU_2005"]=d1$IMD3008_G[d1$IMD1004=="DEU_2005"]#party likability scale up letters after removing CSU
d13$C3009_A[d13$C1004=="DEU_2009"&d13$C2027==9]=d13$C3009_F[d13$C1004=="DEU_2009"&d13$C2027==9]#set party like CDU as CSU in Bavaria
d1$IMD3008_A[d1$IMD1004=="DEU_2009"]=d13$C3009_A[d13$C1004=="DEU_2009"]
d1$IMD3008_F[d1$IMD1004=="DEU_2009"]=NA#party likability scale up letters after removing CSU
d14$D3011_A[d14$D1004=="DEU_2013"&d14$D2028==2]=d14$D3011_E[d14$D1004=="DEU_2013"&d14$D2028==2]#set party like CDU as CSU in Bavaria
d1$IMD3008_A[d1$IMD1004=="DEU_2013"]=d14$D3011_A[d14$D1004=="DEU_2013"]
d1$IMD3008_E[d1$IMD1004=="DEU_2013"]=d1$IMD3008_F[d1$IMD1004=="DEU_2013"]#party likability scale up letters after removing CSU
d1$IMD3008_F[d1$IMD1004=="DEU_2013"]=NA
d5$E3017_A[d5$E1004=="DEU_2017"&d5$E2020==2]=d5$E3017_G[d5$E1004=="DEU_2017"&d5$E2020==2]#set party like CDU as CSU in Bavaria
# d5$E3017_A[d5$E1004=="DEU_2017"&d5$E2020==2]=rowMeans(d5[d5$E1004=="DEU_2017"&d5$E2020==2,names(d5) %in% c("E3017_A","E3017_G")],na.rm=T)#set party like CDU as mean CDU and CSU in Bavaria
d11$A3021_B[d11$A1004=="DEU_1998"&d11$A2019==9]=d11$A3021_C[d11$A1004=="DEU_1998"&d11$A2019==9]#set leader like CDU as CSU in Bavaria
d1$IMD3009_B[d1$IMD1004=="DEU_1998"]=d11$A3021_B[d11$A1004=="DEU_1998"]
d1$IMD3009_C[d1$IMD1004=="DEU_1998"]=d1$IMD3009_D[d1$IMD1004=="DEU_1998"]#leader likability scale up letters after removing CSU
d1$IMD3009_D[d1$IMD1004=="DEU_1998"]=d1$IMD3009_E[d1$IMD1004=="DEU_1998"]
d1$IMD3009_E[d1$IMD1004=="DEU_1998"]=d1$IMD3009_F[d1$IMD1004=="DEU_1998"]
d1$IMD3009_F[d1$IMD1004=="DEU_1998"]=NA
d13$C3010_B[d13$C1004=="DEU_2005"&d13$C2027==9]=d13$C3010_F[d13$C1004=="DEU_2005"&d13$C2027==9]#set leader like CDU as CSU in Bavaria
d1$IMD3009_B[d1$IMD1004=="DEU_2005"]=d13$C3010_B[d13$C1004=="DEU_2005"]
d1$IMD3009_F[d1$IMD1004=="DEU_2005"]=d1$IMD3009_G[d1$IMD1004=="DEU_2005"]#leader likability scale up letters after removing CSU
d13$C3010_A[d13$C1004=="DEU_2009"&d13$C2027==9]=d13$C3010_F[d13$C1004=="DEU_2009"&d13$C2027==9]#set leader like CDU as CSU in Bavaria
d1$IMD3009_A[d1$IMD1004=="DEU_2009"]=d13$C3010_A[d13$C1004=="DEU_2009"]
d1$IMD3009_F[d1$IMD1004=="DEU_2009"]=NA#party likability scale up letters after removing CSU
d14$D3012_A[d14$D1004=="DEU_2013"&d14$D2028==2]=d14$D3012_E[d14$D1004=="DEU_2013"&d14$D2028==2]
d1$IMD3009_A[d1$IMD1004=="DEU_2013"]=d14$D3012_A[d14$D1004=="DEU_2013"]
d1$IMD3009_E[d1$IMD1004=="DEU_2013"]=d1$IMD3009_F[d1$IMD1004=="DEU_2013"]#party likability scale up letters after removing CSU
d1$IMD3009_F[d1$IMD1004=="DEU_2013"]=NA
d5$E3018_A[d5$E1004=="DEU_2017"&d5$E2020==2]=d5$E3018_G[d5$E1004=="DEU_2017"&d5$E2020==2]#set leader like CDU as CSU in Bavaria
d11$A3032_B[d11$A1004=="DEU_1998"&d11$A2019==9]=d11$A3032_C[d11$A1004=="DEU_1998"&d11$A2019==9]#set LRP CDU as CSU in Bavaria
d1$IMD3007_B[d1$IMD1004=="DEU_1998"]=d11$A3032_B[d11$A1004=="DEU_1998"]
d1$IMD3007_C[d1$IMD1004=="DEU_1998"]=d1$IMD3007_D[d1$IMD1004=="DEU_1998"]#LRP scale up letters after removing CSU
d1$IMD3007_D[d1$IMD1004=="DEU_1998"]=d1$IMD3007_E[d1$IMD1004=="DEU_1998"]
d1$IMD3007_E[d1$IMD1004=="DEU_1998"]=d1$IMD3007_F[d1$IMD1004=="DEU_1998"]
d1$IMD3007_F[d1$IMD1004=="DEU_1998"]=NA
d13$C3011_B[d13$C1004=="DEU_2005"&d13$C2027==9]=d13$C3011_F[d13$C1004=="DEU_2005"&d13$C2027==9]#set LRP CDU as CSU in Bavaria
d1$IMD3007_B[d1$IMD1004=="DEU_2005"]=d13$C3011_B[d13$C1004=="DEU_2005"]
d1$IMD3007_F[d1$IMD1004=="DEU_2005"]=d1$IMD3007_G[d1$IMD1004=="DEU_2005"]#LRP scale up letters after removing CSU
d13$C3011_A[d13$C1004=="DEU_2009"&d13$C2027==9]=d13$C3011_F[d13$C1004=="DEU_2009"&d13$C2027==9]#set LRP CDU as CSU in Bavaria
d1$IMD3007_A[d1$IMD1004=="DEU_2009"]=d13$C3011_A[d13$C1004=="DEU_2009"]
d1$IMD3007_F[d1$IMD1004=="DEU_2009"]=NA#LRP scale up letters after removing CSU
d14$D3013_A[d14$D1004=="DEU_2013"&d14$D2028==2]=d14$D3013_E[d14$D1004=="DEU_2013"&d14$D2028==2]#set LRP CDU as CSU in Bavaria
d1$IMD3007_A[d1$IMD1004=="DEU_2013"]=d14$D3013_A[d14$D1004=="DEU_2013"]
d1$IMD3007_E[d1$IMD1004=="DEU_2013"]=d1$IMD3007_F[d1$IMD1004=="DEU_2013"]#LRP scale up letters after removing CSU
d1$IMD3007_F[d1$IMD1004=="DEU_2013"]=NA
d5$E3019_A[d5$E1004=="DEU_2017"&d5$E2020==2]=d5$E3019_G[d5$E1004=="DEU_2017"&d5$E2020==2]#set LRP CDU as CSU in Bavaria

#exclude East Germany in order to have West Germany as congruent electorate (and exclude Eastern Europe):
d11_temp=d11[d11$A1004=="DEU_1998",c("A1005","A2019")]#keep only DEU_1998, RespID and region
d1=merge(d1, d11_temp, by.x=c("IMD1005"), by.y=c("A1005"), all.x=T)
d13_temp=d13[d13$C1004=="DEU_2005",c("C1005","C2027")]#keep only DEU_2005, RespID and region
d1=merge(d1, d13_temp, by.x=c("IMD1005"), by.y=c("C1005"), all.x=T)
colnames(d1)[colnames(d1)=="C2027"]="C2027_"
d13_temp=d13[d13$C1004=="DEU_2009",c("C1005","C2027")]#keep only DEU_1998, RespID and region
d1=merge(d1, d13_temp, by.x=c("IMD1005"), by.y=c("C1005"), all.x=T)
d14_temp=d14[d14$D1004=="DEU_2013",c("D1005","D2028")]#keep only DEU_1998, RespID and region
d1=merge(d1, d14_temp, by.x=c("IMD1005"), by.y=c("D1005"), all.x=T)
d1=d1[!(d1$A2019 %in% c(11,12,13,14,15,16,17)),]#filter out East Germany regions
d1=d1[!(d1$C2027_ %in% c(11,12,13,14,15,16,17)),]
d1=d1[!(d1$C2027 %in% c(11,12,13,14,15,16,17)),]
d1=d1[!(d1$D2028 %in% c(3,4,8,13,14,16)),]
d1=d1[,-which(names(d1) %in% c("A2019","C2027_","C2027","D2028"))]#remove working variables and datasets:
rm(d11_temp,d13_temp,d14_temp)
d5=d5[!(d5$E1004=="DEU_2017" & d5$E2020 %in% c(3,4,8,13,14,16)),]#remove East Germany in CSES5

#Great Britain:#
d1$NP[d1$IMD1004=="GBR_1997"]=5
d1$NP[d1$IMD1004=="GBR_2015"]=6
#Greece:# (2009 different law, but almost same disproportional mechanism, kept the same)
d1$NP[d1$IMD1004=="GRC_2009"]=6
d1$NP[d1$IMD1004=="GRC_2012"]=6
d1$NP[d1$IMD1004=="GRC_2015"]=6
d5$NP[d5$E1004=="GRC_2016"]=6
#Iceland:#
d1$NP[d1$IMD1004=="ISL_1999"]=5
d1$NP[d1$IMD1004=="ISL_2007"]=6
d1$NP[d1$IMD1004=="ISL_2009"]=6
d1$NP[d1$IMD1004=="ISL_2013"]=6
d5$NP[d5$E1004=="ISL_2016"]=6#
d5$NP[d5$E1004=="ISL_2017"]=6#
#Ireland:#
d1$NP[d1$IMD1004=="IRL_2007"]=6
d1$NP[d1$IMD1004=="IRL_2011"]=5
d5$NP[d5$E1004=="IRL_2016"]=6
d1$IMD3008_E[d1$IMD1004=="IRL_2011"]=d1$IMD3008_F[d1$IMD1004=="IRL_2011"]#party likability: E not coded, substituted with F
d1$IMD3008_F[d1$IMD1004=="IRL_2011"]=NA
d1$IMD3009_E[d1$IMD1004=="IRL_2011"]=d1$IMD3009_F[d1$IMD1004=="IRL_2011"]#leader likability: E not coded, substituted with F
d1$IMD3009_F[d1$IMD1004=="IRL_2011"]=NA
d1$IMD3007_E[d1$IMD1004=="IRL_2011"]=d1$IMD3007_F[d1$IMD1004=="IRL_2011"]#LRP: E not coded, substituted with F
d1$IMD3007_F[d1$IMD1004=="IRL_2011"]=NA
#Israel:#
d1$NP[d1$IMD1004=="ISR_1996"]=5
d1$NP[d1$IMD1004=="ISR_2006"]=6
d1$NP[d1$IMD1004=="ISR_2013"]=6
d1$IMD3008_F[d1$IMD1004=="ISR_2013"]=d1$IMD3008_G[d1$IMD1004=="ISR_2013"]#party likability: F not coded, substituted with G
d1$IMD3009_F[d1$IMD1004=="ISR_2013"]=d1$IMD3009_G[d1$IMD1004=="ISR_2013"]#leader likability: F not coded, substituted with G
d1$IMD3007_F[d1$IMD1004=="ISR_2013"]=d1$IMD3007_G[d1$IMD1004=="ISR_2013"]#party likability: F not coded, substituted with G
#Italy:#
d5$CoalCode[d5$E1004=="ITA_2018"]=1
d5$NP[d5$E1004=="ITA_2018"]=6
d5$E5002_B[d5$E1004=="ITA_2018"]=17.8#recoding ITA_2018 percent seats (error)
#Latvia:#
d1$CoalCode[d1$IMD1004=="LVA_2010"]=1
d1$CoalCode[d1$IMD1004=="LVA_2011"]=1
d1$CoalCode[d1$IMD1004=="LVA_2014"]=1
d1$NP[d1$IMD1004=="LVA_2010"]=6
d1$NP[d1$IMD1004=="LVA_2011"]=6
d1$NP[d1$IMD1004=="LVA_2014"]=6
#Mexico:#
d1$CoalCode[d1$IMD1004=="MEX_2000"]=1
d1$CoalCode[d1$IMD1004=="MEX_2006"]=1
d1$CoalCode[d1$IMD1004=="MEX_2009"]=1
d1$NP[d1$IMD1004=="MEX_2000"]=3
d1$NP[d1$IMD1004=="MEX_2006"]=3
d1$NP[d1$IMD1004=="MEX_2009"]=3
d1$IMD5002_A[d1$IMD1004=="MEX_2000"]=44.6#recoding MEX_2000 percent seats (error)
d1$IMD5002_B[d1$IMD1004=="MEX_2000"]=42.2
d1$IMD5002_C[d1$IMD1004=="MEX_2000"]=13.2
d1$IMD5002_A[d1$IMD1004=="MEX_2006"]=41.2#recoding MEX_2006 percent seats (error)
d1$IMD5002_B[d1$IMD1004=="MEX_2006"]=31.6
d1$IMD5002_C[d1$IMD1004=="MEX_2006"]=24.6
d1$IMD5001_A[d1$IMD1004=="MEX_2009"]=36.9#recoding MEX_2009 percent vote (error)
d1$IMD5001_B[d1$IMD1004=="MEX_2009"]=28.0
d1$IMD5001_C[d1$IMD1004=="MEX_2009"]=12.2
d1$IMD5002_A[d1$IMD1004=="MEX_2009"]=47.4#recoding MEX_2006 percent seats (error)
d1$IMD5002_B[d1$IMD1004=="MEX_2009"]=28.6
d1$IMD5002_C[d1$IMD1004=="MEX_2009"]=14.2
d1$IMD3008_D[d1$IMD1004=="MEX_2000"]=NA#party likability MEX: D,E,F not coded
d1$IMD3008_E[d1$IMD1004=="MEX_2000"]=NA
d1$IMD3008_F[d1$IMD1004=="MEX_2000"]=NA
d1$IMD3008_D[d1$IMD1004=="MEX_2006"]=NA
d1$IMD3008_E[d1$IMD1004=="MEX_2006"]=NA
d1$IMD3008_F[d1$IMD1004=="MEX_2006"]=NA
d1$IMD3008_D[d1$IMD1004=="MEX_2009"]=NA
d1$IMD3008_E[d1$IMD1004=="MEX_2009"]=NA
d1$IMD3008_F[d1$IMD1004=="MEX_2009"]=NA
d1$IMD3009_D[d1$IMD1004=="MEX_2000"]=NA#leader likability MEX: D,E,F not coded
d1$IMD3009_E[d1$IMD1004=="MEX_2000"]=NA
d1$IMD3009_F[d1$IMD1004=="MEX_2000"]=NA
d1$IMD3009_D[d1$IMD1004=="MEX_2006"]=NA
d1$IMD3009_E[d1$IMD1004=="MEX_2006"]=NA
d1$IMD3009_F[d1$IMD1004=="MEX_2006"]=NA
d1$IMD3009_D[d1$IMD1004=="MEX_2009"]=NA
d1$IMD3009_E[d1$IMD1004=="MEX_2009"]=NA
d1$IMD3009_F[d1$IMD1004=="MEX_2009"]=NA
d1$IMD3007_D[d1$IMD1004=="MEX_2000"]=NA#LRP MEX: D,E,F not coded
d1$IMD3007_E[d1$IMD1004=="MEX_2000"]=NA
d1$IMD3007_F[d1$IMD1004=="MEX_2000"]=NA
d1$IMD3007_D[d1$IMD1004=="MEX_2006"]=NA
d1$IMD3007_E[d1$IMD1004=="MEX_2006"]=NA
d1$IMD3007_F[d1$IMD1004=="MEX_2006"]=NA
d1$IMD3007_D[d1$IMD1004=="MEX_2009"]=NA
d1$IMD3007_E[d1$IMD1004=="MEX_2009"]=NA
d1$IMD3007_F[d1$IMD1004=="MEX_2009"]=NA
#Montenegro:#
d1$CoalCode[d1$IMD1004=="MNE_2012"]=1
d1$NP[d1$IMD1004=="MNE_2012"]=6
#Netherlands:#
d1$NP[d1$IMD1004=="NLD_1998"]=6
d1$NP[d1$IMD1004=="NLD_2006"]=6
d1$NP[d1$IMD1004=="NLD_2010"]=6
#New Zealand:#
d1$NP[d1$IMD1004=="NZL_1996"]=6
d1$NP[d1$IMD1004=="NZL_2008"]=6
d1$NP[d1$IMD1004=="NZL_2011"]=6
d1$NP[d1$IMD1004=="NZL_2014"]=6
d5$NP[d5$E1004=="NZL_2017"]=6#
d1$IMD5001_A[d1$IMD1004=="NZL_1996"]=33.8#recoding NZL_1996 percent vote
d1$IMD5001_B[d1$IMD1004=="NZL_1996"]=28.2#(they coded LH instead than PL percent vote)
d1$IMD5001_C[d1$IMD1004=="NZL_1996"]=13.3
d1$IMD5001_D[d1$IMD1004=="NZL_1996"]=10.1
d1$IMD5001_E[d1$IMD1004=="NZL_1996"]=6.1
d1$IMD5001_F[d1$IMD1004=="NZL_1996"]=4.3
#Norway:#
d1$NP[d1$IMD1004=="NOR_1997"]=6
d1$NP[d1$IMD1004=="NOR_2005"]=6
d1$NP[d1$IMD1004=="NOR_2009"]=6
d1$NP[d1$IMD1004=="NOR_2013"]=6
d5$NP[d5$E1004=="NOR_2017"]=6#
#Peru:#
d1$CoalCode[d1$IMD1004=="PER_2011"]=1
d1$CoalCode[d1$IMD1004=="PER_2016"]=1
d1$NP[d1$IMD1004=="PER_2001"]=6
d1$NP[d1$IMD1004=="PER_2011"]=5
d1$NP[d1$IMD1004=="PER_2016"]=6
d1$IMD5002_B[d1$IMD1004=="PER_2001"]=23.3#recoding PER_2011 percent seats (error)
d1$IMD3008_F[d1$IMD1004=="PER_2011"]=NA#party likability: F leader not coded
d1$IMD3009_F[d1$IMD1004=="PER_2011"]=NA#leader likability: F leader not coded
d1$IMD3007_F[d1$IMD1004=="PER_2011"]=NA#LRP: F leader not coded
#Poland:#
d1$NP[d1$IMD1004=="POL_1997"]=6
d1$NP[d1$IMD1004=="POL_2005"]=6
d1$NP[d1$IMD1004=="POL_2007"]=4
d1$NP[d1$IMD1004=="POL_2011"]=5
d1$IMD3008_F[d1$IMD1004=="POL_2005"]=d1$IMD3008_H[d1$IMD1004=="POL_2005"]#party likability: F not coded, substituted with H
d1$IMD3009_F[d1$IMD1004=="POL_2005"]=d1$IMD3009_H[d1$IMD1004=="POL_2005"]#leader likability: F not coded, substituted with H
d1$IMD3007_F[d1$IMD1004=="POL_2005"]=d1$IMD3007_H[d1$IMD1004=="POL_2005"]#LRP: F not coded, substituted with H
#Portugal:#
d1$NP[d1$IMD1004=="PRT_2002"]=6
d1$NP[d1$IMD1004=="PRT_2009"]=5
d1$NP[d1$IMD1004=="PRT_2015"]=5
d1$IMD5001_A[d1$IMD1004=="PRT_2002"]=40.2#recoding percent vote (error)
d1$IMD5001_B[d1$IMD1004=="PRT_2002"]=37.8
d1$IMD3008_E[d1$IMD1004=="PRT_2015"]=d1$IMD3008_F[d1$IMD1004=="PRT_2015"]#party likability
d1$IMD3008_F[d1$IMD1004=="PRT_2015"]=NA#scale up letters since _E not coded
d1$IMD3009_E[d1$IMD1004=="PRT_2015"]=d1$IMD3009_F[d1$IMD1004=="PRT_2015"]#leader likability
d1$IMD3009_F[d1$IMD1004=="PRT_2015"]=NA#scale up letters since _E not coded
d1$IMD3007_E[d1$IMD1004=="PRT_2015"]=d1$IMD3007_F[d1$IMD1004=="PRT_2015"]#LRP
d1$IMD3007_F[d1$IMD1004=="PRT_2015"]=NA#scale up letters since _E not coded
#Serbia:#
d1$NP[d1$IMD1004=="SRB_2012"]=6
#Slovakia:#
d1$NP[d1$IMD1004=="SVK_2010"]=6
d1$NP[d1$IMD1004=="SVK_2016"]=6
#Slovenia:#
d1$NP[d1$IMD1004=="SVN_1996"]=6
d1$NP[d1$IMD1004=="SVN_2008"]=6
d1$NP[d1$IMD1004=="SVN_2011"]=6
d1$IMD5002_A[d1$IMD1004=="SVN_2011"]=31.1#recoding SVN_2011 percent seats (error)
d1$IMD5002_B[d1$IMD1004=="SVN_2011"]=28.9
d1$IMD5002_C[d1$IMD1004=="SVN_2011"]=11.1
d1$IMD5002_D[d1$IMD1004=="SVN_2011"]=8.9
d1$IMD5002_E[d1$IMD1004=="SVN_2011"]=6.7
d1$IMD5002_F[d1$IMD1004=="SVN_2011"]=6.7
d1$IMD3008_F[d1$IMD1004=="SVN_1996"]=d1$IMD3008_G[d1$IMD1004=="SVN_1996"]#party likability: F not coded, substituted with G
d1$IMD3009_F[d1$IMD1004=="SVN_1996"]=d1$IMD3009_G[d1$IMD1004=="SVN_1996"]#leader likability: F not coded, substituted with G
d1$IMD3007_F[d1$IMD1004=="SVN_1996"]=d1$IMD3007_G[d1$IMD1004=="SVN_1996"]#LRP: F not coded, substituted with G
#South Africa:#
d1$NP[d1$IMD1004=="ZAF_2009"]=6
d1$NP[d1$IMD1004=="ZAF_2014"]=6
#Spain:#
d1$NP[d1$IMD1004=="ESP_1996"]=5#no data for _F
d1$NP[d1$IMD1004=="ESP_2000"]=5#no data for _F
d1$NP[d1$IMD1004=="ESP_2008"]=6
d11$A3020_D[d11$A1004=="ESP_1996"&d11$A2019!=5]=NA#set party like Convergence and Union as NA outside Cataluna
d1$IMD3008_D[d1$IMD1004=="ESP_1996"]=d11$A3020_D[d11$A1004=="ESP_1996"]
d11$A3020_D[d11$A1004=="ESP_2000"&d11$A2019!=5]=NA
d1$IMD3008_D[d1$IMD1004=="ESP_2000"]=d11$A3020_D[d11$A1004=="ESP_2000"]
d13$C3009_C[d13$C1004=="ESP_2008"&d13$C2027!=5]=NA
d1$IMD3008_C[d1$IMD1004=="ESP_2008"]=d13$C3009_C[d13$C1004=="ESP_2008"]
d11$A3020_E[d11$A1004=="ESP_1996"&!(d11$A2019%in%c(15,16))]=NA#set party like Basque Nationalist Party as NA outside Navarra and Pais Vasco
d1$IMD3008_E[d1$IMD1004=="ESP_1996"]=d11$A3020_E[d11$A1004=="ESP_1996"]
d11$A3020_E[d11$A1004=="ESP_2000"&!(d11$A2019%in%c(15,16))]=NA
d1$IMD3008_E[d1$IMD1004=="ESP_2000"]=d11$A3020_E[d11$A1004=="ESP_2000"]
d13$C3009_D[d13$C1004=="ESP_2008"&!(d13$C2027%in%c(15,16))]=NA
d1$IMD3008_D[d1$IMD1004=="ESP_2008"]=d13$C3009_D[d13$C1004=="ESP_2008"]
d13$C3009_E[d13$C1004=="ESP_2008"&!(d13$C2027%in%c(5,17))]=NA#set party like Republican Left of Catalonia as NA outside Cataluna and Comunidad Valenciana
d1$IMD3008_E[d1$IMD1004=="ESP_2008"]=d13$C3009_E[d13$C1004=="ESP_2008"]
d11$A3021_D[d11$A1004=="ESP_1996"&d11$A2019!=5]=NA#set leader like Convergence and Union as NA outside Cataluna
d1$IMD3009_D[d1$IMD1004=="ESP_1996"]=d11$A3021_D[d11$A1004=="ESP_1996"]
d11$A3021_D[d11$A1004=="ESP_2000"&d11$A2019!=5]=NA
d1$IMD3009_D[d1$IMD1004=="ESP_2000"]=d11$A3021_D[d11$A1004=="ESP_2000"]
d13$C3010_C[d13$C1004=="ESP_2008"&d13$C2027!=5]=NA
d1$IMD3009_C[d1$IMD1004=="ESP_2008"]=d13$C3010_C[d13$C1004=="ESP_2008"]
d11$A3021_E[d11$A1004=="ESP_1996"&!(d11$A2019%in%c(15,16))]=NA#set leader like Basque Nationalist Party as NA outside Navarra and Pais Vasco
d1$IMD3009_E[d1$IMD1004=="ESP_1996"]=d11$A3021_E[d11$A1004=="ESP_1996"]
d11$A3021_E[d11$A1004=="ESP_2000"&!(d11$A2019%in%c(15,16))]=NA
d1$IMD3009_E[d1$IMD1004=="ESP_2000"]=d11$A3021_E[d11$A1004=="ESP_2000"]
d13$C3010_D[d13$C1004=="ESP_2008"&!(d13$C2027%in%c(15,16))]=NA
d1$IMD3009_D[d1$IMD1004=="ESP_2008"]=d13$C3010_D[d13$C1004=="ESP_2008"]
d13$C3010_E[d13$C1004=="ESP_2008"&!(d13$C2027%in%c(5,17))]=NA#set leader like Republican Left of Catalonia as NA outside Cataluna and Comunidad Valenciana
d1$IMD3009_E[d1$IMD1004=="ESP_2008"]=d13$C3010_E[d13$C1004=="ESP_2008"]
d11$A3032_D[d11$A1004=="ESP_1996"&d11$A2019!=5]=NA#set LRP Convergence and Union as NA outside Cataluna
d1$IMD3007_D[d1$IMD1004=="ESP_1996"]=d11$A3032_D[d11$A1004=="ESP_1996"]
d11$A3032_D[d11$A1004=="ESP_2000"&d11$A2019!=5]=NA
d1$IMD3007_D[d1$IMD1004=="ESP_2000"]=d11$A3032_D[d11$A1004=="ESP_2000"]
d13$C3011_C[d13$C1004=="ESP_2008"&d13$C2027!=5]=NA
d1$IMD3007_C[d1$IMD1004=="ESP_2008"]=d13$C3011_C[d13$C1004=="ESP_2008"]
d11$A3032_E[d11$A1004=="ESP_1996"&!(d11$A2019%in%c(15,16))]=NA#set LRP Basque Nationalist Party as NA outside Navarra and Pais Vasco
d1$IMD3007_E[d1$IMD1004=="ESP_1996"]=d11$A3032_E[d11$A1004=="ESP_1996"]
d11$A3032_E[d11$A1004=="ESP_2000"&!(d11$A2019%in%c(15,16))]=NA
d1$IMD3007_E[d1$IMD1004=="ESP_2000"]=d11$A3032_E[d11$A1004=="ESP_2000"]
d13$C3011_D[d13$C1004=="ESP_2008"&!(d13$C2027%in%c(15,16))]=NA
d1$IMD3007_D[d1$IMD1004=="ESP_2008"]=d13$C3011_D[d13$C1004=="ESP_2008"]
d13$C3011_E[d13$C1004=="ESP_2008"&!(d13$C2027%in%c(5,17))]=NA#set LRP Republican Left of Catalonia as NA outside Cataluna and Comunidad Valenciana
d1$IMD3007_E[d1$IMD1004=="ESP_2008"]=d13$C3011_E[d13$C1004=="ESP_2008"]
#Sweden:#
d1$NP[d1$IMD1004=="SWE_1998"]=6
d1$NP[d1$IMD1004=="SWE_2006"]=6
d1$NP[d1$IMD1004=="SWE_2014"]=6
d1$IMD5002_A[d1$IMD1004=="SWE_1998"]=37.5#recoding SWE_1998 percent seats (error)
#Taiwan:#
d1$NP[d1$IMD1004=="TWN_1996"]=3


#D1----
#Election variables:
colnames(d1)[colnames(d1)=="IMD1001"]="Ed"#Election Dataset
d1$Ed[d1$IMD1008_MOD_1==1]="CSES-MODULE-1"
d1$Ed[d1$IMD1008_MOD_2==1]="CSES-MODULE-2"#Note: PRT_2002 is assigned to both modules 1 and 2, I'm assigning it to only module 2 bc of year
d1$Ed[d1$IMD1008_MOD_3==1]="CSES-MODULE-3"
d1$Ed[d1$IMD1008_MOD_4==1]="CSES-MODULE-4"
colnames(d1)[colnames(d1)=="IMD1004"]="Es"#Election Study
d1$Ec=substr(d1$Es, 1, 3)#Election Country
colnames(d1)[colnames(d1)=="IMD1008_YEAR"]="Ey"#Election Year
colnames(d1)[colnames(d1)=="IMD1005"]="Rid"#Respondent id
colnames(d1)[colnames(d1)=="IMD2001_1"]="Rage"#Respondent Age
d1$Rage[d1$Rage>115]=NA
d1$Rwd=1#Respondent Weight Demographic:
d1$Rwd[d1$Es=="CAN_1997"]=d1$IMD1010_1[d1$Es=="CAN_1997"]
d1$Rwd[d1$Es=="DNK_1998"]=d1$IMD1010_2[d1$Es=="DNK_1998"]
d1$Rwd[d1$Es=="ESP_1996"]=d1$IMD1010_2[d1$Es=="ESP_1996"]
d1$Rwd[d1$Es=="ESP_2000"]=d1$IMD1010_2[d1$Es=="ESP_2000"]
d1$Rwd[d1$Es=="ESP_2000"&d1$Rwd==0]=1#missing weight put to neutral weight of 1
d1$Rwd[d1$Es=="GBR_1997"]=d1$IMD1010_1[d1$Es=="GBR_1997"]
d1$Rwd[d1$Es=="NZL_1996"]=d1$IMD1010_1[d1$Es=="NZL_1996"]
d1$Rwd[d1$Es=="CAN_2008"]=d1$IMD1010_1[d1$Es=="CAN_2008"]
d1$Rwd[d1$Es=="DEU_2005"]=d1$IMD1010_2[d1$Es=="DEU_2005"]
d1$Rwd[d1$Es=="DNK_2007"]=d1$IMD1010_1[d1$Es=="DNK_2007"]
d1$Rwd[d1$Es=="ESP_2008"]=d1$IMD1010_2[d1$Es=="ESP_2008"]
d1$Rwd[d1$Es=="IRL_2007"]=d1$IMD1010_2[d1$Es=="IRL_2007"]
d1$Rwd[d1$Es=="NLD_2006"]=d1$IMD1010_2[d1$Es=="NLD_2006"]
d1$Rwd[d1$Es=="NLD_2010"]=d1$IMD1010_2[d1$Es=="NLD_2010"]
d1$Rwd[d1$Es=="NZL_2008"]=d1$IMD1010_2[d1$Es=="NZL_2008"]
d1$Rwd[d1$Es=="AUS_2013"]=d1$IMD1010_2[d1$Es=="AUS_2013"]
d1$Rwd[d1$Es=="AUT_2013"]=d1$IMD1010_2[d1$Es=="AUT_2013"]
d1$Rwd[d1$Es=="CAN_2011"]=d1$IMD1010_1[d1$Es=="CAN_2011"]
d1$Rwd[d1$Es=="CAN_2015"]=d1$IMD1010_1[d1$Es=="CAN_2015"]
d1$Rwd[d1$Es=="FIN_2015"]=d1$IMD1010_2[d1$Es=="FIN_2015"]#also political
d1$Rwd[d1$Es=="GBR_2015"]=d1$IMD1010_2[d1$Es=="GBR_2015"]
d1$Rwd[d1$Es=="GRC_2015"]=d1$IMD1010_2[d1$Es=="GRC_2015"]#also political
d1$Rwd[d1$Es=="NOR_2013"]=d1$IMD1010_2[d1$Es=="NOR_2013"]
d1$Rwd[d1$Es=="NZL_2011"]=d1$IMD1010_2[d1$Es=="NZL_2011"]
d1$Rwd[d1$Es=="NZL_2014"]=d1$IMD1010_2[d1$Es=="NZL_2014"]
d1$Rwd[d1$Es=="PRT_2015"]=d1$IMD1010_2[d1$Es=="PRT_2015"]
#Vote and party/leader features:
colnames(d1)[colnames(d1)=="IMD3001_LH"]="Vt"#Vote Turnout (0=non voted; 1=voted) TURNOUT - CURRENT LOWER HOUSE ELECTION
d1$Vt[d1$Vt>1]=NA#coding missing as NAs
colnames(d1)[colnames(d1)=="IMD3002_LH_PL"]="Vn"#Vote Numerical CURRENT LOWER HOUSE ELECTION: VOTE CHOICE - PARTY LIST
d1$Vn[d1$Es=="AUS_1996"]=d1$IMD3002_LH_DC[d1$Es=="AUS_1996"]#assigning LH_DC to Vn for relevant elections
d1$Vn[d1$Es=="AUS_2004"]=d1$IMD3002_LH_DC[d1$Es=="AUS_2004"]
d1$Vn[d1$Es=="AUS_2007"]=d1$IMD3002_LH_DC[d1$Es=="AUS_2007"]
d1$Vn[d1$Es=="AUS_2013"]=d1$IMD3002_LH_DC[d1$Es=="AUS_2013"]
d1$Vn[d1$Es=="CAN_1997"]=d1$IMD3002_LH_DC[d1$Es=="CAN_1997"]
d1$Vn[d1$Es=="CAN_2008"]=d1$IMD3002_LH_DC[d1$Es=="CAN_2008"]
d1$Vn[d1$Es=="CAN_2011"]=d1$IMD3002_LH_DC[d1$Es=="CAN_2011"]
d1$Vn[d1$Es=="CAN_2015"]=d1$IMD3002_LH_DC[d1$Es=="CAN_2015"]
d1$Vn[d1$Es=="FRA_2007"]=d1$IMD3002_LH_DC[d1$Es=="FRA_2007"]
d1$Vn[d1$Es=="GBR_1997"]=d1$IMD3002_LH_DC[d1$Es=="GBR_1997"]
d1$Vn[d1$Es=="GBR_2015"]=d1$IMD3002_LH_DC[d1$Es=="GBR_2015"]
d1$Vn[d1$Es=="IRL_2007"]=d1$IMD3002_LH_DC[d1$Es=="IRL_2007"]
d1$Vn[d1$Es=="IRL_2011"]=d1$IMD3002_LH_DC[d1$Es=="IRL_2011"]
d1$Vn[d1$Es=="MEX_2000"]=d1$IMD3002_LH_DC[d1$Es=="MEX_2000"]
d1$Vn[d1$Es=="MEX_2006"]=d1$IMD3002_LH_DC[d1$Es=="MEX_2006"]
d1$Vn[d1$Es=="MEX_2009"]=d1$IMD3002_LH_DC[d1$Es=="MEX_2009"]
d1$Vn[d1$Es=="TWN_1996"]=d1$IMD3002_LH_DC[d1$Es=="TWN_1996"]

d1$Va=NA#Vote Alphabetical: code vote from numerical (000001-9000000) to alphabetical modified (1:6):
#d1 assign alphabetical modified codes instead than numerical codes----
d1$Va[d1$Es=="ARG_2015"&d1$Vn==0320001]=1#
d1$Va[d1$Es=="ARG_2015"&d1$Vn==0320002]=2
d1$Va[d1$Es=="ARG_2015"&d1$Vn==0320003]=3
d1$Va[d1$Es=="ARG_2015"&d1$Vn==0320004]=4
d1$Va[d1$Es=="ARG_2015"&d1$Vn==0320005]=5
d1$Va[d1$Es=="ARG_2015"&d1$Vn==0320006]=6

d1$Va[d1$Es=="AUS_1996"&d1$Vn==0360003]=1#
d1$Va[d1$Es=="AUS_1996"&d1$Vn==0360001]=2
d1$Va[d1$Es=="AUS_1996"&d1$Vn==0360002]=3
d1$Va[d1$Es=="AUS_1996"&d1$Vn==0360004]=4

d1$Va[d1$Es=="AUS_2004"&d1$Vn==0360001]=1#scaled up since E not coded
d1$Va[d1$Es=="AUS_2004"&d1$Vn==0360003]=2
d1$Va[d1$Es=="AUS_2004"&d1$Vn==0360005]=3
d1$Va[d1$Es=="AUS_2004"&d1$Vn==0360002]=4
d1$Va[d1$Es=="AUS_2004"&d1$Vn==0360004]=5
d1$Va[d1$Es=="AUS_2004"&d1$Vn==0360008]=6

d1$Va[d1$Es=="AUS_2007"&d1$Vn==0360003]=1#
d1$Va[d1$Es=="AUS_2007"&d1$Vn==0360001]=2
d1$Va[d1$Es=="AUS_2007"&d1$Vn==0360002]=3
d1$Va[d1$Es=="AUS_2007"&d1$Vn==0360005]=4

d1$Va[d1$Es=="AUS_2013"&d1$Vn==0360001]=1#scaled up E since D not coded
d1$Va[d1$Es=="AUS_2013"&d1$Vn==0360003]=2
d1$Va[d1$Es=="AUS_2013"&d1$Vn==0360005]=3
d1$Va[d1$Es=="AUS_2013"&d1$Vn==0360002]=4

d1$Va[d1$Es=="AUT_2013"&d1$Vn==0400001]=1#
d1$Va[d1$Es=="AUT_2013"&d1$Vn==0400002]=2
d1$Va[d1$Es=="AUT_2013"&d1$Vn==0400003]=3
d1$Va[d1$Es=="AUT_2013"&d1$Vn==0400005]=4
d1$Va[d1$Es=="AUT_2013"&d1$Vn==0400004]=5
d1$Va[d1$Es=="AUT_2013"&d1$Vn==0400010]=6

d1$Va[d1$Es=="BGR_2014"&d1$Vn==1000001]=1#
d1$Va[d1$Es=="BGR_2014"&d1$Vn==1000020]=2
d1$Va[d1$Es=="BGR_2014"&d1$Vn==1000003]=3
d1$Va[d1$Es=="BGR_2014"&d1$Vn==1000004]=4
d1$Va[d1$Es=="BGR_2014"&d1$Vn==1000005]=5
d1$Va[d1$Es=="BGR_2014"&d1$Vn==1000006]=6

d1$Va[d1$Es=="CAN_1997"&d1$Vn==1240001]=1#
d1$Va[d1$Es=="CAN_1997"&d1$Vn==1240007]=2
d1$Va[d1$Es=="CAN_1997"&d1$Vn==1240006]=3
d1$Va[d1$Es=="CAN_1997"&d1$Vn==1240003]=4
d1$Va[d1$Es=="CAN_1997"&d1$Vn==1240004]=5

d1$Va[d1$Es=="CAN_2008"&d1$Vn==1240002]=1#
d1$Va[d1$Es=="CAN_2008"&d1$Vn==1240001]=2
d1$Va[d1$Es=="CAN_2008"&d1$Vn==1240004]=3
d1$Va[d1$Es=="CAN_2008"&d1$Vn==1240003]=4
d1$Va[d1$Es=="CAN_2008"&d1$Vn==1240005]=5

d1$Va[d1$Es=="CAN_2011"&d1$Vn==1240002]=1#
d1$Va[d1$Es=="CAN_2011"&d1$Vn==1240003]=2
d1$Va[d1$Es=="CAN_2011"&d1$Vn==1240001]=3
d1$Va[d1$Es=="CAN_2011"&d1$Vn==1240004]=4
d1$Va[d1$Es=="CAN_2011"&d1$Vn==1240005]=5

d1$Va[d1$Es=="CAN_2015"&d1$Vn==1240001]=1#
d1$Va[d1$Es=="CAN_2015"&d1$Vn==1240002]=2
d1$Va[d1$Es=="CAN_2015"&d1$Vn==1240003]=3
d1$Va[d1$Es=="CAN_2015"&d1$Vn==1240004]=4
d1$Va[d1$Es=="CAN_2015"&d1$Vn==1240005]=5

d1$Va[d1$Es=="CZE_1996"&d1$Vn==2030002]=1#
d1$Va[d1$Es=="CZE_1996"&d1$Vn==2030001]=2
d1$Va[d1$Es=="CZE_1996"&d1$Vn==2030004]=3
d1$Va[d1$Es=="CZE_1996"&d1$Vn==2030003]=4
d1$Va[d1$Es=="CZE_1996"&d1$Vn==2030018]=5
d1$Va[d1$Es=="CZE_1996"&d1$Vn==2030011]=6

d1$Va[d1$Es=="CZE_2006"&d1$Vn==2030002]=1#
d1$Va[d1$Es=="CZE_2006"&d1$Vn==2030001]=2
d1$Va[d1$Es=="CZE_2006"&d1$Vn==2030004]=3
d1$Va[d1$Es=="CZE_2006"&d1$Vn==2030003]=4
d1$Va[d1$Es=="CZE_2006"&d1$Vn==2030005]=5
d1$Va[d1$Es=="CZE_2006"&d1$Vn==2030074]=6

d1$Va[d1$Es=="CZE_2010"&d1$Vn==2030001]=1#
d1$Va[d1$Es=="CZE_2010"&d1$Vn==2030002]=2
d1$Va[d1$Es=="CZE_2010"&d1$Vn==2030010]=3
d1$Va[d1$Es=="CZE_2010"&d1$Vn==2030004]=4
d1$Va[d1$Es=="CZE_2010"&d1$Vn==2030014]=5
d1$Va[d1$Es=="CZE_2010"&d1$Vn==2030003]=6

d1$Va[d1$Es=="CZE_2013"&d1$Vn==2030001]=1#
d1$Va[d1$Es=="CZE_2013"&d1$Vn==2030006]=2
d1$Va[d1$Es=="CZE_2013"&d1$Vn==2030004]=3
d1$Va[d1$Es=="CZE_2013"&d1$Vn==2030010]=4
d1$Va[d1$Es=="CZE_2013"&d1$Vn==2030002]=5
d1$Va[d1$Es=="CZE_2013"&d1$Vn==2030013]=6

d1$Va[d1$Es=="DEU_1998"&d1$Vn==2760004]=1#DEU:scale up letters after removing CSU
d1$Va[d1$Es=="DEU_1998"&d1$Vn==2760001]=2##note: no need to change 2760001 since it is already vote for CDU, or CSU in Bavaria
d1$Va[d1$Es=="DEU_1998"&d1$Vn==2760005]=3
d1$Va[d1$Es=="DEU_1998"&d1$Vn==2760006]=4
d1$Va[d1$Es=="DEU_1998"&d1$Vn==2760007]=5

d1$Va[d1$Es=="DEU_2005"&d1$Vn==2760004]=1#(comes from module 3, not module 2)
d1$Va[d1$Es=="DEU_2005"&d1$Vn==2760001]=2
d1$Va[d1$Es=="DEU_2005"&d1$Vn==2760006]=3
d1$Va[d1$Es=="DEU_2005"&d1$Vn==2760007]=4
d1$Va[d1$Es=="DEU_2005"&d1$Vn==2760005]=5
d1$Va[d1$Es=="DEU_2005"&d1$Vn==2760009]=6

d1$Va[d1$Es=="DEU_2009"&d1$Vn==2760001]=1#
d1$Va[d1$Es=="DEU_2009"&d1$Vn==2760004]=2
d1$Va[d1$Es=="DEU_2009"&d1$Vn==2760006]=3
d1$Va[d1$Es=="DEU_2009"&d1$Vn==2760007]=4
d1$Va[d1$Es=="DEU_2009"&d1$Vn==2760005]=5

d1$Va[d1$Es=="DEU_2013"&d1$Vn==2760001]=1#
d1$Va[d1$Es=="DEU_2013"&d1$Vn==2760004]=2
d1$Va[d1$Es=="DEU_2013"&d1$Vn==2760007]=3
d1$Va[d1$Es=="DEU_2013"&d1$Vn==2760005]=4
d1$Va[d1$Es=="DEU_2013"&d1$Vn==2760006]=5

d1$Va[d1$Es=="DNK_1998"&d1$Vn==2080001]=1#
d1$Va[d1$Es=="DNK_1998"&d1$Vn==2080002]=2
d1$Va[d1$Es=="DNK_1998"&d1$Vn==2080007]=3
d1$Va[d1$Es=="DNK_1998"&d1$Vn==2080004]=4
d1$Va[d1$Es=="DNK_1998"&d1$Vn==2080010]=5
d1$Va[d1$Es=="DNK_1998"&d1$Vn==2080011]=6

d1$Va[d1$Es=="DNK_2007"&d1$Vn==2080002]=1#
d1$Va[d1$Es=="DNK_2007"&d1$Vn==2080001]=2
d1$Va[d1$Es=="DNK_2007"&d1$Vn==2080010]=3
d1$Va[d1$Es=="DNK_2007"&d1$Vn==2080004]=4
d1$Va[d1$Es=="DNK_2007"&d1$Vn==2080007]=5
d1$Va[d1$Es=="DNK_2007"&d1$Vn==2080006]=6

d1$Va[d1$Es=="ESP_1996"&d1$Vn==7240002]=1#
d1$Va[d1$Es=="ESP_1996"&d1$Vn==7240001]=2
d1$Va[d1$Es=="ESP_1996"&d1$Vn==7240004]=3
d1$Va[d1$Es=="ESP_1996"&d1$Vn==7240003]=4
d1$Va[d1$Es=="ESP_1996"&d1$Vn==7240005]=5

d1$Va[d1$Es=="ESP_2000"&d1$Vn==7240002]=1#
d1$Va[d1$Es=="ESP_2000"&d1$Vn==7240001]=2
d1$Va[d1$Es=="ESP_2000"&d1$Vn==7240004]=3
d1$Va[d1$Es=="ESP_2000"&d1$Vn==7240003]=4
d1$Va[d1$Es=="ESP_2000"&d1$Vn==7240005]=5

d1$Va[d1$Es=="ESP_2008"&d1$Vn==7240001]=1#
d1$Va[d1$Es=="ESP_2008"&d1$Vn==7240002]=2
d1$Va[d1$Es=="ESP_2008"&d1$Vn==7240003]=3
d1$Va[d1$Es=="ESP_2008"&d1$Vn==7240005]=4
d1$Va[d1$Es=="ESP_2008"&d1$Vn==7240009]=5
d1$Va[d1$Es=="ESP_2008"&d1$Vn==7240004]=6

d1$Va[d1$Es=="FIN_2007"&d1$Vn==2460001]=1#
d1$Va[d1$Es=="FIN_2007"&d1$Vn==2460003]=2
d1$Va[d1$Es=="FIN_2007"&d1$Vn==2460002]=3
d1$Va[d1$Es=="FIN_2007"&d1$Vn==2460006]=4
d1$Va[d1$Es=="FIN_2007"&d1$Vn==2460005]=5
d1$Va[d1$Es=="FIN_2007"&d1$Vn==2460008]=6

d1$Va[d1$Es=="FIN_2011"&d1$Vn==2460003]=1#
d1$Va[d1$Es=="FIN_2011"&d1$Vn==2460002]=2
d1$Va[d1$Es=="FIN_2011"&d1$Vn==2460004]=3
d1$Va[d1$Es=="FIN_2011"&d1$Vn==2460001]=4
d1$Va[d1$Es=="FIN_2011"&d1$Vn==2460006]=5
d1$Va[d1$Es=="FIN_2011"&d1$Vn==2460005]=6

d1$Va[d1$Es=="FIN_2015"&d1$Vn==2460001]=1#
d1$Va[d1$Es=="FIN_2015"&d1$Vn==2460003]=2
d1$Va[d1$Es=="FIN_2015"&d1$Vn==2460004]=3
d1$Va[d1$Es=="FIN_2015"&d1$Vn==2460002]=4
d1$Va[d1$Es=="FIN_2015"&d1$Vn==2460005]=5
d1$Va[d1$Es=="FIN_2015"&d1$Vn==2460006]=6

d1$Va[d1$Es=="FRA_2007"&d1$Vn==2500001]=1#
d1$Va[d1$Es=="FRA_2007"&d1$Vn==2500002]=2
d1$Va[d1$Es=="FRA_2007"&d1$Vn==2500003]=3
d1$Va[d1$Es=="FRA_2007"&d1$Vn==2500007]=4
d1$Va[d1$Es=="FRA_2007"&d1$Vn==2500004]=5
d1$Va[d1$Es=="FRA_2007"&d1$Vn==2500005]=6

d1$Va[d1$Es=="GBR_1997"&d1$Vn==8260002]=1#
d1$Va[d1$Es=="GBR_1997"&d1$Vn==8260001]=2
d1$Va[d1$Es=="GBR_1997"&d1$Vn==8260003]=3
d1$Va[d1$Es=="GBR_1997"&d1$Vn==8260004]=4
d1$Va[d1$Es=="GBR_1997"&d1$Vn==8260006]=5

d1$Va[d1$Es=="GBR_2015"&d1$Vn==8260001]=1#
d1$Va[d1$Es=="GBR_2015"&d1$Vn==8260002]=2
d1$Va[d1$Es=="GBR_2015"&d1$Vn==8260007]=3
d1$Va[d1$Es=="GBR_2015"&d1$Vn==8260003]=4
d1$Va[d1$Es=="GBR_2015"&d1$Vn==8260004]=5
d1$Va[d1$Es=="GBR_2015"&d1$Vn==8260005]=6

d1$Va[d1$Es=="GRC_2009"&d1$Vn==3000001]=1#
d1$Va[d1$Es=="GRC_2009"&d1$Vn==3000002]=2
d1$Va[d1$Es=="GRC_2009"&d1$Vn==3000004]=3
d1$Va[d1$Es=="GRC_2009"&d1$Vn==3000019]=4
d1$Va[d1$Es=="GRC_2009"&d1$Vn==3000003]=5
d1$Va[d1$Es=="GRC_2009"&d1$Vn==3000005]=6

d1$Va[d1$Es=="GRC_2012"&d1$Vn==3000002]=1#
d1$Va[d1$Es=="GRC_2012"&d1$Vn==3000003]=2
d1$Va[d1$Es=="GRC_2012"&d1$Vn==3000001]=3
d1$Va[d1$Es=="GRC_2012"&d1$Vn==3000007]=4
d1$Va[d1$Es=="GRC_2012"&d1$Vn==3000006]=5
d1$Va[d1$Es=="GRC_2012"&d1$Vn==3000012]=6

d1$Va[d1$Es=="GRC_2015"&d1$Vn==3000003]=1#
d1$Va[d1$Es=="GRC_2015"&d1$Vn==3000002]=2
d1$Va[d1$Es=="GRC_2015"&d1$Vn==3000006]=3
d1$Va[d1$Es=="GRC_2015"&d1$Vn==3000013]=4
d1$Va[d1$Es=="GRC_2015"&d1$Vn==3000004]=5
d1$Va[d1$Es=="GRC_2015"&d1$Vn==3000007]=6

d1$Va[d1$Es=="HRV_2007"&d1$Vn==1910001]=1#
d1$Va[d1$Es=="HRV_2007"&d1$Vn==1910002]=2
d1$Va[d1$Es=="HRV_2007"&d1$Vn==1910003]=3
d1$Va[d1$Es=="HRV_2007"&d1$Vn==1910007]=4
d1$Va[d1$Es=="HRV_2007"&d1$Vn==1910006]=5
d1$Va[d1$Es=="HRV_2007"&d1$Vn==1910010]=6

d1$Va[d1$Es=="IRL_2007"&d1$Vn==3720001]=1#
d1$Va[d1$Es=="IRL_2007"&d1$Vn==3720002]=2
d1$Va[d1$Es=="IRL_2007"&d1$Vn==3720003]=3
d1$Va[d1$Es=="IRL_2007"&d1$Vn==3720005]=4
d1$Va[d1$Es=="IRL_2007"&d1$Vn==3720004]=5
d1$Va[d1$Es=="IRL_2007"&d1$Vn==3720006]=6

d1$Va[d1$Es=="IRL_2011"&d1$Vn==3720002]=1#E not coded, substituted with F
d1$Va[d1$Es=="IRL_2011"&d1$Vn==3720003]=2
d1$Va[d1$Es=="IRL_2011"&d1$Vn==3720001]=3
d1$Va[d1$Es=="IRL_2011"&d1$Vn==3720004]=4
d1$Va[d1$Es=="IRL_2011"&d1$Vn==3720005]=5

d1$Va[d1$Es=="ISL_1999"&d1$Vn==3520001]=1#
d1$Va[d1$Es=="ISL_1999"&d1$Vn==3520002]=2
d1$Va[d1$Es=="ISL_1999"&d1$Vn==3520003]=3
d1$Va[d1$Es=="ISL_1999"&d1$Vn==3520005]=4
d1$Va[d1$Es=="ISL_1999"&d1$Vn==3520004]=5

d1$Va[d1$Es=="ISL_2007"&d1$Vn==3520001]=1#
d1$Va[d1$Es=="ISL_2007"&d1$Vn==3520002]=2
d1$Va[d1$Es=="ISL_2007"&d1$Vn==3520005]=3
d1$Va[d1$Es=="ISL_2007"&d1$Vn==3520003]=4
d1$Va[d1$Es=="ISL_2007"&d1$Vn==3520004]=5
d1$Va[d1$Es=="ISL_2007"&d1$Vn==3520013]=6

d1$Va[d1$Es=="ISL_2009"&d1$Vn==3520002]=1#
d1$Va[d1$Es=="ISL_2009"&d1$Vn==3520001]=2
d1$Va[d1$Es=="ISL_2009"&d1$Vn==3520005]=3
d1$Va[d1$Es=="ISL_2009"&d1$Vn==3520003]=4
d1$Va[d1$Es=="ISL_2009"&d1$Vn==3520009]=5
d1$Va[d1$Es=="ISL_2009"&d1$Vn==3520004]=6

d1$Va[d1$Es=="ISL_2013"&d1$Vn==3520001]=1#
d1$Va[d1$Es=="ISL_2013"&d1$Vn==3520003]=2
d1$Va[d1$Es=="ISL_2013"&d1$Vn==3520002]=3
d1$Va[d1$Es=="ISL_2013"&d1$Vn==3520005]=4
d1$Va[d1$Es=="ISL_2013"&d1$Vn==3520007]=5
d1$Va[d1$Es=="ISL_2013"&d1$Vn==3520008]=6

d1$Va[d1$Es=="ISR_1996"&d1$Vn==3760002]=1#
d1$Va[d1$Es=="ISR_1996"&d1$Vn==3760001]=2
d1$Va[d1$Es=="ISR_1996"&d1$Vn==3760003]=3
d1$Va[d1$Es=="ISR_1996"&d1$Vn==3760014]=4
d1$Va[d1$Es=="ISR_1996"&d1$Vn==3760004]=5

d1$Va[d1$Es=="ISR_2006"&d1$Vn==3760016]=1#
d1$Va[d1$Es=="ISR_2006"&d1$Vn==3760002]=2
d1$Va[d1$Es=="ISR_2006"&d1$Vn==3760003]=3
d1$Va[d1$Es=="ISR_2006"&d1$Vn==3760001]=4
d1$Va[d1$Es=="ISR_2006"&d1$Vn==3760021]=5
d1$Va[d1$Es=="ISR_2006"&d1$Vn==3760050]=6

d1$Va[d1$Es=="ISR_2013"&d1$Vn==3760011]=1#F not coded, substituted with G
d1$Va[d1$Es=="ISR_2013"&d1$Vn==3760007]=2
d1$Va[d1$Es=="ISR_2013"&d1$Vn==3760002]=3
d1$Va[d1$Es=="ISR_2013"&d1$Vn==3760008]=4
d1$Va[d1$Es=="ISR_2013"&d1$Vn==3760003]=5
d1$Va[d1$Es=="ISR_2013"&d1$Vn==3760009]=6

d1$Va[d1$Es=="LVA_2010"&d1$Vn==4280002]=1#
d1$Va[d1$Es=="LVA_2010"&d1$Vn==4280001]=2
d1$Va[d1$Es=="LVA_2010"&d1$Vn==4280005]=3
d1$Va[d1$Es=="LVA_2010"&d1$Vn==4280003]=4
d1$Va[d1$Es=="LVA_2010"&d1$Vn==4280004]=5
d1$Va[d1$Es=="LVA_2010"&d1$Vn==4280006]=6

d1$Va[d1$Es=="LVA_2011"&d1$Vn==4280001]=1#
d1$Va[d1$Es=="LVA_2011"&d1$Vn==4280007]=2
d1$Va[d1$Es=="LVA_2011"&d1$Vn==4280002]=3
d1$Va[d1$Es=="LVA_2011"&d1$Vn==4280003]=4
d1$Va[d1$Es=="LVA_2011"&d1$Vn==4280005]=5
d1$Va[d1$Es=="LVA_2011"&d1$Vn==4280036]=6

d1$Va[d1$Es=="LVA_2014"&d1$Vn==4280001]=1#
d1$Va[d1$Es=="LVA_2014"&d1$Vn==4280002]=2
d1$Va[d1$Es=="LVA_2014"&d1$Vn==4280005]=3
d1$Va[d1$Es=="LVA_2014"&d1$Vn==4280003]=4
d1$Va[d1$Es=="LVA_2014"&d1$Vn==4280013]=5
d1$Va[d1$Es=="LVA_2014"&d1$Vn==4280014]=6

d1$Va[d1$Es=="MEX_2000"&d1$Vn==4840022]=1#
d1$Va[d1$Es=="MEX_2000"&d1$Vn==4840001]=2
d1$Va[d1$Es=="MEX_2000"&d1$Vn==4840023]=3

d1$Va[d1$Es=="MEX_2006"&d1$Vn==4840002]=1#
d1$Va[d1$Es=="MEX_2006"&d1$Vn==4840003]=2
d1$Va[d1$Es=="MEX_2006"&d1$Vn==4840001]=3

d1$Va[d1$Es=="MEX_2009"&d1$Vn==4840001]=1#
d1$Va[d1$Es=="MEX_2009"&d1$Vn==4840002]=2
d1$Va[d1$Es=="MEX_2009"&d1$Vn==4840003]=3

d1$Va[d1$Es=="MNE_2012"&d1$Vn==4990011]=1#
d1$Va[d1$Es=="MNE_2012"&d1$Vn==4990012]=2
d1$Va[d1$Es=="MNE_2012"&d1$Vn==4990003]=3
d1$Va[d1$Es=="MNE_2012"&d1$Vn==4990006]=4
d1$Va[d1$Es=="MNE_2012"&d1$Vn==4990009]=5
d1$Va[d1$Es=="MNE_2012"&d1$Vn==4990037]=6

d1$Va[d1$Es=="NLD_1998"&d1$Vn==5280002]=1#
d1$Va[d1$Es=="NLD_1998"&d1$Vn==5280003]=2
d1$Va[d1$Es=="NLD_1998"&d1$Vn==5280001]=3
d1$Va[d1$Es=="NLD_1998"&d1$Vn==5280004]=4
d1$Va[d1$Es=="NLD_1998"&d1$Vn==5280005]=5
d1$Va[d1$Es=="NLD_1998"&d1$Vn==5280006]=6

d1$Va[d1$Es=="NLD_2006"&d1$Vn==5280001]=1#
d1$Va[d1$Es=="NLD_2006"&d1$Vn==5280002]=2
d1$Va[d1$Es=="NLD_2006"&d1$Vn==5280006]=3
d1$Va[d1$Es=="NLD_2006"&d1$Vn==5280003]=4
d1$Va[d1$Es=="NLD_2006"&d1$Vn==5280009]=5
d1$Va[d1$Es=="NLD_2006"&d1$Vn==5280005]=6

d1$Va[d1$Es=="NLD_2010"&d1$Vn==5280003]=1#
d1$Va[d1$Es=="NLD_2010"&d1$Vn==5280002]=2
d1$Va[d1$Es=="NLD_2010"&d1$Vn==5280009]=3
d1$Va[d1$Es=="NLD_2010"&d1$Vn==5280001]=4
d1$Va[d1$Es=="NLD_2010"&d1$Vn==5280006]=5
d1$Va[d1$Es=="NLD_2010"&d1$Vn==5280004]=6

d1$Va[d1$Es=="NOR_1997"&d1$Vn==5780001]=1#
d1$Va[d1$Es=="NOR_1997"&d1$Vn==5780003]=2
d1$Va[d1$Es=="NOR_1997"&d1$Vn==5780002]=3
d1$Va[d1$Es=="NOR_1997"&d1$Vn==5780007]=4
d1$Va[d1$Es=="NOR_1997"&d1$Vn==5780006]=5
d1$Va[d1$Es=="NOR_1997"&d1$Vn==5780004]=6

d1$Va[d1$Es=="NOR_2005"&d1$Vn==5780001]=1#
d1$Va[d1$Es=="NOR_2005"&d1$Vn==5780003]=2
d1$Va[d1$Es=="NOR_2005"&d1$Vn==5780002]=3
d1$Va[d1$Es=="NOR_2005"&d1$Vn==5780004]=4
d1$Va[d1$Es=="NOR_2005"&d1$Vn==5780007]=5
d1$Va[d1$Es=="NOR_2005"&d1$Vn==5780006]=6

d1$Va[d1$Es=="NOR_2009"&d1$Vn==5780001]=1#
d1$Va[d1$Es=="NOR_2009"&d1$Vn==5780003]=2
d1$Va[d1$Es=="NOR_2009"&d1$Vn==5780002]=3
d1$Va[d1$Es=="NOR_2009"&d1$Vn==5780004]=4
d1$Va[d1$Es=="NOR_2009"&d1$Vn==5780006]=5
d1$Va[d1$Es=="NOR_2009"&d1$Vn==5780007]=6

d1$Va[d1$Es=="NOR_2013"&d1$Vn==5780001]=1#
d1$Va[d1$Es=="NOR_2013"&d1$Vn==5780002]=2
d1$Va[d1$Es=="NOR_2013"&d1$Vn==5780003]=3
d1$Va[d1$Es=="NOR_2013"&d1$Vn==5780007]=4
d1$Va[d1$Es=="NOR_2013"&d1$Vn==5780006]=5
d1$Va[d1$Es=="NOR_2013"&d1$Vn==5780008]=6

d1$Va[d1$Es=="NZL_1996"&d1$Vn==5540001]=1#
d1$Va[d1$Es=="NZL_1996"&d1$Vn==5540002]=2
d1$Va[d1$Es=="NZL_1996"&d1$Vn==5540003]=3
d1$Va[d1$Es=="NZL_1996"&d1$Vn==5540009]=4
d1$Va[d1$Es=="NZL_1996"&d1$Vn==5540004]=5
d1$Va[d1$Es=="NZL_1996"&d1$Vn==5540037]=6

d1$Va[d1$Es=="NZL_2008"&d1$Vn==5540001]=1#
d1$Va[d1$Es=="NZL_2008"&d1$Vn==5540002]=2
d1$Va[d1$Es=="NZL_2008"&d1$Vn==5540005]=3
d1$Va[d1$Es=="NZL_2008"&d1$Vn==5540003]=4
d1$Va[d1$Es=="NZL_2008"&d1$Vn==5540004]=5
d1$Va[d1$Es=="NZL_2008"&d1$Vn==5540008]=6

d1$Va[d1$Es=="NZL_2011"&d1$Vn==5540001]=1#
d1$Va[d1$Es=="NZL_2011"&d1$Vn==5540002]=2
d1$Va[d1$Es=="NZL_2011"&d1$Vn==5540005]=3
d1$Va[d1$Es=="NZL_2011"&d1$Vn==5540003]=4
d1$Va[d1$Es=="NZL_2011"&d1$Vn==5540007]=5
d1$Va[d1$Es=="NZL_2011"&d1$Vn==5540008]=6

d1$Va[d1$Es=="NZL_2014"&d1$Vn==5540001]=1#
d1$Va[d1$Es=="NZL_2014"&d1$Vn==5540002]=2
d1$Va[d1$Es=="NZL_2014"&d1$Vn==5540005]=3
d1$Va[d1$Es=="NZL_2014"&d1$Vn==5540003]=4
d1$Va[d1$Es=="NZL_2014"&d1$Vn==5540007]=5
d1$Va[d1$Es=="NZL_2014"&d1$Vn==5540010]=6

d1$Va[d1$Es=="PER_2001"&d1$Vn==6040008]=1#
d1$Va[d1$Es=="PER_2001"&d1$Vn==6040026]=2
d1$Va[d1$Es=="PER_2001"&d1$Vn==6040052]=3
d1$Va[d1$Es=="PER_2001"&d1$Vn==6040023]=4
d1$Va[d1$Es=="PER_2001"&d1$Vn==6040030]=5
d1$Va[d1$Es=="PER_2001"&d1$Vn==6040027]=6

d1$Va[d1$Es=="PER_2011"&d1$Vn==6040050]=1#
d1$Va[d1$Es=="PER_2011"&d1$Vn==6040001]=2
d1$Va[d1$Es=="PER_2011"&d1$Vn==6040008]=3
d1$Va[d1$Es=="PER_2011"&d1$Vn==6040051]=4
d1$Va[d1$Es=="PER_2011"&d1$Vn==6040015]=5

d1$Va[d1$Es=="PER_2016"&d1$Vn==6040001]=1#
d1$Va[d1$Es=="PER_2016"&d1$Vn==6040002]=2
d1$Va[d1$Es=="PER_2016"&d1$Vn==6040003]=3
d1$Va[d1$Es=="PER_2016"&d1$Vn==6040004]=4
d1$Va[d1$Es=="PER_2016"&d1$Vn==6040005]=5
d1$Va[d1$Es=="PER_2016"&d1$Vn==6040006]=6

d1$Va[d1$Es=="POL_1997"&d1$Vn==6160004]=1#
d1$Va[d1$Es=="POL_1997"&d1$Vn==6160009]=2
d1$Va[d1$Es=="POL_1997"&d1$Vn==6160008]=3
d1$Va[d1$Es=="POL_1997"&d1$Vn==6160003]=4
d1$Va[d1$Es=="POL_1997"&d1$Vn==6160010]=5
d1$Va[d1$Es=="POL_1997"&d1$Vn==6160007]=6

d1$Va[d1$Es=="POL_2005"&d1$Vn==6160002]=1#F not coded, substituted with H
d1$Va[d1$Es=="POL_2005"&d1$Vn==6160001]=2
d1$Va[d1$Es=="POL_2005"&d1$Vn==6160024]=3
d1$Va[d1$Es=="POL_2005"&d1$Vn==6160009]=4
d1$Va[d1$Es=="POL_2005"&d1$Vn==6160014]=5
d1$Va[d1$Es=="POL_2005"&d1$Vn==6160068]=6

d1$Va[d1$Es=="POL_2007"&d1$Vn==6160001]=1#
d1$Va[d1$Es=="POL_2007"&d1$Vn==6160002]=2
d1$Va[d1$Es=="POL_2007"&d1$Vn==6160020]=3
d1$Va[d1$Es=="POL_2007"&d1$Vn==6160003]=4

d1$Va[d1$Es=="POL_2011"&d1$Vn==6160001]=1#
d1$Va[d1$Es=="POL_2011"&d1$Vn==6160002]=2
d1$Va[d1$Es=="POL_2011"&d1$Vn==6160017]=3
d1$Va[d1$Es=="POL_2011"&d1$Vn==6160003]=4
d1$Va[d1$Es=="POL_2011"&d1$Vn==6160009]=5

d1$Va[d1$Es=="PRT_2002"&d1$Vn==6200001]=1#
d1$Va[d1$Es=="PRT_2002"&d1$Vn==6200002]=2
d1$Va[d1$Es=="PRT_2002"&d1$Vn==6200003]=3
d1$Va[d1$Es=="PRT_2002"&d1$Vn==6200004]=4
d1$Va[d1$Es=="PRT_2002"&d1$Vn==6200006]=5
d1$Va[d1$Es=="PRT_2002"&d1$Vn==6200007]=6

d1$Va[d1$Es=="PRT_2009"&d1$Vn==6200002]=1#
d1$Va[d1$Es=="PRT_2009"&d1$Vn==6200001]=2
d1$Va[d1$Es=="PRT_2009"&d1$Vn==6200003]=3
d1$Va[d1$Es=="PRT_2009"&d1$Vn==6200006]=4
d1$Va[d1$Es=="PRT_2009"&d1$Vn==6200004]=5

d1$Va[d1$Es=="PRT_2015"&d1$Vn==6200011]=1#
d1$Va[d1$Es=="PRT_2015"&d1$Vn==6200002]=2
d1$Va[d1$Es=="PRT_2015"&d1$Vn==6200006]=3
d1$Va[d1$Es=="PRT_2015"&d1$Vn==6200004]=4
d1$Va[d1$Es=="PRT_2015"&d1$Vn==6200020]=5

d1$Va[d1$Es=="SRB_2012"&d1$Vn==6880001]=1#
d1$Va[d1$Es=="SRB_2012"&d1$Vn==6880002]=2
d1$Va[d1$Es=="SRB_2012"&d1$Vn==6880004]=3
d1$Va[d1$Es=="SRB_2012"&d1$Vn==6880005]=4
d1$Va[d1$Es=="SRB_2012"&d1$Vn==6880006]=5
d1$Va[d1$Es=="SRB_2012"&d1$Vn==6880007]=6

d1$Va[d1$Es=="SVK_2010"&d1$Vn==7030001]=1#
d1$Va[d1$Es=="SVK_2010"&d1$Vn==7030003]=2
d1$Va[d1$Es=="SVK_2010"&d1$Vn==7030002]=3
d1$Va[d1$Es=="SVK_2010"&d1$Vn==7030004]=4
d1$Va[d1$Es=="SVK_2010"&d1$Vn==7030008]=5
d1$Va[d1$Es=="SVK_2010"&d1$Vn==7030006]=6

d1$Va[d1$Es=="SVK_2016"&d1$Vn==7030001]=1#
d1$Va[d1$Es=="SVK_2016"&d1$Vn==7030002]=2
d1$Va[d1$Es=="SVK_2016"&d1$Vn==7030009]=3
d1$Va[d1$Es=="SVK_2016"&d1$Vn==7030006]=4
d1$Va[d1$Es=="SVK_2016"&d1$Vn==7030012]=5
d1$Va[d1$Es=="SVK_2016"&d1$Vn==7030014]=6

d1$Va[d1$Es=="SVN_1996"&d1$Vn==7050004]=1#F not coded, substituted with G
d1$Va[d1$Es=="SVN_1996"&d1$Vn==7050006]=2
d1$Va[d1$Es=="SVN_1996"&d1$Vn==7050001]=3
d1$Va[d1$Es=="SVN_1996"&d1$Vn==7050017]=4
d1$Va[d1$Es=="SVN_1996"&d1$Vn==7050002]=5
d1$Va[d1$Es=="SVN_1996"&d1$Vn==7050008]=6

d1$Va[d1$Es=="SVN_2008"&d1$Vn==7050002]=1#
d1$Va[d1$Es=="SVN_2008"&d1$Vn==7050001]=2
d1$Va[d1$Es=="SVN_2008"&d1$Vn==7050020]=3
d1$Va[d1$Es=="SVN_2008"&d1$Vn==7050003]=4
d1$Va[d1$Es=="SVN_2008"&d1$Vn==7050008]=5
d1$Va[d1$Es=="SVN_2008"&d1$Vn==7050019]=6

d1$Va[d1$Es=="SVN_2011"&d1$Vn==7050016]=1#
d1$Va[d1$Es=="SVN_2011"&d1$Vn==7050001]=2
d1$Va[d1$Es=="SVN_2011"&d1$Vn==7050002]=3
d1$Va[d1$Es=="SVN_2011"&d1$Vn==7050009]=4
d1$Va[d1$Es=="SVN_2011"&d1$Vn==7050003]=5
d1$Va[d1$Es=="SVN_2011"&d1$Vn==7050006]=6

d1$Va[d1$Es=="SWE_1998"&d1$Vn==7520001]=1#
d1$Va[d1$Es=="SWE_1998"&d1$Vn==7520002]=2
d1$Va[d1$Es=="SWE_1998"&d1$Vn==7520007]=3
d1$Va[d1$Es=="SWE_1998"&d1$Vn==7520004]=4
d1$Va[d1$Es=="SWE_1998"&d1$Vn==7520006]=5
d1$Va[d1$Es=="SWE_1998"&d1$Vn==7520003]=6

d1$Va[d1$Es=="SWE_2006"&d1$Vn==7520001]=1#
d1$Va[d1$Es=="SWE_2006"&d1$Vn==7520002]=2
d1$Va[d1$Es=="SWE_2006"&d1$Vn==7520006]=3
d1$Va[d1$Es=="SWE_2006"&d1$Vn==7520003]=4
d1$Va[d1$Es=="SWE_2006"&d1$Vn==7520004]=5
d1$Va[d1$Es=="SWE_2006"&d1$Vn==7520007]=6

d1$Va[d1$Es=="SWE_2014"&d1$Vn==7520001]=1#
d1$Va[d1$Es=="SWE_2014"&d1$Vn==7520002]=2
d1$Va[d1$Es=="SWE_2014"&d1$Vn==7520008]=3
d1$Va[d1$Es=="SWE_2014"&d1$Vn==7520005]=4
d1$Va[d1$Es=="SWE_2014"&d1$Vn==7520006]=5
d1$Va[d1$Es=="SWE_2014"&d1$Vn==7520007]=6

d1$Va[d1$Es=="THA_2007"&d1$Vn==7640014]=1
d1$Va[d1$Es=="THA_2007"&d1$Vn==7640002]=2
d1$Va[d1$Es=="THA_2007"&d1$Vn==7640021]=3
d1$Va[d1$Es=="THA_2007"&d1$Vn==7640008]=4
d1$Va[d1$Es=="THA_2007"&d1$Vn==7640060]=5
d1$Va[d1$Es=="THA_2007"&d1$Vn==7640024]=6

d1$Va[d1$Es=="THA_2011"&d1$Vn==7640001]=1
d1$Va[d1$Es=="THA_2011"&d1$Vn==7640002]=2
d1$Va[d1$Es=="THA_2011"&d1$Vn==7640003]=3
d1$Va[d1$Es=="THA_2011"&d1$Vn==7640004]=4
d1$Va[d1$Es=="THA_2011"&d1$Vn==7640005]=5
d1$Va[d1$Es=="THA_2011"&d1$Vn==7640006]=6

d1$Va[d1$Es=="TWN_1996"&d1$Vn==1580002]=1#
d1$Va[d1$Es=="TWN_1996"&d1$Vn==1580001]=2
d1$Va[d1$Es=="TWN_1996"&d1$Vn==1580004]=3

d1$Va[d1$Es=="ZAF_2009"&d1$Vn==7100001]=1#
d1$Va[d1$Es=="ZAF_2009"&d1$Vn==7100002]=2
d1$Va[d1$Es=="ZAF_2009"&d1$Vn==7100008]=3
d1$Va[d1$Es=="ZAF_2009"&d1$Vn==7100004]=4
d1$Va[d1$Es=="ZAF_2009"&d1$Vn==7100022]=5
d1$Va[d1$Es=="ZAF_2009"&d1$Vn==7100006]=6

d1$Va[d1$Es=="ZAF_2014"&d1$Vn==7100001]=1#
d1$Va[d1$Es=="ZAF_2014"&d1$Vn==7100002]=2
d1$Va[d1$Es=="ZAF_2014"&d1$Vn==7100003]=3
d1$Va[d1$Es=="ZAF_2014"&d1$Vn==7100004]=4
d1$Va[d1$Es=="ZAF_2014"&d1$Vn==7100005]=5
d1$Va[d1$Es=="ZAF_2014"&d1$Vn==7100006]=6

colnames(d1)[colnames(d1) %in% paste0("IMD3008_", LETTERS[1:6])]=paste0("Pl.", 1:6)#party likability [0,10] A-F into 1:6----
d1[,names(d1) %in% paste0("Pl.", 1:6)][d1[,names(d1) %in% paste0("Pl.", 1:6)]>10]=NA#coding missing as NAs

colnames(d1)[colnames(d1) %in% paste0("IMD3009_", LETTERS[1:6])]=paste0("Ll.", 1:6)#leader likability [0,10] A-F into 1:6
d1[,names(d1) %in% paste0("Ll.", 1:6)][d1[,names(d1) %in% paste0("Ll.", 1:6)]>10]=NA#coding missing as NAs

colnames(d1)[colnames(d1) %in% paste0("IMD5002_", LETTERS[1:6])]=paste0("Ps.", 1:6)#party seats A-F into 1:6 PERCENT SEATS - LOWER HOUSE
d1[,names(d1) %in% paste0("Ps.", 1:6)][d1[,names(d1) %in% paste0("Ps.", 1:6)]>100]=NA#coding missing as NAs

colnames(d1)[colnames(d1) %in% paste0("IMD5001_", LETTERS[1:6])]=paste0("Pv.", 1:6)#party votes A-F into 1:6 PERCENT VOTE - LOWER HOUSE
d1[,names(d1) %in% paste0("Pv.", 1:6)][d1[,names(d1) %in% paste0("Pv.", 1:6)]>100]=NA#coding missing as NAs

#EEeconS Election Economy Subjective:
d1$IMD3013_1[d1$IMD3013_1>5]=NA#coding missing as NAs
d1$IMD3013_2[d1$IMD3013_2>5]=NA#coding missing as NAs
d1$IMD3013_3[d1$IMD3013_3>5]=NA#coding missing as NAs
d1$EeconS=NA
d1$EeconS=ifelse(d1$IMD3013_1<3, d1$IMD3013_2,
                 (ifelse(d1$IMD3013_1==3, 3,
                         (ifelse(d1$IMD3013_1>3,d1$IMD3013_3,d1$EeconS)))))

#recode EeconS from CSES scale
# 1. GOTTEN BETTER c(1,2,3,4,5)
# 5. GOTTEN WORSE 
# to TEV scale, noting that the TEV used the conversion scale below for same elections:
# 0. Negative evaluation c(0, 0.25, 0.5, 0.75, 1)
# 1. Positive evaluation
d1$EeconS[d1$EeconS==1]=1
d1$EeconS[d1$EeconS==2]=0.75
d1$EeconS[d1$EeconS==3]=0.5
d1$EeconS[d1$EeconS==4]=0.25
d1$EeconS[d1$EeconS==5]=0

#Elections with EeconS:
table(d1$Es,d1$EeconS)
d1$Eecon=0
d1$Eecon[d1$Es %in% c("ARG_2015","AUS_1996","AUS_2013","AUT_2013","BGR_2014","CAN_1997","CZE_1996","CZE_2013","DEU_1998
                      ","DEU_2013","DNK_1998","ESP_1996","ESP_2000","FIN_2015","GBR_1997","GBR_2015
                      ","GRC_2012","GRC_2015","HUN_1998","IRL_2011","ISL_1999","ISL_2013","ISR_1996","ISR_2013","LVA_2011","LVA_2014","MEX_2000
                      ","MNE_2012","NLD_1998","NOR_1997","NOR_2013","NZL_1996","NZL_2011","NZL_2014
                      ","PER_2016","POL_1997","POL_2011","PRT_2002","PRT_2015","SRB_2012","SVK_2016","SVN_1996","SVN_2011","SWE_1998
                      ","SWE_2014","TWN_1996","ZAF_2014")]=1

#LRR=Left Right Respondent (voter LR self):
d1$LRR=d1$IMD3006
d1$LRR[d1$LRR>10]=NA
#recoding to 0-10 scale: already 0-10 scale
#LRP.i=Left Right Parties:
d1$LRP.1=d1$IMD3007_A#LEFT-RIGHT - PARTY A
d1$LRP.2=d1$IMD3007_B#LEFT-RIGHT - PARTY B
d1$LRP.3=d1$IMD3007_C#LEFT-RIGHT - PARTY C
d1$LRP.4=d1$IMD3007_D#LEFT-RIGHT - PARTY D
d1$LRP.5=d1$IMD3007_E#LEFT-RIGHT - PARTY E
d1$LRP.6=d1$IMD3007_F#LEFT-RIGHT - PARTY F
# d1$LRP.7=as.numeric(NA)
d1$LRP.1[d1$LRP.1>10]=NA
d1$LRP.2[d1$LRP.2>10]=NA
d1$LRP.3[d1$LRP.3>10]=NA
d1$LRP.4[d1$LRP.4>10]=NA
d1$LRP.5[d1$LRP.5>10]=NA
d1$LRP.6[d1$LRP.6>10]=NA
#recoding to 0-10 scale: already 0-10 scale
#Elections with LRR and LRP:
d1$LR=1
d1$LR[d1$Es %in% c("ARG_2015","PER_2016","NZL_2011","MNE_2012","CZE_2006")]=0

#Seats to vote ratios:
d1$P1SVrat=d1$Ps.1/d1$Pv.1#for most voted party
d1$P2SVrat=d1$Ps.2/d1$Pv.2#for second most voted party
d1$P1P2rat=d1$P1SVrat/d1$P2SVrat#ratio of rations first/second party
#Gallagher Disproportionality measure:
d1$GDM=NA
d1$GDM[d1$Es=="AUS_1996"]=10.97
d1$GDM[d1$Es=="AUS_2004"]=8.60
d1$GDM[d1$Es=="AUS_2007"]=10.27
d1$GDM[d1$Es=="AUS_2013"]=9.54
d1$GDM[d1$Es=="AUT_2008"]=2.92
d1$GDM[d1$Es=="AUT_2013"]=3.31
d1$GDM[d1$Es=="BGR_2014"]=2.52
d1$GDM[d1$Es=="CAN_1997"]=13.26
d1$GDM[d1$Es=="CAN_2008"]=10.09
d1$GDM[d1$Es=="CAN_2011"]=12.42
d1$GDM[d1$Es=="CAN_2015"]=12.01
d1$GDM[d1$Es=="CZE_1996"]=5.55
d1$GDM[d1$Es=="CZE_2006"]=5.72
d1$GDM[d1$Es=="CZE_2010"]=8.76
d1$GDM[d1$Es=="CZE_2013"]=6.12
d1$GDM[d1$Es=="DEU_1998"]=3.15#The figure in the effective number of legislative parties (Ns) column is based on 
d1$GDM[d1$Es=="DEU_2005"]=2.16#treating the CDU and CSU as separate parties.
d1$GDM[d1$Es=="DEU_2009"]=3.40#
d1$GDM[d1$Es=="DEU_2013"]=7.83#
d1$GDM[d1$Es=="DNK_1998"]=0.42
d1$GDM[d1$Es=="DNK_2007"]=0.72
d1$GDM[d1$Es=="ESP_1996"]=5.36
d1$GDM[d1$Es=="ESP_2000"]=6.10
d1$GDM[d1$Es=="ESP_2008"]=4.49
d1$GDM[d1$Es=="FIN_2007"]=3.20
d1$GDM[d1$Es=="FIN_2011"]=2.95
d1$GDM[d1$Es=="FIN_2015"]=3.03
d1$GDM[d1$Es=="FRA_2007"]=13.58
d1$GDM[d1$Es=="GBR_1997"]=16.51
d1$GDM[d1$Es=="GBR_2015"]=15.02
d1$GDM[d1$Es=="GRC_2009"]=7.29
d1$GDM[d1$Es=="GRC_2012"]=9.96#Jun
d1$GDM[d1$Es=="GRC_2015"]=9.98#Jan
d1$GDM[d1$Es=="IRL_2007"]=5.85
d1$GDM[d1$Es=="IRL_2011"]=8.69
d1$GDM[d1$Es=="ISL_1999"]=1.06
d1$GDM[d1$Es=="ISL_2007"]=3.49
d1$GDM[d1$Es=="ISL_2009"]=2.58
d1$GDM[d1$Es=="ISL_2013"]=6.23
d1$GDM[d1$Es=="ISR_1996"]=1.65
d1$GDM[d1$Es=="ISR_2006"]=2.49
d1$GDM[d1$Es=="ISR_2013"]=3.09
d1$GDM[d1$Es=="LVA_2010"]=2.80
d1$GDM[d1$Es=="LVA_2011"]=2.76
d1$GDM[d1$Es=="LVA_2014"]=2.30
d1$GDM[d1$Es=="MEX_2000"]=6.70#overall
d1$GDM[d1$Es=="MEX_2006"]=6.34#overall
d1$GDM[d1$Es=="MEX_2009"]=10.46#overall
d1$GDM[d1$Es=="MNE_2012"]=2.28
d1$GDM[d1$Es=="NLD_1998"]=1.28
d1$GDM[d1$Es=="NLD_2006"]=1.03
d1$GDM[d1$Es=="NLD_2010"]=0.81
d1$GDM[d1$Es=="NOR_1997"]=3.44
d1$GDM[d1$Es=="NOR_2005"]=2.67
d1$GDM[d1$Es=="NOR_2009"]=3.01
d1$GDM[d1$Es=="NOR_2013"]=2.56
d1$GDM[d1$Es=="NZL_1996"]=3.43
d1$GDM[d1$Es=="NZL_2008"]=3.84
d1$GDM[d1$Es=="NZL_2011"]=2.38
d1$GDM[d1$Es=="NZL_2014"]=3.72
d1$GDM[d1$Es=="POL_1997"]=10.63
d1$GDM[d1$Es=="POL_2005"]=6.97
d1$GDM[d1$Es=="POL_2007"]=4.67
d1$GDM[d1$Es=="POL_2011"]=5.95
d1$GDM[d1$Es=="PRT_2002"]=4.64
d1$GDM[d1$Es=="PRT_2009"]=5.63
d1$GDM[d1$Es=="PRT_2015"]=5.65
d1$GDM[d1$Es=="ROU_1996"]=6.34
d1$GDM[d1$Es=="SRB_2012"]=6.53
d1$GDM[d1$Es=="SVK_2010"]=7.46
d1$GDM[d1$Es=="SVK_2016"]=6.10
d1$GDM[d1$Es=="SVN_1996"]=3.59
d1$GDM[d1$Es=="SVN_2008"]=3.89
d1$GDM[d1$Es=="SVN_2011"]=3.64
d1$GDM[d1$Es=="SWE_1998"]=0.97
d1$GDM[d1$Es=="SWE_2006"]=3.02
d1$GDM[d1$Es=="SWE_2014"]=2.64
d1$GDM[d1$Es=="TWN_1996"]=#no data!!! check notes
d1$GDM[d1$Es=="ZAF_2009"]=0.30
d1$GDM[d1$Es=="ZAF_2014"]=0.37


#d1 assign CSES numerical ID to alternatives (CSESid.i):----
d1$CSESid.1=NA
d1$CSESid.2=NA
d1$CSESid.3=NA
d1$CSESid.4=NA
d1$CSESid.5=NA
d1$CSESid.6=NA
# d1$CSESid.7=NA

d1$CSESid.1[d1$Es=="ARG_2015"]=0320001
d1$CSESid.2[d1$Es=="ARG_2015"]=0320002
d1$CSESid.3[d1$Es=="ARG_2015"]=0320003
d1$CSESid.4[d1$Es=="ARG_2015"]=0320004
d1$CSESid.5[d1$Es=="ARG_2015"]=0320005
d1$CSESid.6[d1$Es=="ARG_2015"]=0320006

d1$CSESid.1[d1$Es=="AUS_1996"]=0360003
d1$CSESid.2[d1$Es=="AUS_1996"]=0360001
d1$CSESid.3[d1$Es=="AUS_1996"]=0360002
d1$CSESid.4[d1$Es=="AUS_1996"]=0360004

d1$CSESid.1[d1$Es=="AUS_2004"]=0360001#scaled up since E not coded
d1$CSESid.2[d1$Es=="AUS_2004"]=0360003
d1$CSESid.3[d1$Es=="AUS_2004"]=0360005
d1$CSESid.4[d1$Es=="AUS_2004"]=0360002
d1$CSESid.5[d1$Es=="AUS_2004"]=0360004
d1$CSESid.6[d1$Es=="AUS_2004"]=0360008

d1$CSESid.1[d1$Es=="AUS_2007"]=0360003
d1$CSESid.2[d1$Es=="AUS_2007"]=0360001
d1$CSESid.3[d1$Es=="AUS_2007"]=0360002
d1$CSESid.4[d1$Es=="AUS_2007"]=0360005

d1$CSESid.1[d1$Es=="AUS_2013"]=0360001#scaled up E since D not coded
d1$CSESid.2[d1$Es=="AUS_2013"]=0360003
d1$CSESid.3[d1$Es=="AUS_2013"]=0360005
d1$CSESid.4[d1$Es=="AUS_2013"]=0360002

d1$CSESid.1[d1$Es=="AUT_2013"]=0400001
d1$CSESid.2[d1$Es=="AUT_2013"]=0400002
d1$CSESid.3[d1$Es=="AUT_2013"]=0400003
d1$CSESid.4[d1$Es=="AUT_2013"]=0400005
d1$CSESid.5[d1$Es=="AUT_2013"]=0400004
d1$CSESid.6[d1$Es=="AUT_2013"]=0400010

d1$CSESid.1[d1$Es=="BGR_2014"]=1000001
d1$CSESid.2[d1$Es=="BGR_2014"]=1000020
d1$CSESid.3[d1$Es=="BGR_2014"]=1000003
d1$CSESid.4[d1$Es=="BGR_2014"]=1000004
d1$CSESid.5[d1$Es=="BGR_2014"]=1000005
d1$CSESid.6[d1$Es=="BGR_2014"]=1000006

d1$CSESid.1[d1$Es=="CAN_1997"]=1240001
d1$CSESid.2[d1$Es=="CAN_1997"]=1240007
d1$CSESid.3[d1$Es=="CAN_1997"]=1240006
d1$CSESid.4[d1$Es=="CAN_1997"]=1240003
d1$CSESid.5[d1$Es=="CAN_1997"]=1240004

d1$CSESid.1[d1$Es=="CAN_2008"]=1240002
d1$CSESid.2[d1$Es=="CAN_2008"]=1240001
d1$CSESid.3[d1$Es=="CAN_2008"]=1240004
d1$CSESid.4[d1$Es=="CAN_2008"]=1240003
d1$CSESid.5[d1$Es=="CAN_2008"]=1240005

d1$CSESid.1[d1$Es=="CAN_2011"]=1240002
d1$CSESid.2[d1$Es=="CAN_2011"]=1240003
d1$CSESid.3[d1$Es=="CAN_2011"]=1240001
d1$CSESid.4[d1$Es=="CAN_2011"]=1240004
d1$CSESid.5[d1$Es=="CAN_2011"]=1240005

d1$CSESid.1[d1$Es=="CAN_2015"]=1240001
d1$CSESid.2[d1$Es=="CAN_2015"]=1240002
d1$CSESid.3[d1$Es=="CAN_2015"]=1240003
d1$CSESid.4[d1$Es=="CAN_2015"]=1240004
d1$CSESid.5[d1$Es=="CAN_2015"]=1240005

d1$CSESid.1[d1$Es=="CZE_1996"]=2030002
d1$CSESid.2[d1$Es=="CZE_1996"]=2030001
d1$CSESid.3[d1$Es=="CZE_1996"]=2030004
d1$CSESid.4[d1$Es=="CZE_1996"]=2030003
d1$CSESid.5[d1$Es=="CZE_1996"]=2030018
d1$CSESid.6[d1$Es=="CZE_1996"]=2030011

d1$CSESid.1[d1$Es=="CZE_2006"]=2030002
d1$CSESid.2[d1$Es=="CZE_2006"]=2030001
d1$CSESid.3[d1$Es=="CZE_2006"]=2030004
d1$CSESid.4[d1$Es=="CZE_2006"]=2030003
d1$CSESid.5[d1$Es=="CZE_2006"]=2030005
d1$CSESid.6[d1$Es=="CZE_2006"]=2030074

d1$CSESid.1[d1$Es=="CZE_2010"]=2030001
d1$CSESid.2[d1$Es=="CZE_2010"]=2030002
d1$CSESid.3[d1$Es=="CZE_2010"]=2030010
d1$CSESid.4[d1$Es=="CZE_2010"]=2030004
d1$CSESid.5[d1$Es=="CZE_2010"]=2030014
d1$CSESid.6[d1$Es=="CZE_2010"]=2030003

d1$CSESid.1[d1$Es=="CZE_2013"]=2030001
d1$CSESid.2[d1$Es=="CZE_2013"]=2030006
d1$CSESid.3[d1$Es=="CZE_2013"]=2030004
d1$CSESid.4[d1$Es=="CZE_2013"]=2030010
d1$CSESid.5[d1$Es=="CZE_2013"]=2030002
d1$CSESid.6[d1$Es=="CZE_2013"]=2030013

d1$CSESid.1[d1$Es=="DEU_1998"]=2760004#DEU:scale up letters after removing CSU
d1$CSESid.2[d1$Es=="DEU_1998"]=2760001##note: no need to change 2760001 since it is already vote for CDU, or CSU in Bavaria
d1$CSESid.3[d1$Es=="DEU_1998"]=2760005
d1$CSESid.4[d1$Es=="DEU_1998"]=2760006
d1$CSESid.5[d1$Es=="DEU_1998"]=2760007

d1$CSESid.1[d1$Es=="DEU_2005"]=2760004#(comes from module 3, not module 2)
d1$CSESid.2[d1$Es=="DEU_2005"]=2760001
d1$CSESid.3[d1$Es=="DEU_2005"]=2760006
d1$CSESid.4[d1$Es=="DEU_2005"]=2760007
d1$CSESid.5[d1$Es=="DEU_2005"]=2760005
d1$CSESid.6[d1$Es=="DEU_2005"]=2760009

d1$CSESid.1[d1$Es=="DEU_2009"]=2760001
d1$CSESid.2[d1$Es=="DEU_2009"]=2760004
d1$CSESid.3[d1$Es=="DEU_2009"]=2760006
d1$CSESid.4[d1$Es=="DEU_2009"]=2760007
d1$CSESid.5[d1$Es=="DEU_2009"]=2760005

d1$CSESid.1[d1$Es=="DEU_2013"]=2760001
d1$CSESid.2[d1$Es=="DEU_2013"]=2760004
d1$CSESid.3[d1$Es=="DEU_2013"]=2760007
d1$CSESid.4[d1$Es=="DEU_2013"]=2760005
d1$CSESid.5[d1$Es=="DEU_2013"]=2760006

d1$CSESid.1[d1$Es=="DNK_1998"]=2080001
d1$CSESid.2[d1$Es=="DNK_1998"]=2080002
d1$CSESid.3[d1$Es=="DNK_1998"]=2080007
d1$CSESid.4[d1$Es=="DNK_1998"]=2080004
d1$CSESid.5[d1$Es=="DNK_1998"]=2080010
d1$CSESid.6[d1$Es=="DNK_1998"]=2080011

d1$CSESid.1[d1$Es=="DNK_2007"]=2080002
d1$CSESid.2[d1$Es=="DNK_2007"]=2080001
d1$CSESid.3[d1$Es=="DNK_2007"]=2080010
d1$CSESid.4[d1$Es=="DNK_2007"]=2080004
d1$CSESid.5[d1$Es=="DNK_2007"]=2080007
d1$CSESid.6[d1$Es=="DNK_2007"]=2080006

d1$CSESid.1[d1$Es=="ESP_1996"]=7240002
d1$CSESid.2[d1$Es=="ESP_1996"]=7240001
d1$CSESid.3[d1$Es=="ESP_1996"]=7240004
d1$CSESid.4[d1$Es=="ESP_1996"]=7240003
d1$CSESid.5[d1$Es=="ESP_1996"]=7240005

d1$CSESid.1[d1$Es=="ESP_2000"]=7240002
d1$CSESid.2[d1$Es=="ESP_2000"]=7240001
d1$CSESid.3[d1$Es=="ESP_2000"]=7240004
d1$CSESid.4[d1$Es=="ESP_2000"]=7240003
d1$CSESid.5[d1$Es=="ESP_2000"]=7240005

d1$CSESid.1[d1$Es=="ESP_2008"]=7240001
d1$CSESid.2[d1$Es=="ESP_2008"]=7240002
d1$CSESid.3[d1$Es=="ESP_2008"]=7240003
d1$CSESid.4[d1$Es=="ESP_2008"]=7240005
d1$CSESid.5[d1$Es=="ESP_2008"]=7240009
d1$CSESid.6[d1$Es=="ESP_2008"]=7240004

d1$CSESid.1[d1$Es=="FIN_2007"]=2460001
d1$CSESid.2[d1$Es=="FIN_2007"]=2460003
d1$CSESid.3[d1$Es=="FIN_2007"]=2460002
d1$CSESid.4[d1$Es=="FIN_2007"]=2460006
d1$CSESid.5[d1$Es=="FIN_2007"]=2460005
d1$CSESid.6[d1$Es=="FIN_2007"]=2460008

d1$CSESid.1[d1$Es=="FIN_2011"]=2460003
d1$CSESid.2[d1$Es=="FIN_2011"]=2460002
d1$CSESid.3[d1$Es=="FIN_2011"]=2460004
d1$CSESid.4[d1$Es=="FIN_2011"]=2460001
d1$CSESid.5[d1$Es=="FIN_2011"]=2460006
d1$CSESid.6[d1$Es=="FIN_2011"]=2460005

d1$CSESid.1[d1$Es=="FIN_2015"]=2460001
d1$CSESid.2[d1$Es=="FIN_2015"]=2460003
d1$CSESid.3[d1$Es=="FIN_2015"]=2460004
d1$CSESid.4[d1$Es=="FIN_2015"]=2460002
d1$CSESid.5[d1$Es=="FIN_2015"]=2460005
d1$CSESid.6[d1$Es=="FIN_2015"]=2460006

d1$CSESid.1[d1$Es=="FRA_2007"]=2500001
d1$CSESid.2[d1$Es=="FRA_2007"]=2500002
d1$CSESid.3[d1$Es=="FRA_2007"]=2500003
d1$CSESid.4[d1$Es=="FRA_2007"]=2500007
d1$CSESid.5[d1$Es=="FRA_2007"]=2500004
d1$CSESid.6[d1$Es=="FRA_2007"]=2500005

d1$CSESid.1[d1$Es=="GBR_1997"]=8260002
d1$CSESid.2[d1$Es=="GBR_1997"]=8260001
d1$CSESid.3[d1$Es=="GBR_1997"]=8260003
d1$CSESid.4[d1$Es=="GBR_1997"]=8260004
d1$CSESid.5[d1$Es=="GBR_1997"]=8260006

d1$CSESid.1[d1$Es=="GBR_2015"]=8260001
d1$CSESid.2[d1$Es=="GBR_2015"]=8260002
d1$CSESid.3[d1$Es=="GBR_2015"]=8260007
d1$CSESid.4[d1$Es=="GBR_2015"]=8260003
d1$CSESid.5[d1$Es=="GBR_2015"]=8260004
d1$CSESid.6[d1$Es=="GBR_2015"]=8260005

d1$CSESid.1[d1$Es=="GRC_2009"]=3000001
d1$CSESid.2[d1$Es=="GRC_2009"]=3000002
d1$CSESid.3[d1$Es=="GRC_2009"]=3000004
d1$CSESid.4[d1$Es=="GRC_2009"]=3000019
d1$CSESid.5[d1$Es=="GRC_2009"]=3000003
d1$CSESid.6[d1$Es=="GRC_2009"]=3000005

d1$CSESid.1[d1$Es=="GRC_2012"]=3000002
d1$CSESid.2[d1$Es=="GRC_2012"]=3000003
d1$CSESid.3[d1$Es=="GRC_2012"]=3000001
d1$CSESid.4[d1$Es=="GRC_2012"]=3000007
d1$CSESid.5[d1$Es=="GRC_2012"]=3000006
d1$CSESid.6[d1$Es=="GRC_2012"]=3000012

d1$CSESid.1[d1$Es=="GRC_2015"]=3000003
d1$CSESid.2[d1$Es=="GRC_2015"]=3000002
d1$CSESid.3[d1$Es=="GRC_2015"]=3000006
d1$CSESid.4[d1$Es=="GRC_2015"]=3000013
d1$CSESid.5[d1$Es=="GRC_2015"]=3000004
d1$CSESid.6[d1$Es=="GRC_2015"]=3000007

d1$CSESid.1[d1$Es=="HRV_2007"]=1910001
d1$CSESid.2[d1$Es=="HRV_2007"]=1910002
d1$CSESid.3[d1$Es=="HRV_2007"]=1910003
d1$CSESid.4[d1$Es=="HRV_2007"]=1910007
d1$CSESid.5[d1$Es=="HRV_2007"]=1910006
d1$CSESid.6[d1$Es=="HRV_2007"]=1910010

d1$CSESid.1[d1$Es=="IRL_2007"]=3720001
d1$CSESid.2[d1$Es=="IRL_2007"]=3720002
d1$CSESid.3[d1$Es=="IRL_2007"]=3720003
d1$CSESid.4[d1$Es=="IRL_2007"]=3720005
d1$CSESid.5[d1$Es=="IRL_2007"]=3720004
d1$CSESid.6[d1$Es=="IRL_2007"]=3720006

d1$CSESid.1[d1$Es=="IRL_2011"]=3720002#E not coded, substituted with F
d1$CSESid.2[d1$Es=="IRL_2011"]=3720003
d1$CSESid.3[d1$Es=="IRL_2011"]=3720001
d1$CSESid.4[d1$Es=="IRL_2011"]=3720004
d1$CSESid.5[d1$Es=="IRL_2011"]=3720005

d1$CSESid.1[d1$Es=="ISL_1999"]=3520001
d1$CSESid.2[d1$Es=="ISL_1999"]=3520002
d1$CSESid.3[d1$Es=="ISL_1999"]=3520003
d1$CSESid.4[d1$Es=="ISL_1999"]=3520005
d1$CSESid.5[d1$Es=="ISL_1999"]=3520004

d1$CSESid.1[d1$Es=="ISL_2007"]=3520001
d1$CSESid.2[d1$Es=="ISL_2007"]=3520002
d1$CSESid.3[d1$Es=="ISL_2007"]=3520005
d1$CSESid.4[d1$Es=="ISL_2007"]=3520003
d1$CSESid.5[d1$Es=="ISL_2007"]=3520004
d1$CSESid.6[d1$Es=="ISL_2007"]=3520013

d1$CSESid.1[d1$Es=="ISL_2009"]=3520002
d1$CSESid.2[d1$Es=="ISL_2009"]=3520001
d1$CSESid.3[d1$Es=="ISL_2009"]=3520005
d1$CSESid.4[d1$Es=="ISL_2009"]=3520003
d1$CSESid.5[d1$Es=="ISL_2009"]=3520009
d1$CSESid.6[d1$Es=="ISL_2009"]=3520004

d1$CSESid.1[d1$Es=="ISL_2013"]=3520001
d1$CSESid.2[d1$Es=="ISL_2013"]=3520003
d1$CSESid.3[d1$Es=="ISL_2013"]=3520002
d1$CSESid.4[d1$Es=="ISL_2013"]=3520005
d1$CSESid.5[d1$Es=="ISL_2013"]=3520007
d1$CSESid.6[d1$Es=="ISL_2013"]=3520008

d1$CSESid.1[d1$Es=="ISR_1996"]=3760002
d1$CSESid.2[d1$Es=="ISR_1996"]=3760001
d1$CSESid.3[d1$Es=="ISR_1996"]=3760003
d1$CSESid.4[d1$Es=="ISR_1996"]=3760014
d1$CSESid.5[d1$Es=="ISR_1996"]=3760004

d1$CSESid.1[d1$Es=="ISR_2006"]=3760016
d1$CSESid.2[d1$Es=="ISR_2006"]=3760002
d1$CSESid.3[d1$Es=="ISR_2006"]=3760003
d1$CSESid.4[d1$Es=="ISR_2006"]=3760001
d1$CSESid.5[d1$Es=="ISR_2006"]=3760021
d1$CSESid.6[d1$Es=="ISR_2006"]=3760050

d1$CSESid.1[d1$Es=="ISR_2013"]=3760011#F not coded, substituted with G
d1$CSESid.2[d1$Es=="ISR_2013"]=3760007
d1$CSESid.3[d1$Es=="ISR_2013"]=3760002
d1$CSESid.4[d1$Es=="ISR_2013"]=3760008
d1$CSESid.5[d1$Es=="ISR_2013"]=3760003
d1$CSESid.6[d1$Es=="ISR_2013"]=3760009

d1$CSESid.1[d1$Es=="LVA_2010"]=4280002
d1$CSESid.2[d1$Es=="LVA_2010"]=4280001
d1$CSESid.3[d1$Es=="LVA_2010"]=4280005
d1$CSESid.4[d1$Es=="LVA_2010"]=4280003
d1$CSESid.5[d1$Es=="LVA_2010"]=4280004
d1$CSESid.6[d1$Es=="LVA_2010"]=4280006

d1$CSESid.1[d1$Es=="LVA_2011"]=4280001
d1$CSESid.2[d1$Es=="LVA_2011"]=4280007
d1$CSESid.3[d1$Es=="LVA_2011"]=4280002
d1$CSESid.4[d1$Es=="LVA_2011"]=4280003
d1$CSESid.5[d1$Es=="LVA_2011"]=4280005
d1$CSESid.6[d1$Es=="LVA_2011"]=4280036

d1$CSESid.1[d1$Es=="LVA_2014"]=4280001
d1$CSESid.2[d1$Es=="LVA_2014"]=4280002
d1$CSESid.3[d1$Es=="LVA_2014"]=4280005
d1$CSESid.4[d1$Es=="LVA_2014"]=4280003
d1$CSESid.5[d1$Es=="LVA_2014"]=4280013
d1$CSESid.6[d1$Es=="LVA_2014"]=4280014

d1$CSESid.1[d1$Es=="MEX_2000"]=4840022
d1$CSESid.2[d1$Es=="MEX_2000"]=4840001
d1$CSESid.3[d1$Es=="MEX_2000"]=4840023

d1$CSESid.1[d1$Es=="MEX_2006"]=4840002
d1$CSESid.2[d1$Es=="MEX_2006"]=4840003
d1$CSESid.3[d1$Es=="MEX_2006"]=4840001

d1$CSESid.1[d1$Es=="MEX_2009"]=4840001
d1$CSESid.2[d1$Es=="MEX_2009"]=4840002
d1$CSESid.3[d1$Es=="MEX_2009"]=4840003

d1$CSESid.1[d1$Es=="MNE_2012"]=4990011
d1$CSESid.2[d1$Es=="MNE_2012"]=4990012
d1$CSESid.3[d1$Es=="MNE_2012"]=4990003
d1$CSESid.4[d1$Es=="MNE_2012"]=4990006
d1$CSESid.5[d1$Es=="MNE_2012"]=4990009
d1$CSESid.6[d1$Es=="MNE_2012"]=4990037

d1$CSESid.1[d1$Es=="NLD_1998"]=5280002
d1$CSESid.2[d1$Es=="NLD_1998"]=5280003
d1$CSESid.3[d1$Es=="NLD_1998"]=5280001
d1$CSESid.4[d1$Es=="NLD_1998"]=5280004
d1$CSESid.5[d1$Es=="NLD_1998"]=5280005
d1$CSESid.6[d1$Es=="NLD_1998"]=5280006

d1$CSESid.1[d1$Es=="NLD_2006"]=5280001
d1$CSESid.2[d1$Es=="NLD_2006"]=5280002
d1$CSESid.3[d1$Es=="NLD_2006"]=5280006
d1$CSESid.4[d1$Es=="NLD_2006"]=5280003
d1$CSESid.5[d1$Es=="NLD_2006"]=5280009
d1$CSESid.6[d1$Es=="NLD_2006"]=5280005

d1$CSESid.1[d1$Es=="NLD_2010"]=5280003
d1$CSESid.2[d1$Es=="NLD_2010"]=5280002
d1$CSESid.3[d1$Es=="NLD_2010"]=5280009
d1$CSESid.4[d1$Es=="NLD_2010"]=5280001
d1$CSESid.5[d1$Es=="NLD_2010"]=5280006
d1$CSESid.6[d1$Es=="NLD_2010"]=5280004

d1$CSESid.1[d1$Es=="NOR_1997"]=5780001
d1$CSESid.2[d1$Es=="NOR_1997"]=5780003
d1$CSESid.3[d1$Es=="NOR_1997"]=5780002
d1$CSESid.4[d1$Es=="NOR_1997"]=5780007
d1$CSESid.5[d1$Es=="NOR_1997"]=5780006
d1$CSESid.6[d1$Es=="NOR_1997"]=5780004

d1$CSESid.1[d1$Es=="NOR_2005"]=5780001
d1$CSESid.2[d1$Es=="NOR_2005"]=5780003
d1$CSESid.3[d1$Es=="NOR_2005"]=5780002
d1$CSESid.4[d1$Es=="NOR_2005"]=5780004
d1$CSESid.5[d1$Es=="NOR_2005"]=5780007
d1$CSESid.6[d1$Es=="NOR_2005"]=5780006

d1$CSESid.1[d1$Es=="NOR_2009"]=5780001
d1$CSESid.2[d1$Es=="NOR_2009"]=5780003
d1$CSESid.3[d1$Es=="NOR_2009"]=5780002
d1$CSESid.4[d1$Es=="NOR_2009"]=5780004
d1$CSESid.5[d1$Es=="NOR_2009"]=5780006
d1$CSESid.6[d1$Es=="NOR_2009"]=5780007

d1$CSESid.1[d1$Es=="NOR_2013"]=5780001
d1$CSESid.2[d1$Es=="NOR_2013"]=5780002
d1$CSESid.3[d1$Es=="NOR_2013"]=5780003
d1$CSESid.4[d1$Es=="NOR_2013"]=5780007
d1$CSESid.5[d1$Es=="NOR_2013"]=5780006
d1$CSESid.6[d1$Es=="NOR_2013"]=5780008

d1$CSESid.1[d1$Es=="NZL_1996"]=5540001
d1$CSESid.2[d1$Es=="NZL_1996"]=5540002
d1$CSESid.3[d1$Es=="NZL_1996"]=5540003
d1$CSESid.4[d1$Es=="NZL_1996"]=5540009
d1$CSESid.5[d1$Es=="NZL_1996"]=5540004
d1$CSESid.6[d1$Es=="NZL_1996"]=5540037

d1$CSESid.1[d1$Es=="NZL_2008"]=5540001
d1$CSESid.2[d1$Es=="NZL_2008"]=5540002
d1$CSESid.3[d1$Es=="NZL_2008"]=5540005
d1$CSESid.4[d1$Es=="NZL_2008"]=5540003
d1$CSESid.5[d1$Es=="NZL_2008"]=5540004
d1$CSESid.6[d1$Es=="NZL_2008"]=5540008

d1$CSESid.1[d1$Es=="NZL_2011"]=5540001
d1$CSESid.2[d1$Es=="NZL_2011"]=5540002
d1$CSESid.3[d1$Es=="NZL_2011"]=5540005
d1$CSESid.4[d1$Es=="NZL_2011"]=5540003
d1$CSESid.5[d1$Es=="NZL_2011"]=5540007
d1$CSESid.6[d1$Es=="NZL_2011"]=5540008

d1$CSESid.1[d1$Es=="NZL_2014"]=5540001
d1$CSESid.2[d1$Es=="NZL_2014"]=5540002
d1$CSESid.3[d1$Es=="NZL_2014"]=5540005
d1$CSESid.4[d1$Es=="NZL_2014"]=5540003
d1$CSESid.5[d1$Es=="NZL_2014"]=5540007
d1$CSESid.6[d1$Es=="NZL_2014"]=5540010

d1$CSESid.1[d1$Es=="PER_2001"]=6040008
d1$CSESid.2[d1$Es=="PER_2001"]=6040026
d1$CSESid.3[d1$Es=="PER_2001"]=6040052
d1$CSESid.4[d1$Es=="PER_2001"]=6040023
d1$CSESid.5[d1$Es=="PER_2001"]=6040030
d1$CSESid.6[d1$Es=="PER_2001"]=6040027

d1$CSESid.1[d1$Es=="PER_2011"]=6040050
d1$CSESid.2[d1$Es=="PER_2011"]=6040001
d1$CSESid.3[d1$Es=="PER_2011"]=6040008
d1$CSESid.4[d1$Es=="PER_2011"]=6040051
d1$CSESid.5[d1$Es=="PER_2011"]=6040015

d1$CSESid.1[d1$Es=="PER_2016"]=6040001
d1$CSESid.2[d1$Es=="PER_2016"]=6040002
d1$CSESid.3[d1$Es=="PER_2016"]=6040003
d1$CSESid.4[d1$Es=="PER_2016"]=6040004
d1$CSESid.5[d1$Es=="PER_2016"]=6040005
d1$CSESid.6[d1$Es=="PER_2016"]=6040006

d1$CSESid.1[d1$Es=="POL_1997"]=6160004
d1$CSESid.2[d1$Es=="POL_1997"]=6160009
d1$CSESid.3[d1$Es=="POL_1997"]=6160008
d1$CSESid.4[d1$Es=="POL_1997"]=6160003
d1$CSESid.5[d1$Es=="POL_1997"]=6160010
d1$CSESid.6[d1$Es=="POL_1997"]=6160007

d1$CSESid.1[d1$Es=="POL_2005"]=6160002#F not coded, substituted with H
d1$CSESid.2[d1$Es=="POL_2005"]=6160001
d1$CSESid.3[d1$Es=="POL_2005"]=6160024
d1$CSESid.4[d1$Es=="POL_2005"]=6160009
d1$CSESid.5[d1$Es=="POL_2005"]=6160014
d1$CSESid.6[d1$Es=="POL_2005"]=6160068

d1$CSESid.1[d1$Es=="POL_2007"]=6160001
d1$CSESid.2[d1$Es=="POL_2007"]=6160002
d1$CSESid.3[d1$Es=="POL_2007"]=6160020
d1$CSESid.4[d1$Es=="POL_2007"]=6160003

d1$CSESid.1[d1$Es=="POL_2011"]=6160001
d1$CSESid.2[d1$Es=="POL_2011"]=6160002
d1$CSESid.3[d1$Es=="POL_2011"]=6160017
d1$CSESid.4[d1$Es=="POL_2011"]=6160003
d1$CSESid.5[d1$Es=="POL_2011"]=6160009

d1$CSESid.1[d1$Es=="PRT_2002"]=6200001
d1$CSESid.2[d1$Es=="PRT_2002"]=6200002
d1$CSESid.3[d1$Es=="PRT_2002"]=6200003
d1$CSESid.4[d1$Es=="PRT_2002"]=6200004
d1$CSESid.5[d1$Es=="PRT_2002"]=6200006
d1$CSESid.6[d1$Es=="PRT_2002"]=6200007

d1$CSESid.1[d1$Es=="PRT_2009"]=6200002
d1$CSESid.2[d1$Es=="PRT_2009"]=6200001
d1$CSESid.3[d1$Es=="PRT_2009"]=6200003
d1$CSESid.4[d1$Es=="PRT_2009"]=6200006
d1$CSESid.5[d1$Es=="PRT_2009"]=6200004

d1$CSESid.1[d1$Es=="PRT_2015"]=6200011
d1$CSESid.2[d1$Es=="PRT_2015"]=6200002
d1$CSESid.3[d1$Es=="PRT_2015"]=6200006
d1$CSESid.4[d1$Es=="PRT_2015"]=6200004
d1$CSESid.5[d1$Es=="PRT_2015"]=6200020

d1$CSESid.1[d1$Es=="SRB_2012"]=6880001
d1$CSESid.2[d1$Es=="SRB_2012"]=6880002
d1$CSESid.3[d1$Es=="SRB_2012"]=6880004
d1$CSESid.4[d1$Es=="SRB_2012"]=6880005
d1$CSESid.5[d1$Es=="SRB_2012"]=6880006
d1$CSESid.6[d1$Es=="SRB_2012"]=6880007

d1$CSESid.1[d1$Es=="SVK_2010"]=7030001
d1$CSESid.2[d1$Es=="SVK_2010"]=7030003
d1$CSESid.3[d1$Es=="SVK_2010"]=7030002
d1$CSESid.4[d1$Es=="SVK_2010"]=7030004
d1$CSESid.5[d1$Es=="SVK_2010"]=7030008
d1$CSESid.6[d1$Es=="SVK_2010"]=7030006

d1$CSESid.1[d1$Es=="SVK_2016"]=7030001
d1$CSESid.2[d1$Es=="SVK_2016"]=7030002
d1$CSESid.3[d1$Es=="SVK_2016"]=7030009
d1$CSESid.4[d1$Es=="SVK_2016"]=7030006
d1$CSESid.5[d1$Es=="SVK_2016"]=7030012
d1$CSESid.6[d1$Es=="SVK_2016"]=7030014

d1$CSESid.1[d1$Es=="SVN_1996"]=7050004#F not coded, substituted with G
d1$CSESid.2[d1$Es=="SVN_1996"]=7050006
d1$CSESid.3[d1$Es=="SVN_1996"]=7050001
d1$CSESid.4[d1$Es=="SVN_1996"]=7050017
d1$CSESid.5[d1$Es=="SVN_1996"]=7050002
d1$CSESid.6[d1$Es=="SVN_1996"]=7050008

d1$CSESid.1[d1$Es=="SVN_2008"]=7050002
d1$CSESid.2[d1$Es=="SVN_2008"]=7050001
d1$CSESid.3[d1$Es=="SVN_2008"]=7050020
d1$CSESid.4[d1$Es=="SVN_2008"]=7050003
d1$CSESid.5[d1$Es=="SVN_2008"]=7050008
d1$CSESid.6[d1$Es=="SVN_2008"]=7050019

d1$CSESid.1[d1$Es=="SVN_2011"]=7050016
d1$CSESid.2[d1$Es=="SVN_2011"]=7050001
d1$CSESid.3[d1$Es=="SVN_2011"]=7050002
d1$CSESid.4[d1$Es=="SVN_2011"]=7050009
d1$CSESid.5[d1$Es=="SVN_2011"]=7050003
d1$CSESid.6[d1$Es=="SVN_2011"]=7050006

d1$CSESid.1[d1$Es=="SWE_1998"]=7520001
d1$CSESid.2[d1$Es=="SWE_1998"]=7520002
d1$CSESid.3[d1$Es=="SWE_1998"]=7520007
d1$CSESid.4[d1$Es=="SWE_1998"]=7520004
d1$CSESid.5[d1$Es=="SWE_1998"]=7520006
d1$CSESid.6[d1$Es=="SWE_1998"]=7520003

d1$CSESid.1[d1$Es=="SWE_2006"]=7520001
d1$CSESid.2[d1$Es=="SWE_2006"]=7520002
d1$CSESid.3[d1$Es=="SWE_2006"]=7520006
d1$CSESid.4[d1$Es=="SWE_2006"]=7520003
d1$CSESid.5[d1$Es=="SWE_2006"]=7520004
d1$CSESid.6[d1$Es=="SWE_2006"]=7520007

d1$CSESid.1[d1$Es=="SWE_2014"]=7520001
d1$CSESid.2[d1$Es=="SWE_2014"]=7520002
d1$CSESid.3[d1$Es=="SWE_2014"]=7520008
d1$CSESid.4[d1$Es=="SWE_2014"]=7520005
d1$CSESid.5[d1$Es=="SWE_2014"]=7520006
d1$CSESid.6[d1$Es=="SWE_2014"]=7520007

d1$CSESid.1[d1$Es=="THA_2007"]=7640014
d1$CSESid.2[d1$Es=="THA_2007"]=7640002
d1$CSESid.3[d1$Es=="THA_2007"]=7640021
d1$CSESid.4[d1$Es=="THA_2007"]=7640008
d1$CSESid.5[d1$Es=="THA_2007"]=7640060
d1$CSESid.6[d1$Es=="THA_2007"]=7640024

d1$CSESid.1[d1$Es=="THA_2011"]=7640001
d1$CSESid.2[d1$Es=="THA_2011"]=7640002
d1$CSESid.3[d1$Es=="THA_2011"]=7640003
d1$CSESid.4[d1$Es=="THA_2011"]=7640004
d1$CSESid.5[d1$Es=="THA_2011"]=7640005
d1$CSESid.6[d1$Es=="THA_2011"]=7640006

d1$CSESid.1[d1$Es=="TWN_1996"]=1580002
d1$CSESid.2[d1$Es=="TWN_1996"]=1580001
d1$CSESid.3[d1$Es=="TWN_1996"]=1580004

d1$CSESid.1[d1$Es=="ZAF_2009"]=7100001
d1$CSESid.2[d1$Es=="ZAF_2009"]=7100002
d1$CSESid.3[d1$Es=="ZAF_2009"]=7100008
d1$CSESid.4[d1$Es=="ZAF_2009"]=7100004
d1$CSESid.5[d1$Es=="ZAF_2009"]=7100022
d1$CSESid.6[d1$Es=="ZAF_2009"]=7100006

d1$CSESid.1[d1$Es=="ZAF_2014"]=7100001
d1$CSESid.2[d1$Es=="ZAF_2014"]=7100002
d1$CSESid.3[d1$Es=="ZAF_2014"]=7100003
d1$CSESid.4[d1$Es=="ZAF_2014"]=7100004
d1$CSESid.5[d1$Es=="ZAF_2014"]=7100005
d1$CSESid.6[d1$Es=="ZAF_2014"]=7100006


#d1 party identification (Pi):----

d1$Pi=NA

#assign Pi:
d1$Pi[d1$Es=="ARG_2015"&d1$IMD3005_3==0320001]=1#
d1$Pi[d1$Es=="ARG_2015"&d1$IMD3005_3==0320002]=2
d1$Pi[d1$Es=="ARG_2015"&d1$IMD3005_3==0320003]=3
d1$Pi[d1$Es=="ARG_2015"&d1$IMD3005_3==0320004]=4
d1$Pi[d1$Es=="ARG_2015"&d1$IMD3005_3==0320005]=5
d1$Pi[d1$Es=="ARG_2015"&d1$IMD3005_3==0320006]=6

d1$Pi[d1$Es=="AUS_1996"&d1$IMD3005_3==0360003]=1#
d1$Pi[d1$Es=="AUS_1996"&d1$IMD3005_3==0360001]=2
d1$Pi[d1$Es=="AUS_1996"&d1$IMD3005_3==0360002]=3
d1$Pi[d1$Es=="AUS_1996"&d1$IMD3005_3==0360004]=4

d1$Pi[d1$Es=="AUS_2004"&d1$IMD3005_3==0360001]=1#scaled up since E not coded
d1$Pi[d1$Es=="AUS_2004"&d1$IMD3005_3==0360003]=2
d1$Pi[d1$Es=="AUS_2004"&d1$IMD3005_3==0360005]=3
d1$Pi[d1$Es=="AUS_2004"&d1$IMD3005_3==0360002]=4
d1$Pi[d1$Es=="AUS_2004"&d1$IMD3005_3==0360004]=5
d1$Pi[d1$Es=="AUS_2004"&d1$IMD3005_3==0360008]=6

d1$Pi[d1$Es=="AUS_2007"&d1$IMD3005_3==0360003]=1#
d1$Pi[d1$Es=="AUS_2007"&d1$IMD3005_3==0360001]=2
d1$Pi[d1$Es=="AUS_2007"&d1$IMD3005_3==0360002]=3
d1$Pi[d1$Es=="AUS_2007"&d1$IMD3005_3==0360005]=4

d1$Pi[d1$Es=="AUS_2013"&d1$IMD3005_3==0360001]=1#scaled up E since D not coded
d1$Pi[d1$Es=="AUS_2013"&d1$IMD3005_3==0360003]=2
d1$Pi[d1$Es=="AUS_2013"&d1$IMD3005_3==0360005]=3
d1$Pi[d1$Es=="AUS_2013"&d1$IMD3005_3==0360002]=4

d1$Pi[d1$Es=="AUT_2013"&d1$IMD3005_3==0400001]=1#
d1$Pi[d1$Es=="AUT_2013"&d1$IMD3005_3==0400002]=2
d1$Pi[d1$Es=="AUT_2013"&d1$IMD3005_3==0400003]=3
d1$Pi[d1$Es=="AUT_2013"&d1$IMD3005_3==0400005]=4
d1$Pi[d1$Es=="AUT_2013"&d1$IMD3005_3==0400004]=5
d1$Pi[d1$Es=="AUT_2013"&d1$IMD3005_3==0400010]=6

d1$Pi[d1$Es=="BGR_2014"&d1$IMD3005_3==1000001]=1#
d1$Pi[d1$Es=="BGR_2014"&d1$IMD3005_3==1000020]=2
d1$Pi[d1$Es=="BGR_2014"&d1$IMD3005_3==1000003]=3
d1$Pi[d1$Es=="BGR_2014"&d1$IMD3005_3==1000004]=4
d1$Pi[d1$Es=="BGR_2014"&d1$IMD3005_3==1000005]=5
d1$Pi[d1$Es=="BGR_2014"&d1$IMD3005_3==1000006]=6

d1$Pi[d1$Es=="CAN_1997"&d1$IMD3005_3==1240001]=1#
d1$Pi[d1$Es=="CAN_1997"&d1$IMD3005_3==1240007]=2
d1$Pi[d1$Es=="CAN_1997"&d1$IMD3005_3==1240006]=3
d1$Pi[d1$Es=="CAN_1997"&d1$IMD3005_3==1240003]=4
d1$Pi[d1$Es=="CAN_1997"&d1$IMD3005_3==1240004]=5

d1$Pi[d1$Es=="CAN_2008"&d1$IMD3005_3==1240002]=1#
d1$Pi[d1$Es=="CAN_2008"&d1$IMD3005_3==1240001]=2
d1$Pi[d1$Es=="CAN_2008"&d1$IMD3005_3==1240004]=3
d1$Pi[d1$Es=="CAN_2008"&d1$IMD3005_3==1240003]=4
d1$Pi[d1$Es=="CAN_2008"&d1$IMD3005_3==1240005]=5

d1$Pi[d1$Es=="CAN_2011"&d1$IMD3005_3==1240002]=1#
d1$Pi[d1$Es=="CAN_2011"&d1$IMD3005_3==1240003]=2
d1$Pi[d1$Es=="CAN_2011"&d1$IMD3005_3==1240001]=3
d1$Pi[d1$Es=="CAN_2011"&d1$IMD3005_3==1240004]=4
d1$Pi[d1$Es=="CAN_2011"&d1$IMD3005_3==1240005]=5

d1$Pi[d1$Es=="CAN_2015"&d1$IMD3005_3==1240001]=1#
d1$Pi[d1$Es=="CAN_2015"&d1$IMD3005_3==1240002]=2
d1$Pi[d1$Es=="CAN_2015"&d1$IMD3005_3==1240003]=3
d1$Pi[d1$Es=="CAN_2015"&d1$IMD3005_3==1240004]=4
d1$Pi[d1$Es=="CAN_2015"&d1$IMD3005_3==1240005]=5

d1$Pi[d1$Es=="CZE_1996"&d1$IMD3005_3==2030002]=1#
d1$Pi[d1$Es=="CZE_1996"&d1$IMD3005_3==2030001]=2
d1$Pi[d1$Es=="CZE_1996"&d1$IMD3005_3==2030004]=3
d1$Pi[d1$Es=="CZE_1996"&d1$IMD3005_3==2030003]=4
d1$Pi[d1$Es=="CZE_1996"&d1$IMD3005_3==2030018]=5
d1$Pi[d1$Es=="CZE_1996"&d1$IMD3005_3==2030011]=6

d1$Pi[d1$Es=="CZE_2006"&d1$IMD3005_3==2030002]=1#
d1$Pi[d1$Es=="CZE_2006"&d1$IMD3005_3==2030001]=2
d1$Pi[d1$Es=="CZE_2006"&d1$IMD3005_3==2030004]=3
d1$Pi[d1$Es=="CZE_2006"&d1$IMD3005_3==2030003]=4
d1$Pi[d1$Es=="CZE_2006"&d1$IMD3005_3==2030005]=5
d1$Pi[d1$Es=="CZE_2006"&d1$IMD3005_3==2030074]=6

d1$Pi[d1$Es=="CZE_2010"&d1$IMD3005_3==2030001]=1#
d1$Pi[d1$Es=="CZE_2010"&d1$IMD3005_3==2030002]=2
d1$Pi[d1$Es=="CZE_2010"&d1$IMD3005_3==2030010]=3
d1$Pi[d1$Es=="CZE_2010"&d1$IMD3005_3==2030004]=4
d1$Pi[d1$Es=="CZE_2010"&d1$IMD3005_3==2030014]=5
d1$Pi[d1$Es=="CZE_2010"&d1$IMD3005_3==2030003]=6

d1$Pi[d1$Es=="CZE_2013"&d1$IMD3005_3==2030001]=1#
d1$Pi[d1$Es=="CZE_2013"&d1$IMD3005_3==2030006]=2
d1$Pi[d1$Es=="CZE_2013"&d1$IMD3005_3==2030004]=3
d1$Pi[d1$Es=="CZE_2013"&d1$IMD3005_3==2030010]=4
d1$Pi[d1$Es=="CZE_2013"&d1$IMD3005_3==2030002]=5
d1$Pi[d1$Es=="CZE_2013"&d1$IMD3005_3==2030013]=6

d1$Pi[d1$Es=="DEU_1998"&d1$IMD3005_3==2760004]=1#DEU:scale up letters after removing CSU
d1$Pi[d1$Es=="DEU_1998"&d1$IMD3005_3==2760001]=2##note: no need to change 2760001 since it is already vote for CDU, or CSU in Bavaria
d1$Pi[d1$Es=="DEU_1998"&d1$IMD3005_3==2760005]=3
d1$Pi[d1$Es=="DEU_1998"&d1$IMD3005_3==2760006]=4
d1$Pi[d1$Es=="DEU_1998"&d1$IMD3005_3==2760007]=5

d1$Pi[d1$Es=="DEU_2005"&d1$IMD3005_3==2760004]=1#(comes from module 3, not module 2)
d1$Pi[d1$Es=="DEU_2005"&d1$IMD3005_3==2760001]=2
d1$Pi[d1$Es=="DEU_2005"&d1$IMD3005_3==2760006]=3
d1$Pi[d1$Es=="DEU_2005"&d1$IMD3005_3==2760007]=4
d1$Pi[d1$Es=="DEU_2005"&d1$IMD3005_3==2760005]=5
d1$Pi[d1$Es=="DEU_2005"&d1$IMD3005_3==2760009]=6

d1$Pi[d1$Es=="DEU_2009"&d1$IMD3005_3==2760001]=1#
d1$Pi[d1$Es=="DEU_2009"&d1$IMD3005_3==2760004]=2
d1$Pi[d1$Es=="DEU_2009"&d1$IMD3005_3==2760006]=3
d1$Pi[d1$Es=="DEU_2009"&d1$IMD3005_3==2760007]=4
d1$Pi[d1$Es=="DEU_2009"&d1$IMD3005_3==2760005]=5

d1$Pi[d1$Es=="DEU_2013"&d1$IMD3005_3==2760001]=1#
d1$Pi[d1$Es=="DEU_2013"&d1$IMD3005_3==2760004]=2
d1$Pi[d1$Es=="DEU_2013"&d1$IMD3005_3==2760007]=3
d1$Pi[d1$Es=="DEU_2013"&d1$IMD3005_3==2760005]=4
d1$Pi[d1$Es=="DEU_2013"&d1$IMD3005_3==2760006]=5

d1$Pi[d1$Es=="DNK_1998"&d1$IMD3005_3==2080001]=1#
d1$Pi[d1$Es=="DNK_1998"&d1$IMD3005_3==2080002]=2
d1$Pi[d1$Es=="DNK_1998"&d1$IMD3005_3==2080007]=3
d1$Pi[d1$Es=="DNK_1998"&d1$IMD3005_3==2080004]=4
d1$Pi[d1$Es=="DNK_1998"&d1$IMD3005_3==2080010]=5
d1$Pi[d1$Es=="DNK_1998"&d1$IMD3005_3==2080011]=6

d1$Pi[d1$Es=="DNK_2007"&d1$IMD3005_3==2080002]=1#
d1$Pi[d1$Es=="DNK_2007"&d1$IMD3005_3==2080001]=2
d1$Pi[d1$Es=="DNK_2007"&d1$IMD3005_3==2080010]=3
d1$Pi[d1$Es=="DNK_2007"&d1$IMD3005_3==2080004]=4
d1$Pi[d1$Es=="DNK_2007"&d1$IMD3005_3==2080007]=5
d1$Pi[d1$Es=="DNK_2007"&d1$IMD3005_3==2080006]=6

d1$Pi[d1$Es=="ESP_1996"&d1$IMD3005_3==7240002]=1#
d1$Pi[d1$Es=="ESP_1996"&d1$IMD3005_3==7240001]=2
d1$Pi[d1$Es=="ESP_1996"&d1$IMD3005_3==7240004]=3
d1$Pi[d1$Es=="ESP_1996"&d1$IMD3005_3==7240003]=4
d1$Pi[d1$Es=="ESP_1996"&d1$IMD3005_3==7240005]=5

d1$Pi[d1$Es=="ESP_2000"&d1$IMD3005_3==7240002]=1#
d1$Pi[d1$Es=="ESP_2000"&d1$IMD3005_3==7240001]=2
d1$Pi[d1$Es=="ESP_2000"&d1$IMD3005_3==7240004]=3
d1$Pi[d1$Es=="ESP_2000"&d1$IMD3005_3==7240003]=4
d1$Pi[d1$Es=="ESP_2000"&d1$IMD3005_3==7240005]=5

d1$Pi[d1$Es=="ESP_2008"&d1$IMD3005_3==7240001]=1#
d1$Pi[d1$Es=="ESP_2008"&d1$IMD3005_3==7240002]=2
d1$Pi[d1$Es=="ESP_2008"&d1$IMD3005_3==7240003]=3
d1$Pi[d1$Es=="ESP_2008"&d1$IMD3005_3==7240005]=4
d1$Pi[d1$Es=="ESP_2008"&d1$IMD3005_3==7240009]=5
d1$Pi[d1$Es=="ESP_2008"&d1$IMD3005_3==7240004]=6

d1$Pi[d1$Es=="FIN_2007"&d1$IMD3005_3==2460001]=1#
d1$Pi[d1$Es=="FIN_2007"&d1$IMD3005_3==2460003]=2
d1$Pi[d1$Es=="FIN_2007"&d1$IMD3005_3==2460002]=3
d1$Pi[d1$Es=="FIN_2007"&d1$IMD3005_3==2460006]=4
d1$Pi[d1$Es=="FIN_2007"&d1$IMD3005_3==2460005]=5
d1$Pi[d1$Es=="FIN_2007"&d1$IMD3005_3==2460008]=6

d1$Pi[d1$Es=="FIN_2011"&d1$IMD3005_3==2460003]=1#
d1$Pi[d1$Es=="FIN_2011"&d1$IMD3005_3==2460002]=2
d1$Pi[d1$Es=="FIN_2011"&d1$IMD3005_3==2460004]=3
d1$Pi[d1$Es=="FIN_2011"&d1$IMD3005_3==2460001]=4
d1$Pi[d1$Es=="FIN_2011"&d1$IMD3005_3==2460006]=5
d1$Pi[d1$Es=="FIN_2011"&d1$IMD3005_3==2460005]=6

d1$Pi[d1$Es=="FIN_2015"&d1$IMD3005_3==2460001]=1#
d1$Pi[d1$Es=="FIN_2015"&d1$IMD3005_3==2460003]=2
d1$Pi[d1$Es=="FIN_2015"&d1$IMD3005_3==2460004]=3
d1$Pi[d1$Es=="FIN_2015"&d1$IMD3005_3==2460002]=4
d1$Pi[d1$Es=="FIN_2015"&d1$IMD3005_3==2460005]=5
d1$Pi[d1$Es=="FIN_2015"&d1$IMD3005_3==2460006]=6

d1$Pi[d1$Es=="FRA_2007"&d1$IMD3005_3==2500001]=1#
d1$Pi[d1$Es=="FRA_2007"&d1$IMD3005_3==2500002]=2
d1$Pi[d1$Es=="FRA_2007"&d1$IMD3005_3==2500003]=3
d1$Pi[d1$Es=="FRA_2007"&d1$IMD3005_3==2500007]=4
d1$Pi[d1$Es=="FRA_2007"&d1$IMD3005_3==2500004]=5
d1$Pi[d1$Es=="FRA_2007"&d1$IMD3005_3==2500005]=6

d1$Pi[d1$Es=="GBR_1997"&d1$IMD3005_3==8260002]=1#
d1$Pi[d1$Es=="GBR_1997"&d1$IMD3005_3==8260001]=2
d1$Pi[d1$Es=="GBR_1997"&d1$IMD3005_3==8260003]=3
d1$Pi[d1$Es=="GBR_1997"&d1$IMD3005_3==8260004]=4
d1$Pi[d1$Es=="GBR_1997"&d1$IMD3005_3==8260006]=5

d1$Pi[d1$Es=="GBR_2015"&d1$IMD3005_3==8260001]=1#
d1$Pi[d1$Es=="GBR_2015"&d1$IMD3005_3==8260002]=2
d1$Pi[d1$Es=="GBR_2015"&d1$IMD3005_3==8260007]=3
d1$Pi[d1$Es=="GBR_2015"&d1$IMD3005_3==8260003]=4
d1$Pi[d1$Es=="GBR_2015"&d1$IMD3005_3==8260004]=5
d1$Pi[d1$Es=="GBR_2015"&d1$IMD3005_3==8260005]=6

d1$Pi[d1$Es=="GRC_2009"&d1$IMD3005_3==3000001]=1#
d1$Pi[d1$Es=="GRC_2009"&d1$IMD3005_3==3000002]=2
d1$Pi[d1$Es=="GRC_2009"&d1$IMD3005_3==3000004]=3
d1$Pi[d1$Es=="GRC_2009"&d1$IMD3005_3==3000019]=4
d1$Pi[d1$Es=="GRC_2009"&d1$IMD3005_3==3000003]=5
d1$Pi[d1$Es=="GRC_2009"&d1$IMD3005_3==3000005]=6

d1$Pi[d1$Es=="GRC_2012"&d1$IMD3005_3==3000002]=1#
d1$Pi[d1$Es=="GRC_2012"&d1$IMD3005_3==3000003]=2
d1$Pi[d1$Es=="GRC_2012"&d1$IMD3005_3==3000001]=3
d1$Pi[d1$Es=="GRC_2012"&d1$IMD3005_3==3000007]=4
d1$Pi[d1$Es=="GRC_2012"&d1$IMD3005_3==3000006]=5
d1$Pi[d1$Es=="GRC_2012"&d1$IMD3005_3==3000012]=6

d1$Pi[d1$Es=="GRC_2015"&d1$IMD3005_3==3000003]=1#
d1$Pi[d1$Es=="GRC_2015"&d1$IMD3005_3==3000002]=2
d1$Pi[d1$Es=="GRC_2015"&d1$IMD3005_3==3000006]=3
d1$Pi[d1$Es=="GRC_2015"&d1$IMD3005_3==3000013]=4
d1$Pi[d1$Es=="GRC_2015"&d1$IMD3005_3==3000004]=5
d1$Pi[d1$Es=="GRC_2015"&d1$IMD3005_3==3000007]=6

d1$Pi[d1$Es=="HRV_2007"&d1$IMD3005_3==1910001]=1#
d1$Pi[d1$Es=="HRV_2007"&d1$IMD3005_3==1910002]=2
d1$Pi[d1$Es=="HRV_2007"&d1$IMD3005_3==1910003]=3
d1$Pi[d1$Es=="HRV_2007"&d1$IMD3005_3==1910007]=4
d1$Pi[d1$Es=="HRV_2007"&d1$IMD3005_3==1910006]=5
d1$Pi[d1$Es=="HRV_2007"&d1$IMD3005_3==1910010]=6

d1$Pi[d1$Es=="IRL_2007"&d1$IMD3005_3==3720001]=1#
d1$Pi[d1$Es=="IRL_2007"&d1$IMD3005_3==3720002]=2
d1$Pi[d1$Es=="IRL_2007"&d1$IMD3005_3==3720003]=3
d1$Pi[d1$Es=="IRL_2007"&d1$IMD3005_3==3720005]=4
d1$Pi[d1$Es=="IRL_2007"&d1$IMD3005_3==3720004]=5
d1$Pi[d1$Es=="IRL_2007"&d1$IMD3005_3==3720006]=6

d1$Pi[d1$Es=="IRL_2011"&d1$IMD3005_3==3720002]=1#E not coded, substituted with F
d1$Pi[d1$Es=="IRL_2011"&d1$IMD3005_3==3720003]=2
d1$Pi[d1$Es=="IRL_2011"&d1$IMD3005_3==3720001]=3
d1$Pi[d1$Es=="IRL_2011"&d1$IMD3005_3==3720004]=4
d1$Pi[d1$Es=="IRL_2011"&d1$IMD3005_3==3720005]=5

d1$Pi[d1$Es=="ISL_1999"&d1$IMD3005_3==3520001]=1#
d1$Pi[d1$Es=="ISL_1999"&d1$IMD3005_3==3520002]=2
d1$Pi[d1$Es=="ISL_1999"&d1$IMD3005_3==3520003]=3
d1$Pi[d1$Es=="ISL_1999"&d1$IMD3005_3==3520005]=4
d1$Pi[d1$Es=="ISL_1999"&d1$IMD3005_3==3520004]=5

d1$Pi[d1$Es=="ISL_2007"&d1$IMD3005_3==3520001]=1#
d1$Pi[d1$Es=="ISL_2007"&d1$IMD3005_3==3520002]=2
d1$Pi[d1$Es=="ISL_2007"&d1$IMD3005_3==3520005]=3
d1$Pi[d1$Es=="ISL_2007"&d1$IMD3005_3==3520003]=4
d1$Pi[d1$Es=="ISL_2007"&d1$IMD3005_3==3520004]=5
d1$Pi[d1$Es=="ISL_2007"&d1$IMD3005_3==3520013]=6

d1$Pi[d1$Es=="ISL_2009"&d1$IMD3005_3==3520002]=1#
d1$Pi[d1$Es=="ISL_2009"&d1$IMD3005_3==3520001]=2
d1$Pi[d1$Es=="ISL_2009"&d1$IMD3005_3==3520005]=3
d1$Pi[d1$Es=="ISL_2009"&d1$IMD3005_3==3520003]=4
d1$Pi[d1$Es=="ISL_2009"&d1$IMD3005_3==3520009]=5
d1$Pi[d1$Es=="ISL_2009"&d1$IMD3005_3==3520004]=6

d1$Pi[d1$Es=="ISL_2013"&d1$IMD3005_3==3520001]=1#
d1$Pi[d1$Es=="ISL_2013"&d1$IMD3005_3==3520003]=2
d1$Pi[d1$Es=="ISL_2013"&d1$IMD3005_3==3520002]=3
d1$Pi[d1$Es=="ISL_2013"&d1$IMD3005_3==3520005]=4
d1$Pi[d1$Es=="ISL_2013"&d1$IMD3005_3==3520007]=5
d1$Pi[d1$Es=="ISL_2013"&d1$IMD3005_3==3520008]=6

d1$Pi[d1$Es=="ISR_1996"&d1$IMD3005_3==3760002]=1#
d1$Pi[d1$Es=="ISR_1996"&d1$IMD3005_3==3760001]=2
d1$Pi[d1$Es=="ISR_1996"&d1$IMD3005_3==3760003]=3
d1$Pi[d1$Es=="ISR_1996"&d1$IMD3005_3==3760014]=4
d1$Pi[d1$Es=="ISR_1996"&d1$IMD3005_3==3760004]=5

d1$Pi[d1$Es=="ISR_2006"&d1$IMD3005_3==3760016]=1#
d1$Pi[d1$Es=="ISR_2006"&d1$IMD3005_3==3760002]=2
d1$Pi[d1$Es=="ISR_2006"&d1$IMD3005_3==3760003]=3
d1$Pi[d1$Es=="ISR_2006"&d1$IMD3005_3==3760001]=4
d1$Pi[d1$Es=="ISR_2006"&d1$IMD3005_3==3760021]=5
d1$Pi[d1$Es=="ISR_2006"&d1$IMD3005_3==3760050]=6

d1$Pi[d1$Es=="ISR_2013"&d1$IMD3005_3==3760011]=1#F not coded, substituted with G
d1$Pi[d1$Es=="ISR_2013"&d1$IMD3005_3==3760007]=2
d1$Pi[d1$Es=="ISR_2013"&d1$IMD3005_3==3760002]=3
d1$Pi[d1$Es=="ISR_2013"&d1$IMD3005_3==3760008]=4
d1$Pi[d1$Es=="ISR_2013"&d1$IMD3005_3==3760003]=5
d1$Pi[d1$Es=="ISR_2013"&d1$IMD3005_3==3760009]=6

d1$Pi[d1$Es=="LVA_2010"&d1$IMD3005_3==4280002]=1#
d1$Pi[d1$Es=="LVA_2010"&d1$IMD3005_3==4280001]=2
d1$Pi[d1$Es=="LVA_2010"&d1$IMD3005_3==4280005]=3
d1$Pi[d1$Es=="LVA_2010"&d1$IMD3005_3==4280003]=4
d1$Pi[d1$Es=="LVA_2010"&d1$IMD3005_3==4280004]=5
d1$Pi[d1$Es=="LVA_2010"&d1$IMD3005_3==4280006]=6

d1$Pi[d1$Es=="LVA_2011"&d1$IMD3005_3==4280001]=1#
d1$Pi[d1$Es=="LVA_2011"&d1$IMD3005_3==4280007]=2
d1$Pi[d1$Es=="LVA_2011"&d1$IMD3005_3==4280002]=3
d1$Pi[d1$Es=="LVA_2011"&d1$IMD3005_3==4280003]=4
d1$Pi[d1$Es=="LVA_2011"&d1$IMD3005_3==4280005]=5
d1$Pi[d1$Es=="LVA_2011"&d1$IMD3005_3==4280036]=6

d1$Pi[d1$Es=="LVA_2014"&d1$IMD3005_3==4280001]=1#
d1$Pi[d1$Es=="LVA_2014"&d1$IMD3005_3==4280002]=2
d1$Pi[d1$Es=="LVA_2014"&d1$IMD3005_3==4280005]=3
d1$Pi[d1$Es=="LVA_2014"&d1$IMD3005_3==4280003]=4
d1$Pi[d1$Es=="LVA_2014"&d1$IMD3005_3==4280013]=5
d1$Pi[d1$Es=="LVA_2014"&d1$IMD3005_3==4280014]=6

d1$Pi[d1$Es=="MEX_2000"&d1$IMD3005_3==4840022]=1#
d1$Pi[d1$Es=="MEX_2000"&d1$IMD3005_3==4840001]=2
d1$Pi[d1$Es=="MEX_2000"&d1$IMD3005_3==4840023]=3

d1$Pi[d1$Es=="MEX_2006"&d1$IMD3005_3==4840002]=1#
d1$Pi[d1$Es=="MEX_2006"&d1$IMD3005_3==4840003]=2
d1$Pi[d1$Es=="MEX_2006"&d1$IMD3005_3==4840001]=3

d1$Pi[d1$Es=="MEX_2009"&d1$IMD3005_3==4840001]=1#
d1$Pi[d1$Es=="MEX_2009"&d1$IMD3005_3==4840002]=2
d1$Pi[d1$Es=="MEX_2009"&d1$IMD3005_3==4840003]=3

d1$Pi[d1$Es=="MNE_2012"&d1$IMD3005_3==4990011]=1#
d1$Pi[d1$Es=="MNE_2012"&d1$IMD3005_3==4990012]=2
d1$Pi[d1$Es=="MNE_2012"&d1$IMD3005_3==4990003]=3
d1$Pi[d1$Es=="MNE_2012"&d1$IMD3005_3==4990006]=4
d1$Pi[d1$Es=="MNE_2012"&d1$IMD3005_3==4990009]=5
d1$Pi[d1$Es=="MNE_2012"&d1$IMD3005_3==4990037]=6

d1$Pi[d1$Es=="NLD_1998"&d1$IMD3005_3==5280002]=1#
d1$Pi[d1$Es=="NLD_1998"&d1$IMD3005_3==5280003]=2
d1$Pi[d1$Es=="NLD_1998"&d1$IMD3005_3==5280001]=3
d1$Pi[d1$Es=="NLD_1998"&d1$IMD3005_3==5280004]=4
d1$Pi[d1$Es=="NLD_1998"&d1$IMD3005_3==5280005]=5
d1$Pi[d1$Es=="NLD_1998"&d1$IMD3005_3==5280006]=6

d1$Pi[d1$Es=="NLD_2006"&d1$IMD3005_3==5280001]=1#
d1$Pi[d1$Es=="NLD_2006"&d1$IMD3005_3==5280002]=2
d1$Pi[d1$Es=="NLD_2006"&d1$IMD3005_3==5280006]=3
d1$Pi[d1$Es=="NLD_2006"&d1$IMD3005_3==5280003]=4
d1$Pi[d1$Es=="NLD_2006"&d1$IMD3005_3==5280009]=5
d1$Pi[d1$Es=="NLD_2006"&d1$IMD3005_3==5280005]=6

d1$Pi[d1$Es=="NLD_2010"&d1$IMD3005_3==5280003]=1#
d1$Pi[d1$Es=="NLD_2010"&d1$IMD3005_3==5280002]=2
d1$Pi[d1$Es=="NLD_2010"&d1$IMD3005_3==5280009]=3
d1$Pi[d1$Es=="NLD_2010"&d1$IMD3005_3==5280001]=4
d1$Pi[d1$Es=="NLD_2010"&d1$IMD3005_3==5280006]=5
d1$Pi[d1$Es=="NLD_2010"&d1$IMD3005_3==5280004]=6

d1$Pi[d1$Es=="NOR_1997"&d1$IMD3005_3==5780001]=1#
d1$Pi[d1$Es=="NOR_1997"&d1$IMD3005_3==5780003]=2
d1$Pi[d1$Es=="NOR_1997"&d1$IMD3005_3==5780002]=3
d1$Pi[d1$Es=="NOR_1997"&d1$IMD3005_3==5780007]=4
d1$Pi[d1$Es=="NOR_1997"&d1$IMD3005_3==5780006]=5
d1$Pi[d1$Es=="NOR_1997"&d1$IMD3005_3==5780004]=6

d1$Pi[d1$Es=="NOR_2005"&d1$IMD3005_3==5780001]=1#
d1$Pi[d1$Es=="NOR_2005"&d1$IMD3005_3==5780003]=2
d1$Pi[d1$Es=="NOR_2005"&d1$IMD3005_3==5780002]=3
d1$Pi[d1$Es=="NOR_2005"&d1$IMD3005_3==5780004]=4
d1$Pi[d1$Es=="NOR_2005"&d1$IMD3005_3==5780007]=5
d1$Pi[d1$Es=="NOR_2005"&d1$IMD3005_3==5780006]=6

d1$Pi[d1$Es=="NOR_2009"&d1$IMD3005_3==5780001]=1#
d1$Pi[d1$Es=="NOR_2009"&d1$IMD3005_3==5780003]=2
d1$Pi[d1$Es=="NOR_2009"&d1$IMD3005_3==5780002]=3
d1$Pi[d1$Es=="NOR_2009"&d1$IMD3005_3==5780004]=4
d1$Pi[d1$Es=="NOR_2009"&d1$IMD3005_3==5780006]=5
d1$Pi[d1$Es=="NOR_2009"&d1$IMD3005_3==5780007]=6

d1$Pi[d1$Es=="NOR_2013"&d1$IMD3005_3==5780001]=1#
d1$Pi[d1$Es=="NOR_2013"&d1$IMD3005_3==5780002]=2
d1$Pi[d1$Es=="NOR_2013"&d1$IMD3005_3==5780003]=3
d1$Pi[d1$Es=="NOR_2013"&d1$IMD3005_3==5780007]=4
d1$Pi[d1$Es=="NOR_2013"&d1$IMD3005_3==5780006]=5
d1$Pi[d1$Es=="NOR_2013"&d1$IMD3005_3==5780008]=6

d1$Pi[d1$Es=="NZL_1996"&d1$IMD3005_3==5540001]=1#
d1$Pi[d1$Es=="NZL_1996"&d1$IMD3005_3==5540002]=2
d1$Pi[d1$Es=="NZL_1996"&d1$IMD3005_3==5540003]=3
d1$Pi[d1$Es=="NZL_1996"&d1$IMD3005_3==5540009]=4
d1$Pi[d1$Es=="NZL_1996"&d1$IMD3005_3==5540004]=5
d1$Pi[d1$Es=="NZL_1996"&d1$IMD3005_3==5540037]=6

d1$Pi[d1$Es=="NZL_2008"&d1$IMD3005_3==5540001]=1#
d1$Pi[d1$Es=="NZL_2008"&d1$IMD3005_3==5540002]=2
d1$Pi[d1$Es=="NZL_2008"&d1$IMD3005_3==5540005]=3
d1$Pi[d1$Es=="NZL_2008"&d1$IMD3005_3==5540003]=4
d1$Pi[d1$Es=="NZL_2008"&d1$IMD3005_3==5540004]=5
d1$Pi[d1$Es=="NZL_2008"&d1$IMD3005_3==5540008]=6

d1$Pi[d1$Es=="NZL_2011"&d1$IMD3005_3==5540001]=1#
d1$Pi[d1$Es=="NZL_2011"&d1$IMD3005_3==5540002]=2
d1$Pi[d1$Es=="NZL_2011"&d1$IMD3005_3==5540005]=3
d1$Pi[d1$Es=="NZL_2011"&d1$IMD3005_3==5540003]=4
d1$Pi[d1$Es=="NZL_2011"&d1$IMD3005_3==5540007]=5
d1$Pi[d1$Es=="NZL_2011"&d1$IMD3005_3==5540008]=6

d1$Pi[d1$Es=="NZL_2014"&d1$IMD3005_3==5540001]=1#
d1$Pi[d1$Es=="NZL_2014"&d1$IMD3005_3==5540002]=2
d1$Pi[d1$Es=="NZL_2014"&d1$IMD3005_3==5540005]=3
d1$Pi[d1$Es=="NZL_2014"&d1$IMD3005_3==5540003]=4
d1$Pi[d1$Es=="NZL_2014"&d1$IMD3005_3==5540007]=5
d1$Pi[d1$Es=="NZL_2014"&d1$IMD3005_3==5540010]=6

d1$Pi[d1$Es=="PER_2001"&d1$IMD3005_3==6040008]=1#
d1$Pi[d1$Es=="PER_2001"&d1$IMD3005_3==6040026]=2
d1$Pi[d1$Es=="PER_2001"&d1$IMD3005_3==6040052]=3
d1$Pi[d1$Es=="PER_2001"&d1$IMD3005_3==6040023]=4
d1$Pi[d1$Es=="PER_2001"&d1$IMD3005_3==6040030]=5
d1$Pi[d1$Es=="PER_2001"&d1$IMD3005_3==6040027]=6

d1$Pi[d1$Es=="PER_2011"&d1$IMD3005_3==6040050]=1#
d1$Pi[d1$Es=="PER_2011"&d1$IMD3005_3==6040001]=2
d1$Pi[d1$Es=="PER_2011"&d1$IMD3005_3==6040008]=3
d1$Pi[d1$Es=="PER_2011"&d1$IMD3005_3==6040051]=4
d1$Pi[d1$Es=="PER_2011"&d1$IMD3005_3==6040015]=5

d1$Pi[d1$Es=="PER_2016"&d1$IMD3005_3==6040001]=1#
d1$Pi[d1$Es=="PER_2016"&d1$IMD3005_3==6040002]=2
d1$Pi[d1$Es=="PER_2016"&d1$IMD3005_3==6040003]=3
d1$Pi[d1$Es=="PER_2016"&d1$IMD3005_3==6040004]=4
d1$Pi[d1$Es=="PER_2016"&d1$IMD3005_3==6040005]=5
d1$Pi[d1$Es=="PER_2016"&d1$IMD3005_3==6040006]=6

d1$Pi[d1$Es=="POL_1997"&d1$IMD3005_3==6160004]=1#
d1$Pi[d1$Es=="POL_1997"&d1$IMD3005_3==6160009]=2
d1$Pi[d1$Es=="POL_1997"&d1$IMD3005_3==6160008]=3
d1$Pi[d1$Es=="POL_1997"&d1$IMD3005_3==6160003]=4
d1$Pi[d1$Es=="POL_1997"&d1$IMD3005_3==6160010]=5
d1$Pi[d1$Es=="POL_1997"&d1$IMD3005_3==6160007]=6

d1$Pi[d1$Es=="POL_2005"&d1$IMD3005_3==6160002]=1#F not coded, substituted with H
d1$Pi[d1$Es=="POL_2005"&d1$IMD3005_3==6160001]=2
d1$Pi[d1$Es=="POL_2005"&d1$IMD3005_3==6160024]=3
d1$Pi[d1$Es=="POL_2005"&d1$IMD3005_3==6160009]=4
d1$Pi[d1$Es=="POL_2005"&d1$IMD3005_3==6160014]=5
d1$Pi[d1$Es=="POL_2005"&d1$IMD3005_3==6160068]=6

d1$Pi[d1$Es=="POL_2007"&d1$IMD3005_3==6160001]=1#
d1$Pi[d1$Es=="POL_2007"&d1$IMD3005_3==6160002]=2
d1$Pi[d1$Es=="POL_2007"&d1$IMD3005_3==6160020]=3
d1$Pi[d1$Es=="POL_2007"&d1$IMD3005_3==6160003]=4

d1$Pi[d1$Es=="POL_2011"&d1$IMD3005_3==6160001]=1#
d1$Pi[d1$Es=="POL_2011"&d1$IMD3005_3==6160002]=2
d1$Pi[d1$Es=="POL_2011"&d1$IMD3005_3==6160017]=3
d1$Pi[d1$Es=="POL_2011"&d1$IMD3005_3==6160003]=4
d1$Pi[d1$Es=="POL_2011"&d1$IMD3005_3==6160009]=5

d1$Pi[d1$Es=="PRT_2002"&d1$IMD3005_3==6200001]=1#
d1$Pi[d1$Es=="PRT_2002"&d1$IMD3005_3==6200002]=2
d1$Pi[d1$Es=="PRT_2002"&d1$IMD3005_3==6200003]=3
d1$Pi[d1$Es=="PRT_2002"&d1$IMD3005_3==6200004]=4
d1$Pi[d1$Es=="PRT_2002"&d1$IMD3005_3==6200006]=5
d1$Pi[d1$Es=="PRT_2002"&d1$IMD3005_3==6200007]=6

d1$Pi[d1$Es=="PRT_2009"&d1$IMD3005_3==6200002]=1#
d1$Pi[d1$Es=="PRT_2009"&d1$IMD3005_3==6200001]=2
d1$Pi[d1$Es=="PRT_2009"&d1$IMD3005_3==6200003]=3
d1$Pi[d1$Es=="PRT_2009"&d1$IMD3005_3==6200006]=4
d1$Pi[d1$Es=="PRT_2009"&d1$IMD3005_3==6200004]=5

d1$Pi[d1$Es=="PRT_2015"&d1$IMD3005_3==6200011]=1#
d1$Pi[d1$Es=="PRT_2015"&d1$IMD3005_3==6200002]=2
d1$Pi[d1$Es=="PRT_2015"&d1$IMD3005_3==6200006]=3
d1$Pi[d1$Es=="PRT_2015"&d1$IMD3005_3==6200004]=4
d1$Pi[d1$Es=="PRT_2015"&d1$IMD3005_3==6200020]=5

d1$Pi[d1$Es=="SRB_2012"&d1$IMD3005_3==6880001]=1#
d1$Pi[d1$Es=="SRB_2012"&d1$IMD3005_3==6880002]=2
d1$Pi[d1$Es=="SRB_2012"&d1$IMD3005_3==6880004]=3
d1$Pi[d1$Es=="SRB_2012"&d1$IMD3005_3==6880005]=4
d1$Pi[d1$Es=="SRB_2012"&d1$IMD3005_3==6880006]=5
d1$Pi[d1$Es=="SRB_2012"&d1$IMD3005_3==6880007]=6

d1$Pi[d1$Es=="SVK_2010"&d1$IMD3005_3==7030001]=1#
d1$Pi[d1$Es=="SVK_2010"&d1$IMD3005_3==7030003]=2
d1$Pi[d1$Es=="SVK_2010"&d1$IMD3005_3==7030002]=3
d1$Pi[d1$Es=="SVK_2010"&d1$IMD3005_3==7030004]=4
d1$Pi[d1$Es=="SVK_2010"&d1$IMD3005_3==7030008]=5
d1$Pi[d1$Es=="SVK_2010"&d1$IMD3005_3==7030006]=6

d1$Pi[d1$Es=="SVK_2016"&d1$IMD3005_3==7030001]=1#
d1$Pi[d1$Es=="SVK_2016"&d1$IMD3005_3==7030002]=2
d1$Pi[d1$Es=="SVK_2016"&d1$IMD3005_3==7030009]=3
d1$Pi[d1$Es=="SVK_2016"&d1$IMD3005_3==7030006]=4
d1$Pi[d1$Es=="SVK_2016"&d1$IMD3005_3==7030012]=5
d1$Pi[d1$Es=="SVK_2016"&d1$IMD3005_3==7030014]=6

d1$Pi[d1$Es=="SVN_1996"&d1$IMD3005_3==7050004]=1#F not coded, substituted with G
d1$Pi[d1$Es=="SVN_1996"&d1$IMD3005_3==7050006]=2
d1$Pi[d1$Es=="SVN_1996"&d1$IMD3005_3==7050001]=3
d1$Pi[d1$Es=="SVN_1996"&d1$IMD3005_3==7050017]=4
d1$Pi[d1$Es=="SVN_1996"&d1$IMD3005_3==7050002]=5
d1$Pi[d1$Es=="SVN_1996"&d1$IMD3005_3==7050008]=6

d1$Pi[d1$Es=="SVN_2008"&d1$IMD3005_3==7050002]=1#
d1$Pi[d1$Es=="SVN_2008"&d1$IMD3005_3==7050001]=2
d1$Pi[d1$Es=="SVN_2008"&d1$IMD3005_3==7050020]=3
d1$Pi[d1$Es=="SVN_2008"&d1$IMD3005_3==7050003]=4
d1$Pi[d1$Es=="SVN_2008"&d1$IMD3005_3==7050008]=5
d1$Pi[d1$Es=="SVN_2008"&d1$IMD3005_3==7050019]=6

d1$Pi[d1$Es=="SVN_2011"&d1$IMD3005_3==7050016]=1#
d1$Pi[d1$Es=="SVN_2011"&d1$IMD3005_3==7050001]=2
d1$Pi[d1$Es=="SVN_2011"&d1$IMD3005_3==7050002]=3
d1$Pi[d1$Es=="SVN_2011"&d1$IMD3005_3==7050009]=4
d1$Pi[d1$Es=="SVN_2011"&d1$IMD3005_3==7050003]=5
d1$Pi[d1$Es=="SVN_2011"&d1$IMD3005_3==7050006]=6

d1$Pi[d1$Es=="SWE_1998"&d1$IMD3005_3==7520001]=1#
d1$Pi[d1$Es=="SWE_1998"&d1$IMD3005_3==7520002]=2
d1$Pi[d1$Es=="SWE_1998"&d1$IMD3005_3==7520007]=3
d1$Pi[d1$Es=="SWE_1998"&d1$IMD3005_3==7520004]=4
d1$Pi[d1$Es=="SWE_1998"&d1$IMD3005_3==7520006]=5
d1$Pi[d1$Es=="SWE_1998"&d1$IMD3005_3==7520003]=6

d1$Pi[d1$Es=="SWE_2006"&d1$IMD3005_3==7520001]=1#
d1$Pi[d1$Es=="SWE_2006"&d1$IMD3005_3==7520002]=2
d1$Pi[d1$Es=="SWE_2006"&d1$IMD3005_3==7520006]=3
d1$Pi[d1$Es=="SWE_2006"&d1$IMD3005_3==7520003]=4
d1$Pi[d1$Es=="SWE_2006"&d1$IMD3005_3==7520004]=5
d1$Pi[d1$Es=="SWE_2006"&d1$IMD3005_3==7520007]=6

d1$Pi[d1$Es=="SWE_2014"&d1$IMD3005_3==7520001]=1#
d1$Pi[d1$Es=="SWE_2014"&d1$IMD3005_3==7520002]=2
d1$Pi[d1$Es=="SWE_2014"&d1$IMD3005_3==7520008]=3
d1$Pi[d1$Es=="SWE_2014"&d1$IMD3005_3==7520005]=4
d1$Pi[d1$Es=="SWE_2014"&d1$IMD3005_3==7520006]=5
d1$Pi[d1$Es=="SWE_2014"&d1$IMD3005_3==7520007]=6

d1$Pi[d1$Es=="THA_2007"&d1$IMD3005_3==7640014]=1
d1$Pi[d1$Es=="THA_2007"&d1$IMD3005_3==7640002]=2
d1$Pi[d1$Es=="THA_2007"&d1$IMD3005_3==7640021]=3
d1$Pi[d1$Es=="THA_2007"&d1$IMD3005_3==7640008]=4
d1$Pi[d1$Es=="THA_2007"&d1$IMD3005_3==7640060]=5
d1$Pi[d1$Es=="THA_2007"&d1$IMD3005_3==7640024]=6

d1$Pi[d1$Es=="THA_2011"&d1$IMD3005_3==7640001]=1
d1$Pi[d1$Es=="THA_2011"&d1$IMD3005_3==7640002]=2
d1$Pi[d1$Es=="THA_2011"&d1$IMD3005_3==7640003]=3
d1$Pi[d1$Es=="THA_2011"&d1$IMD3005_3==7640004]=4
d1$Pi[d1$Es=="THA_2011"&d1$IMD3005_3==7640005]=5
d1$Pi[d1$Es=="THA_2011"&d1$IMD3005_3==7640006]=6

d1$Pi[d1$Es=="TWN_1996"&d1$IMD3005_3==1580002]=1#
d1$Pi[d1$Es=="TWN_1996"&d1$IMD3005_3==1580001]=2
d1$Pi[d1$Es=="TWN_1996"&d1$IMD3005_3==1580004]=3

d1$Pi[d1$Es=="ZAF_2009"&d1$IMD3005_3==7100001]=1#
d1$Pi[d1$Es=="ZAF_2009"&d1$IMD3005_3==7100002]=2
d1$Pi[d1$Es=="ZAF_2009"&d1$IMD3005_3==7100008]=3
d1$Pi[d1$Es=="ZAF_2009"&d1$IMD3005_3==7100004]=4
d1$Pi[d1$Es=="ZAF_2009"&d1$IMD3005_3==7100022]=5
d1$Pi[d1$Es=="ZAF_2009"&d1$IMD3005_3==7100006]=6

d1$Pi[d1$Es=="ZAF_2014"&d1$IMD3005_3==7100001]=1#
d1$Pi[d1$Es=="ZAF_2014"&d1$IMD3005_3==7100002]=2
d1$Pi[d1$Es=="ZAF_2014"&d1$IMD3005_3==7100003]=3
d1$Pi[d1$Es=="ZAF_2014"&d1$IMD3005_3==7100004]=4
d1$Pi[d1$Es=="ZAF_2014"&d1$IMD3005_3==7100005]=5
d1$Pi[d1$Es=="ZAF_2014"&d1$IMD3005_3==7100006]=6


#D5----
#Election variables:
colnames(d5)[colnames(d5)=="E1001"]="Ed"#Election Dataset
colnames(d5)[colnames(d5)=="E1004"]="Es"#Election Study
d5$Es[d5$Es=="GRC_2015"]="GRC_2016"#temporarily modifying 2015s year to code all the variables with
#unique election-year, it will be modified before the analysis
d5$Ec=substr(d5$Es, 1, 3)#Election Country
colnames(d5)[colnames(d5)=="E1008"]="Ey"#Election Year
d5$Ey[d5$Es=="GRC_2016"]=2016#temporarily modifying 2015s year to code all the variables with
#unique election-year, it will be modified before the analysis
colnames(d5)[colnames(d5)=="E1005"]="Rid"#Respondent id
colnames(d5)[colnames(d5)=="E2001_Y"]="Rage"#Respondent Age (E2001_Y is year of birth, get age detracting it to Ey)
d5$Rage[d5$Rage>2020]=NA
d5$Rage=d5$Ey-d5$Rage
d5$Rage[d5$Rage>115]=NA
d5$Rwd=1#Respondent Weight Demographic:
d5$Rwd[d5$Es=="AUS_2019"]=d5$E1010_1[d5$Es=="AUS_2019"]
d5$Rwd[d5$Es=="AUT_2017"]=d5$E1010_2[d5$Es=="AUT_2017"]
d5$Rwd[d5$Es=="DEU_2017"]=d5$E1010_2[d5$Es=="DEU_2017"]
d5$Rwd[d5$Es=="GRC_2016"]=d5$E1010_2[d5$Es=="GRC_2016"]
d5$Rwd[d5$Es=="IRL_2016"]=d5$E1010_2[d5$Es=="IRL_2016"]
d5$Rwd[d5$Es=="ITA_2018"]=d5$E1010_2[d5$Es=="ITA_2018"]
d5$Rwd[d5$Es=="NOR_2017"]=d5$E1010_2[d5$Es=="NOR_2017"]
d5$Rwd[d5$Es=="NZL_2017"]=d5$E1010_2[d5$Es=="NZL_2017"]
#Vote and party/leader features:
colnames(d5)[colnames(d5)=="E3012_LH"]="Vt"#Vote Turnout (0=non voted; 1=voted) CURRENT LOWER HOUSE ELECTION: DID RESPONDENT CAST A BALLOT
d5$Vt[d5$Vt>1]=NA#coding missing as NAs
colnames(d5)[colnames(d5)=="E3013_LH_PL"]="Vn"#Vote Numerical CURRENT LOWER HOUSE ELECTION: VOTE CHOICE - PARTY LIST
d5$Vn[d5$Es=="IRL_2016"]=d5$E3013_LH_DC[d5$Es=="IRL_2016"]#assigning LH_DC to Vn for relevant elections
d5$Vn[d5$Es=="AUS_2019"]=d5$E3013_LH_DC[d5$Es=="AUS_2019"]
d5$Va=NA#Vote Alphabetical: code vote from numerical (000001-999987) to alphabetical modified (1:6):
#d5 assign alphabetical modified codes instead than numerical codes----

d5$Va[d5$Es=="AUS_2019"&d5$Vn==036001]=1##
d5$Va[d5$Es=="AUS_2019"&d5$Vn==036002]=2
d5$Va[d5$Es=="AUS_2019"&d5$Vn==036003]=3
d5$Va[d5$Es=="AUS_2019"&d5$Vn==036004]=4
d5$Va[d5$Es=="AUS_2019"&d5$Vn==036005]=5
d5$Va[d5$Es=="AUS_2019"&d5$Vn==036006]=6

d5$Va[d5$Es=="AUT_2017"&d5$Vn==040001]=1#
d5$Va[d5$Es=="AUT_2017"&d5$Vn==040002]=2
d5$Va[d5$Es=="AUT_2017"&d5$Vn==040003]=3
d5$Va[d5$Es=="AUT_2017"&d5$Vn==040004]=4
d5$Va[d5$Es=="AUT_2017"&d5$Vn==040005]=5
d5$Va[d5$Es=="AUT_2017"&d5$Vn==040006]=6

d5$Va[d5$Es=="DEU_2017"&d5$Vn==276009]=1#vote for either CDU,CSU (same comments than for 2760001)
d5$Va[d5$Es=="DEU_2017"&d5$Vn==276002]=2
d5$Va[d5$Es=="DEU_2017"&d5$Vn==276003]=3
d5$Va[d5$Es=="DEU_2017"&d5$Vn==276004]=4
d5$Va[d5$Es=="DEU_2017"&d5$Vn==276005]=5
d5$Va[d5$Es=="DEU_2017"&d5$Vn==276006]=6

d5$Va[d5$Es=="GRC_2016"&d5$Vn==300001]=1#
d5$Va[d5$Es=="GRC_2016"&d5$Vn==300002]=2
d5$Va[d5$Es=="GRC_2016"&d5$Vn==300003]=3
d5$Va[d5$Es=="GRC_2016"&d5$Vn==300004]=4
d5$Va[d5$Es=="GRC_2016"&d5$Vn==300005]=5
d5$Va[d5$Es=="GRC_2016"&d5$Vn==300006]=6

d5$Va[d5$Es=="IRL_2016"&d5$Vn==372001]=1#
d5$Va[d5$Es=="IRL_2016"&d5$Vn==372002]=2
d5$Va[d5$Es=="IRL_2016"&d5$Vn==372003]=3
d5$Va[d5$Es=="IRL_2016"&d5$Vn==372004]=4
d5$Va[d5$Es=="IRL_2016"&d5$Vn==372005]=5
d5$Va[d5$Es=="IRL_2016"&d5$Vn==372006]=6

d5$Va[d5$Es=="ISL_2016"&d5$Vn==352001]=1##
d5$Va[d5$Es=="ISL_2016"&d5$Vn==352002]=2
d5$Va[d5$Es=="ISL_2016"&d5$Vn==352003]=3
d5$Va[d5$Es=="ISL_2016"&d5$Vn==352004]=4
d5$Va[d5$Es=="ISL_2016"&d5$Vn==352005]=5
d5$Va[d5$Es=="ISL_2016"&d5$Vn==352006]=6

d5$Va[d5$Es=="ISL_2017"&d5$Vn==352101]=1##
d5$Va[d5$Es=="ISL_2017"&d5$Vn==352102]=2
d5$Va[d5$Es=="ISL_2017"&d5$Vn==352103]=3
d5$Va[d5$Es=="ISL_2017"&d5$Vn==352104]=4
d5$Va[d5$Es=="ISL_2017"&d5$Vn==352105]=5
d5$Va[d5$Es=="ISL_2017"&d5$Vn==352106]=6

d5$Va[d5$Es=="ITA_2018"&d5$Vn==380001]=1#Movimento 5 Stelle
d5$Va[d5$Es=="ITA_2018"&d5$Vn==380002]=2#Partito Democratico (PD)
d5$Va[d5$Es=="ITA_2018"&d5$Vn==380003]=3#Lega (LN)
d5$Va[d5$Es=="ITA_2018"&d5$Vn==380004]=4
d5$Va[d5$Es=="ITA_2018"&d5$Vn==380005]=5
d5$Va[d5$Es=="ITA_2018"&d5$Vn==380006]=6

d5$Va[d5$Es=="NOR_2017"&d5$Vn==578001]=1##
d5$Va[d5$Es=="NOR_2017"&d5$Vn==578002]=2
d5$Va[d5$Es=="NOR_2017"&d5$Vn==578003]=3
d5$Va[d5$Es=="NOR_2017"&d5$Vn==578004]=4
d5$Va[d5$Es=="NOR_2017"&d5$Vn==578005]=5
d5$Va[d5$Es=="NOR_2017"&d5$Vn==578006]=6

d5$Va[d5$Es=="NZL_2017"&d5$Vn==554001]=1##
d5$Va[d5$Es=="NZL_2017"&d5$Vn==554002]=2
d5$Va[d5$Es=="NZL_2017"&d5$Vn==554003]=3
d5$Va[d5$Es=="NZL_2017"&d5$Vn==554004]=4
d5$Va[d5$Es=="NZL_2017"&d5$Vn==554005]=5
d5$Va[d5$Es=="NZL_2017"&d5$Vn==554006]=6

colnames(d5)[colnames(d5) %in% paste0("E3017_", LETTERS[1:6])]=paste0("Pl.", 1:6)#party likability [0,10] A-F into 1:6----
d5[,names(d5) %in% paste0("Pl.", 1:6)][d5[,names(d5) %in% paste0("Pl.", 1:6)]>10]=NA#coding missing as NAs

colnames(d5)[colnames(d5) %in% paste0("E3018_", LETTERS[1:6])]=paste0("Ll.", 1:6)#leader likability [0,10] A-F into 1:6
d5[,names(d5) %in% paste0("Ll.", 1:6)][d5[,names(d5) %in% paste0("Ll.", 1:6)]>10]=NA#coding missing as NAs

colnames(d5)[colnames(d5) %in% paste0("E5002_", LETTERS[1:6])]=paste0("Ps.", 1:6)#party seats A-F into 1:6 PERCENT SEATS - LOWER HOUSE
d5[,names(d5) %in% paste0("Ps.", 1:6)][d5[,names(d5) %in% paste0("Ps.", 1:6)]>100]=NA#coding missing as NAs

colnames(d5)[colnames(d5) %in% paste0("E5001_", LETTERS[1:6])]=paste0("Pv.", 1:6)#party vote A-F into 1:6 PERCENT VOTE - LOWER HOUSE
d5[,names(d5) %in% paste0("Pv.", 1:6)][d5[,names(d5) %in% paste0("Pv.", 1:6)]>100]=NA#coding missing as NAs

#EEeconS Election Economy Subjective:
colnames(d5)[colnames(d5)=="E3011"]="EeconS"#Election economy state Subjective (no data from CSES#3!!!)
d5$EeconS[d5$EeconS>5]=NA#coding missing as NAs

#recode EeconS from CSES scale
# 1. GOTTEN BETTER c(1,2,3,4,5)
# 5. GOTTEN WORSE 
# to TEV scale, noting that the TEV used the conversion scale below for same elections:
# 0. Negative evaluation c(0, 0.25, 0.5, 0.75, 1)
# 1. Positive evaluation
d5$EeconS[d5$EeconS==1]=1
d5$EeconS[d5$EeconS==2]=0.75
d5$EeconS[d5$EeconS==3]=0.5
d5$EeconS[d5$EeconS==4]=0.25
d5$EeconS[d5$EeconS==5]=0

#Elections with EeconS:
table(d5$Es,d5$EeconS)
d5$Eecon=1

#LRR=Left Right Respondent (voter LR self):
d5$LRR=d5$E3020
d5$LRR[d5$LRR>10]=NA
#recoding to 0-10 scale: already 0-10 scale
#LRP.i=Left Right Parties:
d5$LRP.1=d5$E3019_A#LEFT-RIGHT - PARTY A
d5$LRP.2=d5$E3019_B#LEFT-RIGHT - PARTY B
d5$LRP.3=d5$E3019_C#LEFT-RIGHT - PARTY C
d5$LRP.4=d5$E3019_D#LEFT-RIGHT - PARTY D
d5$LRP.5=d5$E3019_E#LEFT-RIGHT - PARTY E
d5$LRP.6=d5$E3019_F#LEFT-RIGHT - PARTY F
# d5$LRP.7=as.numeric(NA)
d5$LRP.1[d5$LRP.1>10]=NA
d5$LRP.2[d5$LRP.2>10]=NA
d5$LRP.3[d5$LRP.3>10]=NA
d5$LRP.4[d5$LRP.4>10]=NA
d5$LRP.5[d5$LRP.5>10]=NA
d5$LRP.6[d5$LRP.6>10]=NA
#recoding to 0-10 scale: already 0-10 scale
#Elections with LRR and LRP:
d5$LR=1

#Seats to vote ratios:
d5$P1SVrat=d5$Ps.1/d5$Pv.1#for most voted party
d5$P2SVrat=d5$Ps.2/d5$Pv.2#for second most voted party
d5$P1P2rat=d5$P1SVrat/d5$P2SVrat#ratio of rations first/second party
#Gallagher Disproportionality measure:
d5$GDM=NA
d5$GDM[d5$Es=="AUT_2017"]=3.73
d5$GDM[d5$Es=="DEU_2017"]=1.95#
d5$GDM[d5$Es=="GRC_2016"]=9.60#Sep
d5$GDM[d5$Es=="IRL_2016"]=5.62
d5$GDM[d5$Es=="ITA_2018"]=5.70

d5$GDM[d5$Es=="AUS_2019"]=12.35#
d5$GDM[d5$Es=="ISL_2016"]=4.46#
d5$GDM[d5$Es=="ISL_2017"]=1.91#
d5$GDM[d5$Es=="NOR_2017"]=3.01#
d5$GDM[d5$Es=="NZL_2017"]=2.73#

#d5 assign CSES numerical ID to alternatives (CSESid.i):----
d5$CSESid.1=NA
d5$CSESid.2=NA
d5$CSESid.3=NA
d5$CSESid.4=NA
d5$CSESid.5=NA
d5$CSESid.6=NA
# d5$CSESid.7=NA

d5$CSESid.1[d5$Es=="AUS_2019"]=036001
d5$CSESid.2[d5$Es=="AUS_2019"]=036002
d5$CSESid.3[d5$Es=="AUS_2019"]=036003
d5$CSESid.4[d5$Es=="AUS_2019"]=036004
d5$CSESid.5[d5$Es=="AUS_2019"]=036005
d5$CSESid.6[d5$Es=="AUS_2019"]=036006

d5$CSESid.1[d5$Es=="AUT_2017"]=040001
d5$CSESid.2[d5$Es=="AUT_2017"]=040002
d5$CSESid.3[d5$Es=="AUT_2017"]=040003
d5$CSESid.4[d5$Es=="AUT_2017"]=040004
d5$CSESid.5[d5$Es=="AUT_2017"]=040005
d5$CSESid.6[d5$Es=="AUT_2017"]=040006

d5$CSESid.1[d5$Es=="DEU_2017"]=276009#vote for either CDU,CSU (same comments than for 2760001)
d5$CSESid.2[d5$Es=="DEU_2017"]=276002
d5$CSESid.3[d5$Es=="DEU_2017"]=276003
d5$CSESid.4[d5$Es=="DEU_2017"]=276004
d5$CSESid.5[d5$Es=="DEU_2017"]=276005
d5$CSESid.6[d5$Es=="DEU_2017"]=276006

d5$CSESid.1[d5$Es=="GRC_2016"]=300001
d5$CSESid.2[d5$Es=="GRC_2016"]=300002
d5$CSESid.3[d5$Es=="GRC_2016"]=300003
d5$CSESid.4[d5$Es=="GRC_2016"]=300004
d5$CSESid.5[d5$Es=="GRC_2016"]=300005
d5$CSESid.6[d5$Es=="GRC_2016"]=300006

d5$CSESid.1[d5$Es=="IRL_2016"]=372001
d5$CSESid.2[d5$Es=="IRL_2016"]=372002
d5$CSESid.3[d5$Es=="IRL_2016"]=372003
d5$CSESid.4[d5$Es=="IRL_2016"]=372004
d5$CSESid.5[d5$Es=="IRL_2016"]=372005
d5$CSESid.6[d5$Es=="IRL_2016"]=372006

d5$CSESid.1[d5$Es=="ISL_2016"]=352001
d5$CSESid.2[d5$Es=="ISL_2016"]=352002
d5$CSESid.3[d5$Es=="ISL_2016"]=352003
d5$CSESid.4[d5$Es=="ISL_2016"]=352004
d5$CSESid.5[d5$Es=="ISL_2016"]=352005
d5$CSESid.6[d5$Es=="ISL_2016"]=352006

d5$CSESid.1[d5$Es=="ISL_2017"]=352101
d5$CSESid.2[d5$Es=="ISL_2017"]=352102
d5$CSESid.3[d5$Es=="ISL_2017"]=352103
d5$CSESid.4[d5$Es=="ISL_2017"]=352104
d5$CSESid.5[d5$Es=="ISL_2017"]=352105
d5$CSESid.6[d5$Es=="ISL_2017"]=352106

d5$CSESid.1[d5$Es=="ITA_2018"]=380001
d5$CSESid.2[d5$Es=="ITA_2018"]=380002
d5$CSESid.3[d5$Es=="ITA_2018"]=380003
d5$CSESid.4[d5$Es=="ITA_2018"]=380004
d5$CSESid.5[d5$Es=="ITA_2018"]=380005
d5$CSESid.6[d5$Es=="ITA_2018"]=380006

d5$CSESid.1[d5$Es=="NOR_2017"]=578001
d5$CSESid.2[d5$Es=="NOR_2017"]=578002
d5$CSESid.3[d5$Es=="NOR_2017"]=578003
d5$CSESid.4[d5$Es=="NOR_2017"]=578004
d5$CSESid.5[d5$Es=="NOR_2017"]=578005
d5$CSESid.6[d5$Es=="NOR_2017"]=578006

d5$CSESid.1[d5$Es=="NZL_2017"]=554001
d5$CSESid.2[d5$Es=="NZL_2017"]=554002
d5$CSESid.3[d5$Es=="NZL_2017"]=554003
d5$CSESid.4[d5$Es=="NZL_2017"]=554004
d5$CSESid.5[d5$Es=="NZL_2017"]=554005
d5$CSESid.6[d5$Es=="NZL_2017"]=554006


#d5 party identification (Pi):----

d5$Pi=NA

#assign Pi:
d5$Pi[d5$Es=="AUS_2019"&d5$E3024_3==036001]=1##
d5$Pi[d5$Es=="AUS_2019"&d5$E3024_3==036002]=2
d5$Pi[d5$Es=="AUS_2019"&d5$E3024_3==036003]=3
d5$Pi[d5$Es=="AUS_2019"&d5$E3024_3==036004]=4
d5$Pi[d5$Es=="AUS_2019"&d5$E3024_3==036005]=5
d5$Pi[d5$Es=="AUS_2019"&d5$E3024_3==036006]=6

d5$Pi[d5$Es=="AUT_2017"&d5$E3024_3==040001]=1#
d5$Pi[d5$Es=="AUT_2017"&d5$E3024_3==040002]=2
d5$Pi[d5$Es=="AUT_2017"&d5$E3024_3==040003]=3
d5$Pi[d5$Es=="AUT_2017"&d5$E3024_3==040004]=4
d5$Pi[d5$Es=="AUT_2017"&d5$E3024_3==040005]=5
d5$Pi[d5$Es=="AUT_2017"&d5$E3024_3==040006]=6

d5$Pi[d5$Es=="DEU_2017"&d5$E3024_3==276009]=1#vote for either CDU,CSU (same comments than for 2760001)
d5$Pi[d5$Es=="DEU_2017"&d5$E3024_3==276002]=2
d5$Pi[d5$Es=="DEU_2017"&d5$E3024_3==276003]=3
d5$Pi[d5$Es=="DEU_2017"&d5$E3024_3==276004]=4
d5$Pi[d5$Es=="DEU_2017"&d5$E3024_3==276005]=5
d5$Pi[d5$Es=="DEU_2017"&d5$E3024_3==276006]=6

d5$Pi[d5$Es=="GRC_2016"&d5$E3024_3==300001]=1#
d5$Pi[d5$Es=="GRC_2016"&d5$E3024_3==300002]=2
d5$Pi[d5$Es=="GRC_2016"&d5$E3024_3==300003]=3
d5$Pi[d5$Es=="GRC_2016"&d5$E3024_3==300004]=4
d5$Pi[d5$Es=="GRC_2016"&d5$E3024_3==300005]=5
d5$Pi[d5$Es=="GRC_2016"&d5$E3024_3==300006]=6

d5$Pi[d5$Es=="IRL_2016"&d5$E3024_3==372001]=1#
d5$Pi[d5$Es=="IRL_2016"&d5$E3024_3==372002]=2
d5$Pi[d5$Es=="IRL_2016"&d5$E3024_3==372003]=3
d5$Pi[d5$Es=="IRL_2016"&d5$E3024_3==372004]=4
d5$Pi[d5$Es=="IRL_2016"&d5$E3024_3==372005]=5
d5$Pi[d5$Es=="IRL_2016"&d5$E3024_3==372006]=6

d5$Pi[d5$Es=="ISL_2016"&d5$E3024_3==352001]=1##
d5$Pi[d5$Es=="ISL_2016"&d5$E3024_3==352002]=2
d5$Pi[d5$Es=="ISL_2016"&d5$E3024_3==352003]=3
d5$Pi[d5$Es=="ISL_2016"&d5$E3024_3==352004]=4
d5$Pi[d5$Es=="ISL_2016"&d5$E3024_3==352005]=5
d5$Pi[d5$Es=="ISL_2016"&d5$E3024_3==352006]=6

d5$Pi[d5$Es=="ISL_2017"&d5$E3024_3==352101]=1##
d5$Pi[d5$Es=="ISL_2017"&d5$E3024_3==352102]=2
d5$Pi[d5$Es=="ISL_2017"&d5$E3024_3==352103]=3
d5$Pi[d5$Es=="ISL_2017"&d5$E3024_3==352104]=4
d5$Pi[d5$Es=="ISL_2017"&d5$E3024_3==352105]=5
d5$Pi[d5$Es=="ISL_2017"&d5$E3024_3==352106]=6

d5$Pi[d5$Es=="ITA_2018"&d5$E3024_3==380001]=1#Movimento 5 Stelle
d5$Pi[d5$Es=="ITA_2018"&d5$E3024_3==380002]=2#Partito Democratico (PD)
d5$Pi[d5$Es=="ITA_2018"&d5$E3024_3==380003]=3#Lega (LN)
d5$Pi[d5$Es=="ITA_2018"&d5$E3024_3==380004]=4
d5$Pi[d5$Es=="ITA_2018"&d5$E3024_3==380005]=5
d5$Pi[d5$Es=="ITA_2018"&d5$E3024_3==380006]=6

d5$Pi[d5$Es=="NOR_2017"&d5$E3024_3==578001]=1##
d5$Pi[d5$Es=="NOR_2017"&d5$E3024_3==578002]=2
d5$Pi[d5$Es=="NOR_2017"&d5$E3024_3==578003]=3
d5$Pi[d5$Es=="NOR_2017"&d5$E3024_3==578004]=4
d5$Pi[d5$Es=="NOR_2017"&d5$E3024_3==578005]=5
d5$Pi[d5$Es=="NOR_2017"&d5$E3024_3==578006]=6

d5$Pi[d5$Es=="NZL_2017"&d5$E3024_3==554001]=1##
d5$Pi[d5$Es=="NZL_2017"&d5$E3024_3==554002]=2
d5$Pi[d5$Es=="NZL_2017"&d5$E3024_3==554003]=3
d5$Pi[d5$Es=="NZL_2017"&d5$E3024_3==554004]=4
d5$Pi[d5$Es=="NZL_2017"&d5$E3024_3==554005]=5
d5$Pi[d5$Es=="NZL_2017"&d5$E3024_3==554006]=6


#MERGE DATASETS----

colnames(d1)
colnames(d5)
d=rbind(d1[,c(279,11,6,2,1,280,36,277,278,47,281,85:90,94:99,284,285:290,282,294,295,51,296:301,127:132,136:141,302)]
        ,d5[,c(529,11,5,1,6,530,55,527,528,111,531,136:141,145:150,533,534:539,108,543,544,115,545:550,234:239,243:248,551)])
colnames(d)
rm(d1,d11,d12,d13,d14,d5)


#CMPid:----

#read party Id match CSES-CMP:
dID=read_csv("CSES to CMP by party ID.csv")

dID=dID[,c(3,4,7,8)]#filter data and set variable names:
colnames(dID)[colnames(dID)=="num_id"]="CSESid"#"Vn"
colnames(dID)[colnames(dID)=="manifestoid"]="CMPid"
colnames(dID)[colnames(dID)=="elecyear"]="Ey"
colnames(dID)[colnames(dID)=="country"]="Ec"

#code party family alternative i (CMPid.i):
d=merge(d, dID, by.x=c("Ey","Ec","CSESid.1"), by.y=c("Ey","Ec","CSESid"), all.x=T)#alternative 1
colnames(d)[colnames(d)=="CMPid"]="CMPid.1"
d=merge(d, dID, by.x=c("Ey","Ec","CSESid.2"), by.y=c("Ey","Ec","CSESid"), all.x=T)#alternative 2
colnames(d)[colnames(d)=="CMPid"]="CMPid.2"
d=merge(d, dID, by.x=c("Ey","Ec","CSESid.3"), by.y=c("Ey","Ec","CSESid"), all.x=T)#alternative 3
colnames(d)[colnames(d)=="CMPid"]="CMPid.3"
d=merge(d, dID, by.x=c("Ey","Ec","CSESid.4"), by.y=c("Ey","Ec","CSESid"), all.x=T)#alternative 4
colnames(d)[colnames(d)=="CMPid"]="CMPid.4"
d=merge(d, dID, by.x=c("Ey","Ec","CSESid.5"), by.y=c("Ey","Ec","CSESid"), all.x=T)#alternative 5
colnames(d)[colnames(d)=="CMPid"]="CMPid.5"
d=merge(d, dID, by.x=c("Ey","Ec","CSESid.6"), by.y=c("Ey","Ec","CSESid"), all.x=T)#alternative 6
colnames(d)[colnames(d)=="CMPid"]="CMPid.6"


#CMPrile:----

#read CMP dataset:
dCMP=read.xlsx("MPDataset_MPDS2020b.xlsx")

#filter CMP dataset and set variable names:
colnames(dCMP)[colnames(dCMP)=="party"]="CMPid_temp"
colnames(dCMP)[colnames(dCMP)=="rile"]="CMPrile"

dCMP$Ec=NA
dCMP$Ec=as.character(dCMP$Ec)
dCMP$Ec[dCMP$country==63]="AUS"
dCMP$Ec[dCMP$country==42]="AUT"
dCMP$Ec[dCMP$country==62]="CAN"
dCMP$Ec[dCMP$country==41]="DEU"
dCMP$Ec[dCMP$country==13]="DNK"
dCMP$Ec[dCMP$country==33]="ESP"
dCMP$Ec[dCMP$country==14]="FIN"
dCMP$Ec[dCMP$country==51]="GBR"
dCMP$Ec[dCMP$country==34]="GRC"
dCMP$Ec[dCMP$country==53]="IRL"
dCMP$Ec[dCMP$country==15]="ISL"
dCMP$Ec[dCMP$country==72]="ISR"
dCMP$Ec[dCMP$country==32]="ITA"
dCMP$Ec[dCMP$country==22]="NLD"
dCMP$Ec[dCMP$country==12]="NOR"
dCMP$Ec[dCMP$country==64]="NZL"
dCMP$Ec[dCMP$country==35]="PRT"
dCMP$Ec[dCMP$country==11]="SWE"
dCMP$Ec[dCMP$country==83]="EST"#
dCMP$Ec[dCMP$country==86]="HUN"
dCMP$Ec[dCMP$country==88]="LTU"
dCMP$Ec[dCMP$country==97]="SVN"

#GRC_2012: there was an election one month before in May, need to exclude it from dCMP
# to merge in dt only the parties in the June election (which is the one coded by CSES):
dCMP=dCMP[!(dCMP$Ec=="GRC"&dCMP$date==201205),]

dCMP$Ey=as.numeric(substr(dCMP$date,1,4))
dCMP$Ey[dCMP$Ec=="GRC"&dCMP$date==201509]=2016#set temporary Ey for double elections

dCMP$Es=paste0(dCMP$Ec,"_",dCMP$Ey)

#keep only variables of interest:
dCMP=dCMP[,c(175,176,7,168)]

#code pervote alternative i (CMPrile.i):
d=merge(d, dCMP, by.x=c("Ec","Ey","CMPid.1"), by.y=c("Ec","Ey","CMPid_temp"), all.x=T)
colnames(d)[colnames(d)=="CMPrile"]="CMPrile.1"
d=merge(d, dCMP, by.x=c("Ec","Ey","CMPid.2"), by.y=c("Ec","Ey","CMPid_temp"), all.x=T)
colnames(d)[colnames(d)=="CMPrile"]="CMPrile.2"
d=merge(d, dCMP, by.x=c("Ec","Ey","CMPid.3"), by.y=c("Ec","Ey","CMPid_temp"), all.x=T)
colnames(d)[colnames(d)=="CMPrile"]="CMPrile.3"
d=merge(d, dCMP, by.x=c("Ec","Ey","CMPid.4"), by.y=c("Ec","Ey","CMPid_temp"), all.x=T)
colnames(d)[colnames(d)=="CMPrile"]="CMPrile.4"
d=merge(d, dCMP, by.x=c("Ec","Ey","CMPid.5"), by.y=c("Ec","Ey","CMPid_temp"), all.x=T)
colnames(d)[colnames(d)=="CMPrile"]="CMPrile.5"
d=merge(d, dCMP, by.x=c("Ec","Ey","CMPid.6"), by.y=c("Ec","Ey","CMPid_temp"), all.x=T)
colnames(d)[colnames(d)=="CMPrile"]="CMPrile.6"


#CMPper415:----

#read CMP dataset:
dCMP=read.xlsx("MPDataset_MPDS2020b.xlsx")

#filter CMP dataset and set variable names:
colnames(dCMP)[colnames(dCMP)=="party"]="CMPid_temp"
colnames(dCMP)[colnames(dCMP)=="per415"]="CMPper415"

dCMP$Ec=NA
dCMP$Ec=as.character(dCMP$Ec)
dCMP$Ec[dCMP$country==63]="AUS"
dCMP$Ec[dCMP$country==42]="AUT"
dCMP$Ec[dCMP$country==62]="CAN"
dCMP$Ec[dCMP$country==41]="DEU"
dCMP$Ec[dCMP$country==13]="DNK"
dCMP$Ec[dCMP$country==33]="ESP"
dCMP$Ec[dCMP$country==14]="FIN"
dCMP$Ec[dCMP$country==51]="GBR"
dCMP$Ec[dCMP$country==34]="GRC"
dCMP$Ec[dCMP$country==53]="IRL"
dCMP$Ec[dCMP$country==15]="ISL"
dCMP$Ec[dCMP$country==72]="ISR"
dCMP$Ec[dCMP$country==32]="ITA"
dCMP$Ec[dCMP$country==22]="NLD"
dCMP$Ec[dCMP$country==12]="NOR"
dCMP$Ec[dCMP$country==64]="NZL"
dCMP$Ec[dCMP$country==35]="PRT"
dCMP$Ec[dCMP$country==11]="SWE"
dCMP$Ec[dCMP$country==83]="EST"#
dCMP$Ec[dCMP$country==86]="HUN"
dCMP$Ec[dCMP$country==88]="LTU"
dCMP$Ec[dCMP$country==97]="SVN"

#GRC_2012: there was an election one month before in May, need to exclude it from dCMP
# to merge in dt only the parties in the June election (which is the one coded by CSES):
dCMP=dCMP[!(dCMP$Ec=="GRC"&dCMP$date==201205),]

dCMP$Ey=as.numeric(substr(dCMP$date,1,4))
dCMP$Ey[dCMP$Ec=="GRC"&dCMP$date==201509]=2016#set temporary Ey for double elections

dCMP$Es=paste0(dCMP$Ec,"_",dCMP$Ey)

#keep only variables of interest:
dCMP=dCMP[,c(175,176,7,59)]

#code pervote alternative i (CMPper415.i):
d=merge(d, dCMP, by.x=c("Ec","Ey","CMPid.1"), by.y=c("Ec","Ey","CMPid_temp"), all.x=T)
colnames(d)[colnames(d)=="CMPper415"]="CMPper415.1"
d=merge(d, dCMP, by.x=c("Ec","Ey","CMPid.2"), by.y=c("Ec","Ey","CMPid_temp"), all.x=T)
colnames(d)[colnames(d)=="CMPper415"]="CMPper415.2"
d=merge(d, dCMP, by.x=c("Ec","Ey","CMPid.3"), by.y=c("Ec","Ey","CMPid_temp"), all.x=T)
colnames(d)[colnames(d)=="CMPper415"]="CMPper415.3"
d=merge(d, dCMP, by.x=c("Ec","Ey","CMPid.4"), by.y=c("Ec","Ey","CMPid_temp"), all.x=T)
colnames(d)[colnames(d)=="CMPper415"]="CMPper415.4"
d=merge(d, dCMP, by.x=c("Ec","Ey","CMPid.5"), by.y=c("Ec","Ey","CMPid_temp"), all.x=T)
colnames(d)[colnames(d)=="CMPper415"]="CMPper415.5"
d=merge(d, dCMP, by.x=c("Ec","Ey","CMPid.6"), by.y=c("Ec","Ey","CMPid_temp"), all.x=T)
colnames(d)[colnames(d)=="CMPper415"]="CMPper415.6"


#VPf:----

#read CMP data and party Id match CSES-CMP:
dID=read_csv("CSES to CMP by party ID.csv")
dCMP=read.xlsx("MPDataset_MPDS2020b.xlsx")

dID=dID[,c(7,8)]#filter data and set variable names:
colnames(dID)[colnames(dID)=="num_id"]="CSESid"#"Vn"
colnames(dID)[colnames(dID)=="manifestoid"]="CMPid"
dID=dID[!duplicated(dID$CSESid),]

dCMP=dCMP[,c(7,10)]#filter data and set variable names:
colnames(dCMP)[colnames(dCMP)=="party"]="CMPid"
colnames(dCMP)[colnames(dCMP)=="parfam"]="Pf"
dCMP=dCMP[!duplicated(dCMP$CMPid),]

#code party family vote choice (VPf):
d=merge(d, dID, by.x=c("Vn"), by.y=c("CSESid"), all.x=T)
d=merge(d, dCMP, by = c("CMPid"), all.x=T)
# sort(unique(d__$Vn[is.na(d__$Pf)]))#the missing Pf are in general for voters for parties that are not the top 6 and hence will be dropped out later on
colnames(d)[colnames(d)=="Pf"]="VPf"
d=d[,-which(names(d) %in% c("CMPid"))]

rm(dID,dCMP)


#PrePost:----
d$PrePost=as.numeric(1)

# 1. Post-election study
# 2. Pre-election study
# 3. Panel, with some pre-election and some post-election data


#Ed election dataset:----
d$Ed="CSES"


#SET DOUBLE ELECTIONS BACK TO ORIGINAL VALUES:----

#settign GRC 2015 elections back to correct year and election study:
d$Es[d$Ec=="GRC"&d$Ey==2015]="GRC_2015j"#January election
d$Es[d$Ec=="GRC"&d$Ey==2016]="GRC_2015s"#September election
d$Ey[d$Ec=="GRC"&d$Ey==2016]=2015


#SAVE DATA:----
colnames(d)
d=d[,c(2:3,16:24,1,25:46,72:73,47:58,60:65,15:4,66:71,59)]
save(d, file = "d29_CSES.RData")

